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Abstract 

In this thesis, we consider the massive field theories in 1+1 dimensions known as 
affine Toda quantum field theories. These have the special property that they possess 
an infinite number of conserved quantities, a feature which greatly simplifies their 
study, and makes extracting exact information about them a tractable problem. 
We consider these theories both in the full space (the bulk) and in the half space 
bounded by an impenetrable boundary at x = 0. In particular, we consider their 
fundamental objects: the scattering matrices in the bulk, and the reflection factors 
at the boundary, both of which can be found in a closed form. 

In Chapter 1, we provide a general introduction to the topic before going on, in 
Chapter 2, to consider the simplest ATFT — the sine-Gordon model — with a bound- 
ary. We begin by studying the classical limit, finding quite a clear picture of the 
boundary structure we can expect in the quantum case, which is introduced in Chap- 
ter 3. We obtain the bound-state structure for all integrable boundary conditions, 
as well as the corresponding reflection factors. This structure turns out to be much 
richer than had hitherto been imagined. 

We then consider more general ATFTs in the bulk. The sine-Gordon model is 
based on a± , but there is an ATFT for any semi-simple Lie algebra. This under- 
lying structure is known to show up in their S-matrices, but the path back to the 
parameters in the Lagrangian is still unclear. We investigate this, our main result 
being the discovery of a "generalised bootstrap" equation which explicitly encodes 
the Lie algebra into the S-matrix. This leads to a number of new S-matrix identities, 
as well as a generalisation of the idea that the conserved charges of the theory form 
an eigenvector of the Cartan matrix. 

Finally, our results are summarised in Chapter 5, and possible directions for 
further study are highlighted. 



Note added 



The work which follows was presented as a thesis for the degree of Doctor of Philosophy at the 
University of Durham in September 2000. Since then, Bajnok et al. [|62], 



64] have also studied 



the problem of the boundary sine-Gordon model with arbitrary integrable boundary conditions. 
Their results provide independent confirmation of many of those reported here. We are also 
grateful to Gabor Takacs and Zoli Bajnok, in particular, for pointing out a number of typos 
Since they apply equally to the content of this thesis, we record them here. 



m 



• We should perhaps have made it clear that the lemmas presented in Chapter 3 rested on an 
implicit assumption that, when the the incoming particle decays into two other particles, 
neither is the outgoing particle. Examples such as Figure 3.7 are therefore technically 
exceptions to the lemmas as stated. These cases were not discussed explicitly in the text 
as they are only relevant when the rapidity of the incoming particle is independent of the 
boundary parameter. 

• In equation (3.24) [equation (18) in Q], the two 1/2X terms should read //A. 

• Just above Table 3.1 [Table 1 in (lj], 1 should read a\' m . 

• In Table 3.7 [Table 7 in fij], the order of the pole at uq is not correct; it should be k. 



The hep-th version of (T| has been updated to incorporate the last three amendments. 
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Chapter 1 



Integrable Quantum Field Theory in 
Two Dimensions 

"/ got really fascinated by these 1+1 models . . . and how mysteriously they jump out 
at you and work and you don 't know why. I am trying to understand all this better. " 

— Richard Feynman 

1.1 Introduction 

'"The time has come, ' the walrus said, 'to talk of many things. Of ships and shoes and 
sealing-wax, and cabbages and kings. 

— Lewis Carroll 

Why study a theory in two dimensions, when the real world has at least four? The difficulty at 
least at present, is that realistic four-dimensional field theories are incredibly complicated, even 
before the further dimensions proposed by superstring theories are considered. Perturbative 
solutions can be found, but exact, non-perturbative, results are few. In a number of cases, 
sufficiently accurate perturbative answers are enough, but many physical phenomena cannot be 
properly understood this way. Uncovering exact solutions to quantum field theories is considered 
by many to be one of the great remaining challenges to particle theorists. 

In view of this, it is perhaps logical to think of taking a step backwards, to consider a simpler 
model which still exhibits some of the same features. Understanding the issues involved a few 
at a time in this way presents a more manageable problem, like climbing a mountain in stages 
rather than hoping to stride straight to the top. A theory with two dimensions, one of space 
and one of time, is a useful starting point. 

The main focus of this thesis will be the scattering matrices (S-matrices) of these theories, 
which describe how the final state of the system is related to the initial state. In general this is 
a messy object to deal with, since there are potentially matrices for the evolution of any number 
of particles into any other number. However, the theories considered here are "simplified" in 
one further respect in that they are integrable theories. This has three main effects: 

• there is no net particle production, so the number of particles in the initial and final states 
must always be the same; 

• the outgoing particles must have the same masses and velocities as the incoming ones; 
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• a general S-matrix, for the scattering of n particles, can always be factorised into a product 
of two-particle S-matrices. 

In essence, this means that, once the S-matrix for the scattering of two particles into two particles 
has been calculated, everything else follows, making a characterisation of the theory in terms of 
these matrices a more attractive and tractable proposition. 

The method for obtaining explicit expressions for these S-matrices is in fact surprisingly 
straightforward (if not necessarily simple to put into practice). It invokes four consistency 
requirements which, between them, provide strong constraints on the S-matrix. This method, 
first formulated in the 1960s, was initially developed to help explain the strong nuclear force 
(see e.g. [||, [|). In the 1970s, it was discovered that these axioms could also be applied to some 
two dimensional quantum field theories, and proved in certain cases to be powerful enough to 
completely determine the S-matrix up to an overall factor. 

A theory is integrable if it has an infinite number of symmetries; the particular theories we 
will be studying, the affine Toda field theories (ATFTs), acquire these through being based on an 
infinite-dimensional Lie algebra. We will study these both in the full 2-dimensional space, and 
in a half space (or half line) defined by introducing an impenetrable boundary at x = 0. As well 
as the usual S-matrices, this requires the introduction of boundary factors to describe particles 
scattering against the boundary. The particles can either simply reflect from this "wall" or bind 
to it, and we will be concerned with the bound state structure this introduces. 

The situation with a boundary is the less well-understood of the two, so we shall study only 
the simplest ATFT, the sine-Gordon model. Even for this case, only the ground and first few 
excited states have been explored. We present a complete description of these matrices, for any 
wall which leaves the resultant theory integrable. 

Without a boundary, the picture is much clearer, and S-matrices have been found for all real- 
coupling ATFTs. However, despite the manifest Lie algebraic structure of the theory, its path 
from the Lagrangian to the S-matrix is still not precisely known, and remains an open problem. 
We will, however, present a convenient way of encoding the Lie algebra into the matrix, with 
the aim of making the task a little easier. This process will also throw up a number of new 
relationships between elements of the S-matrix. 

Apart from their interest in connection with higher-dimensional theories, two-dimensional 
integrable models have an increasing number of applications in their own right. They are, for 
example, useful in studying impurity problems in an interacting ID electron gas Q or edge 
excitations in fractional quantum Hall states [Q, H|. (A recent review can be found in J10||.) 

In the following sections all this will be put on a more formal basis, paving the way for the 
discussion of the ATFTs which will occupy our attention for the remainder of the thesis. 

1.2 Exact S-matrices 

"The three rules of the Librarians of Time and Space are: 1) Silence; 2) Books must be 
returned no later than the date last shown; and 
3) Do not interfere with the nature of causality." 

— Terry Pratchett, Guards! Guards! 



Much of the discussion in this section is based on [11, 12, 13]. Before proceeding to specifics, 
a gentle introduction to two-dimensional field theory is perhaps appropriate. 
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Let us begin by considering a general Euclidean field theory with one space and one time 
dimension (x l ,x 2 ) = (x,t) defined (in the Lagrangian approach) by the classical action 

/+oo r+oo 
dx / dta((p,d^ip), (1.1) 
-oo J — oo 

where tp(x, t) is some set of fundamental fields and the action density a(ip, d^ip) is a local function 
of these fields and the derivatives d^ip = d(p/dx^ with \i = 1,2. For simplicity we shall also 
use light-cone coordinates, so that, in place of (p ^ 1 ) for the two-momentum, we will take 
(p,p) = (p° + p 1 ,^ — p l ). 

We will be considering the particles only on mass-shell (i.e. real, rather than virtual, parti- 
cles), which means that their two-momenta satisfy the mass-shell condition 

PaPa = (Pa) 2 -(p 1 a) 2 = m 2 a , (a = 1, 2, . . . ,7l). (1.2) 

The two momenta can be conveniently parametrised in terms of their rapidity 9: 

(p° a ,p 1 a ) = (m a e e «,m a e- e «). (1.3) 

Suppose our theory contains n different types of particle A a ,a = 1,2,... , n, with masses 
m a . The asymptotic particle states are generated by the "particle creation operators" A a (6): 

\A ai {9 1 )A a2 {9 2 ) • • • A an {e n )) = A ai (0i)4«(02) • • • A an (9 n )\0). (1.4) 

Looking into the far past, we shall call the state an in state if there are no further interactions 
as t — > — oo. This means that the fastest particle must be on the left, and the slowest on the 
right, with all the others in order in-between, i.e. 9\ > 62 > ■ ■ ■ > n . Similarly, if there are no 
further interactions as t — > 00, the state will be called an out state, and the rapidities must be 
in the reverse order. 

The S-matrix can now be introduced as a mapping between the in-state basis and the out- 
state basis. It is useful to consider the A x (9)s as non-commuting symbols, giving them an exis- 
tence outside the ket vectors, so that we can write the above state simply as A ai (9i)A a2 (92) • • • 
A an (9 n ). Considering an m-particle in-state, we have 

A ai {9 1 )A a2 {9 2 )...A an {9 n ) = 

00 

E E ^::i:( i---^^i---C)Ai(^M6 2 (^)...^ m (C) (i-5) 

m=18' 1 <...<e' m 

where a sum on b± ... b n is implied, and the sum on the 9[ will, in general, turn out to involve a 
number of integrals. The rapidities will also be constrained by momentum conservation. 

For a general theory, we can proceed no further, and introducing the S-matrix would appear 
only to have complicated matters. However, for an integrable theory, the whole situation becomes 
dramatically simpler. The name derives from the classical formulation of such theories, which 
can be cast as partial differential equations; these were said to be integrable if it was possible 
to find an explicit solution. It was found that a solution was only possible if there were an 
infinite number of symmetries constraining the behaviour of the equation, and preventing it 
from becoming chaotic. The same applies here: possessing so many symmetries constrains the 
S-matrix sufficiently to allow an exact solution to be found. 

Energy- momentum is always a conserved quantity, and its operator, P, is said to have 
(Lorentz) spin 1 as it transforms under a Lorentz boost L a : 9 — > 9' = 9 + a as P — > P' = e a P. 
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This means that a boost of 2in — a complete rotation — has no effect. On a one-particle state, 
the action of P = (P, P) would be 

P\A a {9)) = m a e e \A a (8)), P\A a (9)) = m a e~ e \A a (9)) . (1.6) 

In general, there can also be other conserved quantities, Q s , which transform in higher represen- 
tations of the 1+1-dimensional Lorentz group as Q s — > Q^, = e as Q s and have spin s since they 
rotate s times under a boost of 2in. This time, the effect of Q s = (Q s ,Q-s) on a one-particle 
state is 

Q s \A a (6)) = q^e se \A a (9)}, Q- s \A a (6)} = q^e^A^)) . (1.7) 



In an integrable theory, there are an infinite number of these conserved quantities (or 
"charges"). It might, at first, appear that such theories are quite improbable. In the theo- 
ries to be considered later, these symmetries are due to an underlying group structure which 
happens to be infinite-dimensional. 

We will concentrate on local conserved charges, which are those whose operators are integrals 
of strictly local densities, meaning that their action on multi-particle states is additive: 

QslAaM . . . A an (e n )) = (?i?e s01 + . . . + q^e^lAaM . . . A an (6 n )). (1.8) 



Just as, above, momentum conservation constrained the sum over the rapidities, so all the 
other conserved quantities provide additional constraints, leaving an infinite number of equations 
to be solved, of the form 

flS?** + ■■■ + = + . . . + (L9) 

Since these must all hold for all possible sets of in-momenta, the only possible solution is the 
trivial one, i.e. n = m, and Oi = 9[, qif = q^ for all i. 

This establishes the fact that there can be no particle production in an integrable theory, 
and that the sets of incoming and outgoing momenta must be equal, thus reducing the workload 
involved in dealing with the S-matrix just to the n — ► n cases. There is, however, one further 
property of integrable theories which makes them even easier to deal with: factorisation. This 
states that, for any n — > n S-matrix, the trajectories of the particles involved can be shifted 
forwards or backwards in space so as to split the vertex into a product of ^n(n — 1) two-particle 



vertices, as shown in figure 1.2. 



The origin of this property can be seen in the fact that, while the momentum operator, for 
example, will act on a given state simply by shifting the position of all the particles by a fixed 
amount, higher-spin operators will, in general, change the positions by an amount depending 
on the initial momentum of the particle. This argument was first proposed by Shankar and 
Witten [14] and elaborated by Parke |]l^]. In rough form it goes as follows. 

The first step is to note that, since we are dealing with a local, causal field theory, the 
particles in any process are sufficiently well separated (at least most of the time) to be consid- 
ered individually, and it is reasonable to consider the effect of the conserved charges particle 
by particle. If we consider a single-particle state, with position approximately x\ and spatial 
momentum approximately pi, the position space wavefunction will be 

i>(x)<x / dpe~ a2(p ~ pl)2 e ip(x - Xl) . (1.10) 
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For simplicity, rather than considering a general spin-s operator, we will now try acting on 
this with P s , the spin s operator which acts as (P) s , i.e. as s copies of the spatial part of the 
two-momentum operator. Applying e~ iaPs to the above wavefunction gives 



/+oo 
dpe -a a ( P -pi) a e ip(x-xi) e -iap« i (1.11) 
-oo 



Since most of the value of the integral is due to the region around p « p%, we can Taylor expand 
the extra factor in powers of (p — Pi) to find new values for the position and momentum. For 
a general momentum-dependent phase factor e~ l ^ p \ this leaves the momentum unchanged but 
shifts the position by x\ = x\+4>'(pi). Here, this gives a position shift of sap^ -1 . For momentum 
itself, this is just a but, for higher spins, the shift must depend on the initial momentum. This, 
as Parke showed, is a general property of higher-spin operators. 

Applying a suitable operator, Q s , near an n — > n vertex thus separates the particles and 
splits up the vertex into \n(n — 1) 2 — > 2 vertices. However, since the operator is related to a 
conserved charge, the amplitude for both processes must be the same, giving us factorisability. 



In addition, applying Q- s rather than Q s causes a mirror-image split, leading to figure 1.2 



This relies, of course, on the fact that, after splitting up the vertex, the particle trajectories 
must still cross somewhere, because we are restricted to only two dimensions. In higher numbers 
of dimensions, it is quite easy to imagine splitting up such a vertex so that the particles never 
meet at all, leading to a trivial S-matrix. In fact, Coleman and Mandula proved the so-called 
Coleman-Mandula theorem, showing that, for any theory with more than one space dimension 
and a conserved charge with spin 2 or more, the S-matrix must be trivial. 

All this shows that an integrable theory in 1+1 dimensions is rather special, in that it has: 



• no particle production; 

• equality of the sets of initial and final momenta; 

• factorisability of the n — > n S-matrix into a product of 2 — > 2 S-matrices. 



With these results, it is clear that the 2 — * 2 S-matrix is the fundamental object of the 
theory, and that, once it has been found, the full S-matrix is only a step away. The 2 — ► 2 
process is just 

A ai {ei)A a2 {9 2 ) = - 6 2 )A b2 (6 2 )A bl (6 1 ), (1.12) 

and is shown graphically in figure |1.1| . (In this, and in all subsequent diagrams, time is taken 
to run up the page, and space from left to right.) Note that momentum conservation demands 
m ai = and m a2 = rrib 2 , so that a\ ^ b% or a 2 ^ b 2 are only possible if there is a degenerate 
mass spectrum. 

The two-particle S-matrix has n 4 elements, but these are not all independent, and are in fact 
strongly constrained. Firstly, it is generally assumed that parity, charge conjugation and time 
reversal (P, C and T) symmetries hold. These impose the conditions 

SV(9) = S$(e)=S%(9) = S%(9). (1.13) 

In addition, they must satisfy four general axioms: the Yang-Baxter equation; a unitarity con- 
dition; analyticity and crossing symmetry; and the bootstrap condition. These are powerful 
demands; using just the first three allows the S-matrix to be pinned down up to the so-called 
"CDD ambiguity" : 

Sg(0)^Sg(0)$(0), (1-14) 
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02 h 



Figure 1.1: S-matrix 



a l CL2 (Z3 dl CL2 a% 




Figure 1.2: Yang-Baxter equation 

where the "CDD factor" satisfies 

$(0) = $(t7r - 0) , $(0)$(-0) = i, (i.i5) 

but is otherwise arbitrary. This can often be further restricted by the bootstrap. 
1.2.1 Yang-Baxter (or "factorisation") equation 

The requirement of factorisability, and, in particular, the ability of operators associated to 
conserved charges to shift trajectories around, is only consistent if figure is true. This gives 
rise to the condition 

S^%(0)^(6 + e')S h c f c l{9') = S:iZ(0')S c a \%(9 + 9>)S b c ^(9). (1.16) 

Formally, this is an associativity condition on the algebra of the Ai(9)s: moving from an 
m-state A ai (9i)A a2 {92) ■ ■ ■ A an {9 n ) to an out-state by a series of pair transpositions, the result 
is independent of their order if and only if the Yang-Baxter equation holds. For three particles, 
there are only two ways of doing this (shown in figure |1.2[ ). For the left-hand diagram, we find 

A^OA^Ks^) = [szzmcMAcMjA^) 

= Sl\% (9) Ac (0 2 )[S b c l C al (0 + 0')4 3 (0 3 )A bl (9 1 )] (1.17) 

= S£g(0)S§2(0 + , )fc b 3 3 (0 , )4 3 (0 3 )4 2 (02)]4 1 (0i) , 
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Physical values 

D 

m xxxxxxxx • 
A {mi — rrij) 2 (nii + rrij) 2 

Cut 



C 




Cut 



Figure 1.3: The complex s-plane 



where we have set 9 = 62 — 0\ and 9' = 9 3 — 92- Doing the same for the right-hand diagram 
yields a relation between the same in- and out-states with a different product of S-matrices. 
Since these equations should be equivalent, the products of S-matrices can be equated, to give 
( |1.16l ). That no further conditions should arise in considering larger numbers of particles can 
be seen through the fact that the Yang-Baxter equation allows any trajectory to be moved past 
any given vertex (by considering just the local area around the vertex). Thus, with repeated 
applications, trajectories can be moved arbitrarily, showing that all possible factorisations are 
equivalent. 

That the Yang-Baxter equation is an associativity condition can most easily be seen when we 
have a non-degenerate mass spectrum. In this case, we can define operators O a b which transpose 
the symbols A a and Ab, and add a suitable S-matrix factor. The Yang-Baxter equation then 
becomes just 

12 (0 13 23 ) = 23 (0 13 12 ), (1.18) 

which is indeed an associativity condition. The Yang-Baxter equation is the extension of this to 
a degenerate spectrum. 



1.2.2 Unitarity, analyticity and crossing symmetry 

The origin of these demands can best be seen by switching to Mandelstam variables, 

s = (pi +P2) 2 , t = { Pl -p 3 ) 2 , u={pi-p A f (1.19) 

with s + t + u = m i- H ere > pi and P2 are the momenta of the incoming particles, with p 3 
and P4 those of the outgoers. Only one of these is independent, so we shall make the standard 



choice and consider s. Making use of (1.3), this can be re- written as 



+ m 2 + 2mimj cosh(6q - 9 2 ). (1.20) 



For a real, physical process, all rapidities are real, so s must be real and satisfy s > (mi+rrij) 2 



However, it is usual to assume that the S-matrix S(s) is an analytic function]^, and can so be 
continued into the complex plane to be single- valued, at least after suitable cuts have been made. 



As it turns out, this can be achieved with two cuts, as shown in figure 1.3 



The cut plane is the physical sheet of the Riemann surface for S; continuing through one 
of the cuts leads to one of the other, unphysical, sheets. Making the cuts in this way, S is 
single- valued, meromorphic and real-analytic^]. Note also that S(s) is real on the axis between 
the cuts, i.e. for (m; — rrij) 2 < s < (m^ + rrij) 2 . 



x It has been suggested M that this is connected to the causality of the theory. 
2 S takes complex-conjugate values at complex conjugate points. 
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Figure 1.4: Crossing 



Unitarity demands that S(s)S^(s) = 1 for physical values of s (just above the right-hand cut). 
This is a matrix equation, so there is an implicit sum over a complete set of asymptotic states 
living between S and Generally, as s increases, states involving more and more particles 
become available, bringing the 2 — > n /S-matrix into play, for n = 3, 4, . . . Here, however, there 
is no particle production, so this cannot happen and we are left with 



)OTV)]* = W 



:i.2i) 



for all physical s + , with * denoting the complex conjugate. Considering s 
place it just above the cut, real analyticity allows this to be re- written as 



as s + ie (e — ► 0) to 



rn rm 



n.22) 



with s = s — ie, just below the cut. (We have skipped many of the details in the interests of 
simplicity. For a more rigorous explanation, see [11 1 or Jl2[.) 

The other important constraint comes from the fundamentally relativistic property of cross- 
ing. If the interaction is assumed to take place at t = 0, "crossing" one of the participating 
particles involves inverting its path in time, so that incoming particles become outgoing and 
vice versa. In general, if one of the incoming particles to an interaction is crossed to become 
outgoing while one of the outgoers is simultaneously crossed to become incoming, the amplitude 
for another physical process is obtained. 



In our case, this amounts to saying that we can look at figure 1.1 from the side, with the 



forward momentum taken as t rather than s. Normally, t 
can be written as 



(Pi - P3) but, here, p 2 = P3, so it 



t= (Pl~ P2f 



2p 2 + 2p 2 2 



(Pi +P2] 



2m? + 2mh 

6 J 



:i.23) 



The amplitude for this process can be found by analytically continuing from the original am- 



plitude to a region where t is physical, i.e. t is real and t > (rrtj + rrij) 2 . From ( 1.23j ), this 
corresponds to s < (mj — rrij) 2 . Physical amplitudes come from approaching this from above in 



t and hence from below in s; a suitable path for continuation is shown in figure 1.4. As a result, 
we have 



S!j(s+) = sf(2rn 2 + 2m 2 j 



I- 



[1.24) 



This picture becomes substantially simpler if we shift back to the rapidity difference 9 through 
the transformation 



cosh 
log 



_ t | s - mf - m 2 
2mii7ij 



2mirrij 



(1.25) 



s — rrii 



m 2 + \J {s - (mi + nij) 2 }{s - (mi - rrij) 2 } 
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Figure 1.5: The 9 plane 



This maps the physical sheet to the "physical strip" < Im 9 < tt, with the unphysical sheets 
being mapped onto the unphysical strips nn < Im 9 < (n + 1)9. Also, the two branch points go 
to and in, with the cuts opening up as shown in figure [L~5| . 

Analytically continuing to the entire plane, and re- writing in terms of 9, the demands of 
analyticity and crossing symmetry become 

S%»(e)S%%(-6) = % 1 tf* (1.26) 

and 

Saltm = S^liin - 9) (1.27) 

respectively. (Note that, for physical 9, [S{6)\* = S{—9).) These results can be combined into 
the "cross-unitarity equation" 

S***(i* - 0)Sg*(i* + 9) = 5 b a \5 b a l . (1.28) 



This unitarity result can also be understood in terms of the algebra of the A symbols. Since 
we are assuming that the S-matrix is analytic, and so can be defined for all complex 9, it seems 
reasonable to demand that ( |1.12| ) still makes sense if we interchange 9\ and 02- The equation 
now relates in-states to out-states, rather than the other way round. If the original equation is 
then applied to what is now an in-state on the rhs, we find 

AMAaM = £ S££(0i - 9 2 )S c b l c b l(9 2 - 9 1 )A Cl (9 1 )A C2 (e 2 ), (1.29) 

which relates an out-state to a sum of other out-states. However, the out-states form an asymp- 
totically complete basis (as do the m-states) and so cannot be broken down, leading us to 
identify the states on each side of the equation and thus yielding ( 1.26j ). 



If the time and space dimensions could be treated on an equal footing (e.g. by working 
in Euclidean rather than Minkowski space) the crossing symmetry result would have become 
Sa\ b a 2 (9) = S*' 2 '?- 1 (ir — 9) , making it clear that it amounted to allowing figure |El] just to be rotated 
on the page. In Minkowski space, this is still true; rotating the diagram is just not as trivial an 
operation. 
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Figure 1.6: Bound state formation 



1.2.3 The Bootstrap Principle 

It is normally assumed that at least some of the simple poles on the physical strip indicate the 



presence of bound states, either in the forward (s) or crossed (t) channel, as shown in figure 1.6. 
Note that this is consistent with there being no particle production provided such poles do not 
appear for physical values of 6. In fact, poles corresponding to bound states only appear for 
purely imaginary 9, with resonance states possible at complex 8. Note also that simple poles 
do not need to correspond to bound states, a fact that will become important later and will be 



discussed in Section 1.2.4 



There are various reasons why this is taken to be true, such as: 

• in quantum mechanics, if there is a pole in the S-matrix for scattering a particle off a 
potential then the wavefunction for the particle bound to the potential can be constructed; 

• tree- level Feynman diagrams. 

In many other ways, however, it has to be taken as an axiom, without a rigorous basis. 



The "fusing angle" for ij — ► k is denoted as U& (as shown in figure 1.6) and indicates 



that Si? will have a simple pole at iUjj for the forward (s-channel) process, and ir — iUfij for 
the crossed (t-channel) version. The intermediate particle, k, is on-shell and so survives for a 
macroscopic length of time. The "bootstrap principle" (or "nuclear democracy") then states 
that k should be expected to be one of the other asymptotic one-particle states of the model. 

This has proved to be immensely useful in discovering the full structure of models once at 
least the fundamental particles — those from which all other particles can be built up as bound 
states — are known. By looking at the interactions of all known particles, adding any new states 
that show up as bound states of the known ones, and repeating the process until everything can 
be accounted for, all the particles in the theory can be found. That is not to say that the problem 
becomes trivial — the fundamental particles must still be discovered by other means — but it is 
simplified greatly. 

Of course, it is not enough just to discover the presence of new states; we need to know their 



Such poles are always simple, though this is not necessarily the case in field theory. 
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Figure 1.7: The mass triangle 



properties as well, such as their mass and, in particular, their S-matrices with other particles. 
Indeed, it is only through poles in these S-matrices that further new bound states can appear. 



For the forward-channel process, as particle k is on-shell, s 



mi, so we have 



ml 



mf + rrij + 2rriimj cos Uh. 



[1.30) 



This is a well-known trigonometric formula, and implies that can be represented as the 



exterior angle of a triangle of sides rrii , rrij and , as shown in figure |L7j . This also shows that 
the three fusing angles satisfy 



:i.31) 



as might be expected from looking at figure 1.6. 

In addition, extending the dictates of factorisation (i.e. having to allow trajectories to be 



moved past a vertex) to the case where a bound state is formed yields figure 1.8, and the 
corresponding "bootstrap equation" 



J a\a2 caz 



(o) = f b c i C2 s c a \%(e + i u a a %)s: 



MC3 , 



IU 



), 



(1.32) 



where f£ b is the "three-particle coupling", as shown in figure |L9 , At the pole where a bound 
state is formed, the S-matrix can be considered as a pair of such couplings, giving 



sL 3 («0 



fk fi'f 
JijJ k 



iu: 



(1.33) 



This is a great help to the aspiring state-hunter, as treating all the relevant S-matrix elements 
at the point where the new bound state is expected to be formed as a set of simultaneous 
equations for the fs allows them to be found and substituted into ( |l,32j ) to give the S-matrices 
involving the new state. 

Another useful relation comes from equating the action of one of the conserved charges, Q s 
on the state before fusing — 1^4^ (6*i )^4j (6*2)) — and after — \A^(9s)). The action is given by (|1.8|) 
and leads to the "conserved charge bootstrap" 



>) 
k 



■isU 



(1.34) 



It is interesting to note, as was pointed out in [17], that if we take the logarithmic derivative 
of the S-matrix, 



fab{0) = -i—lnS ab (6), 



[1.35) 
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Figure 1.9: Three-particle coupling 



expanded according to 



fc=i 

and insert it into the logarithmic derivative of the bootstrap equation, we recover 

(s) (s) —isJI b (s) i<v~ 



(1.36) 



[1.37) 



showing that the rows and columns of <pW provide solutions for the conserved charge bootstrap 



1.2.4 The Coleman-Thun mechanism 



If all poles were simple, and inevitably corresponded to the creation of a bound state, as in 
quantum mechanics, the story would now be complete. However, this is not the case; not only 
do some theories give rise to double, triple, or higher order poles, but not all simple poles have 
a natural interpretation in terms of bound states. 

The solution to this problem was discovered by Coleman and Thun in 1978 [18], in terms 
of anomalous threshold singularities. For a given Feynman diagram, if the external momenta 
are such that one or more of the internal propagators are simultaneously on-shell (i.e. can be 
considered as real, rather than virtual, particles) then it turns out that the loop integrals give 
rise to a singularity in the amplitude. The bound states considered above are simple examples 
of this, with one propagator (the bound state particle) on-shell. 
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In three or more dimensions, all the singularities which do not correspond to bound states 
show up as branch points, but in 1+1-dimensions, they can appear as poles. The practical upshot 
of this is that such poles should be considered as being due not to the tree-level diagrams we 
have been looking at so far, but to more complicated diagrams which are nonetheless composed 



entirely of on-shell particles, such as the one shown in figure 1.10. This diagram, if it was 



possible, would be expected to produce a double pole in the appropriate S-matrix element. 

A useful "rule of thumb" is that the order of pole a diagram gives is equal to the number 
of "degrees of freedom", e.g. the number of internal lengths in the diagram which can be 
independently adjusted without destroying it. For example, in the bound state diagram, the 
only internal length was the bound state line, but this could be made as long or short as 
desired without problems. Similarly, in figure PH the upper or lower triangles can be scaled 



independently. 

The origin of this rule lies in the fact that, when the Feynman integral of a diagram with 
P internal propagators and L loops is calculated, it turns out to give a pole of order p = 
P — 2L. (Further details can be found in ||.) We now need to apply Euler's well-known 
formula vertices — edges + faces = 1 for any closed diagram, i.e. V — P + L = 1, to get 
p = 2V — P — 2. Each of the V vertices is of three-point typeQ, and each propagator is attached 
to two vertices, except for the four external ones (which are not counted in P), so P = 3 ^ 4 . 
Thus, p = L + 1 = \V = Ap. 

The easiest way to proceed from here is to consider this purely as a problem of topology, and 
start with the diagram without external legs (i.e. with 4 2-point vertices and V — 4 3-point ones), 
then successively remove 2-point vertices and their attached propagators. Since the position of 
these vertices is dictated by the other vertices present, this cannot change the number of degrees 
of freedom. Once this procedure has been exhausted, we can continue by removing the 1-point 
vertices (together with their attached propagator), at the cost of one degree of freedom per 
vertex. Proceeding in this way, we eventually end up with a diagram containing only 3- or 
0-point vertices. A closed network of 3-point vertices can have had no propagators or vertices 
removed from it during the above process, and so, if present, must have existed as a disconnected 
set in the original diagram. Since this is not possible, and since such a network would, in any 
case, not permit momentum to be conserved at each vertex, we must in fact have only 0-point 
vertices, i.e. isolated points. We still have an arbitrary choice of origin to make, and so will 
choose to locate it at one of the vertices. Each of the remaining points can then be moved freely 
and independently, giving the diagram two degrees of freedom per vertex. 



4 By counting the number of faces as the number of loops, we have implicitly taken the points where two particles 
collide but do not bind not to count as vertices. This is different to the usual interpretation of Euler's formula, 
but not inconsistent with it. By taking the diagram to exist in three dimensions, a topological transformation 
can be applied to remove the "extra" vertices and faces. 
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Removing a 2-point vertices and b 1-point ones leaves V — a — b vertices and P — 2a — b = 
propagators. Allowing for the b degrees of freedom which were lost along the way, this implies 
that the original diagram had 2(V — a — b — 1) + b = 2V — 2a — b + 2 degrees of freedom. 
Using the fact that there can be no propagators left, this is just 2V — P — 2 = p. For later 
reference, note that this argument depends only on the fact that no initial vertex is of any more 
than 3-point type, and not on the fact that all vertices are of this type, as the first results do. 
This means that, although calculating the order of a diagram just by halving the number of 
vertices is probably the easiest approach in the bulk, using the number of degrees of freedom is 
a more generally applicable method. (Note, also, that it makes no reference to the integrability 
or otherwise of the theory.) 

Through this method, a pole of any order can be explained in terms of a sufficiently exotic 
on-shell (or Landau) diagram. The one remaining problem is that the only such diagram which 
could ever explain a simple pole is the formation of a bound state. If we are to argue that this 
does not always happen, we have to find a process to take its place. 

Perhaps the most obvious way that the order of the diagram could be reduced would be if it 
so happened that one of the "internal" S-matrix elements had a zero just at the right rapidity. 
However, even if this does not happen, the order can still be lower than expected. 

The explanation for this is that there is not necessarily just one diagram which can be drawn 



to fit a given pair of incoming particles. For example, in figure 1.10 , the theory might allow 
a different set of particles to be used for the internal lines, e.g. substituting anti-particle for 
particle on each line in the upper or lower triangle, without disturbing the diagram. It is even 
possible that an entirely different diagram could be drawn to fit the same external lines. In such 
a case, all possible diagrams must be added together with appropriate relative weights. If a 
cancellation occurs between the different diagrams, then the overall order of the pole produced 
is lower than what would be expected for any of the diagrams individually. For our example, if 
this sum came to zero, then they would collectively contribute a simple, rather than a double, 
pole. 



1.3 Boundary field theory 

The theories we have been considering so far have lived on the "full line" stretching to infinity 
in both directions. Many interesting new features arise if we insert a "wall" at x = to restrict 
the world to the "half line" between zero and negative infinity. Far away from the wall, particles 
behave in exactly the same way as before but, when the approach the wall, two things can 
happen. Either they will reflect off it, or they will bind to it, forming a boundary bound state. 
The introduction of a wall is thus not just a simple matter of geometry, and a boundary analogue 
of the S-matrix, termed the "reflection factor" must also be introduced. 



This idea was first introduced by Cherednik [19 1 , though it took 10 years or so for the topic 



to be put on a footing comparable with the bulk theory. This was achieved by Ghoshal and 



Zamolodchikov [13], as well as Fring and Koberle |20| and Sasaki [21]. A good review of the 



topic can be found in 22]. 



In algebraic terms, Ghoshal and Zamolodchikov imagined the ground state of such a theory — 
|0) b — as being formally created from the ground state of the bulk theory by a "boundary creating 
operator" B, creating an infinitely heavy and impenetrable particle B sitting at x = 0. Thus 

\0)b = B\0). (1.38) 

While this is a purely formal object, it makes analogy with the bulk theory straightforward. Far 
from the wall, everything is exactly the same as for the corresponding bulk theory, allowing the 



1.3 Boundary field theory 



15 




same set of asymptotic particle states, so an in-state of the boundary theory is just 

A ai (O^Aa, (0 2 ) ...A aN (9 N )\0) B = A ai (6 1 )A a2 (9 2 ) . . . A aN (9 N )B\0) , (1.39) 

with 0i > 6> 2 > • • • > N > 0. 

By analogy with the bulk S-matrix, they then introduced a reflection factor to interpolate 
between the in- and out-states through the relations 

A a (0)B = R b a (9)A b (-9)B , (1.40) 

illustrated in figure |1.11| . 

Following the previous discussion, we will consider the boundary version of integrable theo- 
ries. This means that the introduction of a suitable "wall" will involve modifying the action by 
adding a boundary potential term which will restrict the particles to the half-line, but also leave 
the theory integrable, allowing us to still have the useful features of factorisability and lack of 
particle production. Importantly, this means that only the 1 — ► 1 reflection factor will need to 
be considered. 

Assuming such a wall can be built, the logical next step is to search for boundary analogues 
of the four conditions placed on the S-matrix above. Of the three symmetries enjoyed by the 
S-matrix — P, C and T — only time-reversal symmetry inevitably remains, demanding R^(9) = 
R b a {9). The presence of charge conjugation symmetry is generally permitted by some choices of 
boundary condition, but is not inevitable, as we shall see later. Finally, parity symmetry must 
inevitably be broken by the introduction of a wall of any type. The four S-matrix conditions, 
however, all have analogues for the boundary, and are sufficient to specify the reflection factor 
up to a boundary CDD ambiguity which satisfies the same constraints as for the bulk. 



1.3.1 Boundary Yang-Baxter equation 

The demands of factorisation again require that trajectories should be able to be moved past 
boundary vertices, i.e. the points where particles interact with the boundary. This is shown in 



figure 1.12, or algebraically as 



Kl(e2)si\%(9 l + 9 2 )Rii{e 1 )s b iX 



S c a \Z(0i - 02)R d c l(9 1 )S d c f/ i (9 1 + 9 2 )R%(9 2 ). (1.41) 
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Figure 1.12: Boundary Yang-Baxter equation 



1.3.2 Boundary unitarity condition 

This is again a straightforward generalisation of the bulk requirement, and results in the condi- 
tion 



R c a (9)R b c (-9) = <t 



(1.42) 



Algebraically, this results from the demand that the reflection factor should also be analytic, 
and so (|1,40 ) should make sense for negative 9. The argument then proceeds in the same way 
as for the S-matrix. 



1.3.3 Boundary crossing symmetry condition 

This time, trying to find a boundary analogue is somewhat more tricky, and in fact it turns out 
to be easier to find a "boundary cross-unitarity" condition 



K ab {9) = S%(29)K b ' a '( 



b'a' 



K a \9) = 4 ( ? 



where 



In terms of the reflection factor, this can also be written as 

R a b (9) = S%(20)R$(i7r-9). 



(1.43) 



:i.44) 



(1.45) 



1.3.4 Boundary bootstrap 

With the introduction of the boundary, there are now two types of bound state to consider: 
bulk "bound state" particles, and "boundary bound states". The second type arise due to an 
incoming particle binding to the boundary, changing its state, as shown in figure [1.13 . For a 
particle a changing the boundary from state a to state (3, we can define a boundary fusing angle 



with a corresponding pole in the reflection factor at iu 



(3 

aa ■ 
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Figure 1.13: Boundary bound state Figure 1.14: Process involving a bulk 

and a boundary coupling 

This also leads to the introduction of a set of boundary couplings g c . Again, the reflection 
factor at a boundary fusing angle can be considered as a pair of boundary couplings, giving 



z - lUaa 



Alternatively, if the particle c can be formed as the bound state of two equal-mass particles in 



the bulk theory, we would expect the process shown in figure 1.14, giving 



29 



K ab (6) w - / c a c . (1.47) 



All this allows us to play a similar game to before to determine the boundary spectrum. 
Assuming that all boundary states other than the lowest (vacuum) state can be formed by the 
binding of a bulk particle to the boundary, and that we can somehow construct reflection factors 
for the vacuum boundary state for all the bulk particles, we can search their pole structures for 
evidence of further boundary states. Constructing a new set of reflection factors for these states, 
searching again, and repeating until all the poles in all the reflection factors can be accounted 
for without introducing further boundary states, we can hopefully obtain the entire spectrum. 
As before, this relies on introducing no more states than are necessary to complete the process, 
which might, in theory, mean some are missed. However, that has so far never been found to 
happen in practice. 

For the bulk bound states, factorisation demands that we be able to move the boundary 



interaction past the bound state formation vertex as shown in figure L15 , leading to 

ffR d c (9) = fl^R%{6 + irf ad )S b b l a a2 {26 + «4 - i<4)<(0 - it&). (1.48) 



A similar demand for the boundary bound states leads to figure |l.lq , and the corresponding 
requirement 

g% a R c d p(9) = gfsf e (9 + iuP h )Rf ga (e)Sg(9 - iufa) . (1.49) 
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1.3.5 The boundary Coleman-Thun mechanism 

Though the discussion is essentially analogous to that of the bulk, the Coleman-Thun mechanism 
becomes increasingly complicated with a boundary present plf . This is because, as well as the 
processes which were possible in the bulk, a new set become possible involving the boundary 
reflection factors. It is even possible to formulate on-shell processes which involve cancellations 
between bulk S-matrix elements and boundary reflection factors. One important result which 
does, however, remain true is that the naive order of an on-shell diagram is equal to the number 
of degrees of freedom. This (or alternatively using p = 2V — P— 2) is perhaps the most useful way 
of proceeding, now that there will be a mixture of 3-point bulk vertices and 2-point boundary 
ones present. 

There are two types of process: ones which involve the boundary vertices ("boundary de- 
pendent") and those where the only interaction between the particles and the boundary is to 
reflect from it ("boundary independent"). Reflection factors, in general, have a factor which is 
independent of any boundary parameters present, but which nonetheless contains simple poles. 
Without the Coleman-Thun mechanism, such poles would have no explanation, since any pole 
which was due to the formation of a bound state with the boundary would be expected to depend 
on the properties of the boundary. 
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Figure 1.17 shows two possible boundary independent processes. In many models where two 
equal-mass particles can form a bound state at relative rapidity a, it would be expected that 
the reflection factor would have a pole at ^(ir — a), explained by the left-hand diagram. The 
right-hand diagram shows a more involved process, which relies on having a suitable bulk vertex. 
The important point to note here is that, to make the triangle close, the angle of incidence on 
the boundary cannot depend on any of the boundary parameters. This means that none of the 
boundary-dependent poles come into play, and so there is always just a simple reflection from 
the boundary. 



Some common boundary dependent processes are shown in figure 1.18. If an incoming 
particle with rapidity 5 forms a boundary bound state, then there will always be a pole in 
the reflection factor at 5 for the same particle on the new state, explicable by the left-hand 
diagram. The boundary initially emits the particle that helped to create it, being reduced to 
the original state in the process. The incoming particle then re-creates the new state. The 
other two diagrams simply rely on there being suitable boundary and bulk vertices to make 
them close. They are naively second order, but could be reduced to first order if the boundary 
reflection factor had a zero at the appropriate rapidity. 

The rightmost diagram is the most interesting, since it can be reduced to first order either 
by a zero of the reflection factor, a zero of the S-matrix element or (depending on the theory) 
cancellation between diagrams for different arrangements of the internal loop. It is this last, in 
particular, which shows how delicate the relationships between the S-matrix and the reflection 
factors need to be to effect the result. 

Another point to make about boundary poles is that they can go from describing a bound 
state to being due to a Coleman-Thun process by a tuning of the boundary parameters. Often, 



at the point where this happens, a process like figure LIS becomes possible. This is a modified 
version of the right-hand diagram in figure |1.17| , where the boundary parameter has been ad- 
justed to make the particle reflect from the boundary at a pole. As the parameters are tuned 
on through this point, the diagram then collapses into a CT process such as the middle diagram 



of figure 1.18 . While there is no general proof, it appears to be true for the sine-Gordon model 
at least that such a collision of boundary-independent and boundary-dependent processes must 
happen for a pole to cease to be due to a bound state. 
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Figure 1.18: Some common boundary dependent Coleman-Thun processes 




Figure 1.19: Coleman-Thun process possible 
only at special boundary parameter values. 



1.4 Summary 

This chapter has provided a brief overview of the world of 1+1-dimensional integrable quantum 
field theory, and some of its most interesting features. The restrictions imposed by integrabil- 
ity make the axiomatic approach immensely powerful, allowing exact S-matrices to be found; 
this is the only arena where such results are possible at present, underlining its importance 
in uncovering non-perturbative results and pointing the way for tackling more realistic field 
theories. 



Chapter 2 

Classical sine-Gordon Theory 



"All these have never yet been seen — 
But scientists who ought to know, 
Assure us that they must be so... 
Oh! let us never, never doubt 
What nobody is sure about!" 

— Hilaire Belloc 

2.1 Introduction 

"First, establish a firm base." 

— Sun Tzu 

In the next chapter, we will study the effect of introducing a boundary into the sine-Gordon 
theory (which, as noted in the introduction, is the simplest ATFT). Before plunging ahead with 
the full quantum theory, however, it is worthwhile to take a look at the classical limit. This 
exhibits essentially the same features as the quantum theory, but in a form that makes it much 
easier to gain a direct understanding of what is going on. 

To take a step even further back, the first section discusses the classical theory without a 
boundary, attempting to motivate the idea that it possesses an infinite number of conserved 
charges, and so is integrable, with all the simplifying features that entails. While not being 
a proof, it will offer a means of calculating as many conserved quantities as desired. It will 
also help to show how "special" the sine-Gordon theory really is: it is one of only two possible 
integrable field theories with a single scalar field. 

Since it is not at all obvious that the introduction of a boundary should preserve many of 
these conserved quantities (let alone the infinite number required to maintain integrability) the 
restrictions integrability imposes on the possible boundary conditions will then be examined, 
and the most general integrable boundary condition found. 

To complete this chapter, and present a physical picture to take into the next, the first few 
classical boundary bound states will be constructed by the method of images. The idea — which 
is perhaps familiar from its use in electromagnetism — is that a given process on the half-line 
can be described by the theory on the full line with a set of "image" particles placed behind the 
boundary. This can indeed be done, for any integrable boundary condition. Lastly — and with 
the benefit of hindsight — we will use this to make some predictions for the full quantum theory, 
smoothing the path ahead. 
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2.2 The bulk theory 



The classical action for the theory on the whole line is 

/oo poo 2 

dx / dt \{d^f - -jf (cos(/^) - 1) , (2.1) 
-oo J — oo P 

where mo sets the mass scale and (5 is the coupling constant. The particular form of the potential 
term gives the theory its integrable properties, so let us, for the moment, consider a more general 
theory with a potential — &V(<p) so that we can investigate how special it really is. The following 
argument was first made by Ghoshal and Zamolodchikov O]; the form given here is taken from 
©■ 

To simplify the notation, it helps to use light-cone co-ordinates, defined through d± = 
\{dt ± d x ). The equation of motion — dA = — then becomes d+d-(p = —V'(<p). 

To construct conserved densities, imagine that there exist two quantities, T and 0, such that 
d—T = d+Q. Rewriting this using x and t, we find 

d t (T-Q) = d x (T + e), (2.2) 

d 



dt 



J dx(T - 6) 



[T + @r oo = 0, (2.3) 



showing that J dx(T — O) is conserved. The search now is for suitable quantities T; here, we 
will focus on polynomials in d+<p,d+ip, . . . , and go order-by-order in the total number of +- 
derivatives. This number will be denoted as s + 1, with T s and Q s standing for Ts and 0s with 
s + 1 +-derivativesQ. The conserved charge will then be annotated as 



Q s 



Jt s+1 -G s _iGfe, (2.4) 



where s can now be seen to stand for the spin of the charge. 
Three other points are worth noting: 



• total d + derivatives can be dropped; 

• a polynomial in which no term has its highest derivative factor occurring linearly can never 
be a total 9+ derivative; 

• for each T s+ i, there is a corresponding T_ s _i, obtained by interchanging d + and cL 
throughout. 



Looking first at s = 1, we find: 

d-T 2 = 2{d + ip)d + d.ip 

= -2(d + <p)V(<p) (2.5) 
= d + [-2V(<p)}, 

showing that (<9+(^) 2 and — 2V(ip) provide a suitable pair for any potential V. This, in fact, is 
not surprising, since Q\ + Q_i is just energy, and Q\ — Q-\ is momentum, two quantities that 
are always conserved. 

Focusing simply on polynomial functions, these will turn out to be unique. 
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There is no solution for s = 2, and the first nontrivial result appears at s = 3, where 

T4= (f) + (2.6) 

provides a solution for any real or imaginary /?, but only if V'" = (3 2 V' . This has the solutions 

(3 = : V = A + B(ip-ipo) 2 , (2.7) 
/3^0 : F = A + Bf&' + Ce-to', (2.8) 

for any constants A, B and C. If /3 = 0, this corresponds to either the massive or massless 
free field theory, depending on whether or not B is non-zero. For (3 ^ 0, we get the (massless) 
Liouville model if B or C is zero. Otherwise, it is the sinh-Gordon model if (5 is real, or sine- 
Gordon if it is imaginary. 

If we were to proceed with this, we would find only one more model with any conserved 
charges above s = 1, namely the Bullough-Dodd model, which appears at s = 5. However, 
the sine-Gordon model would turn out to have conserved charges at all odd s, which is the 
crucial pointQ. (In fact, since Parke's argument shows that any model with a conserved charge 
above s = 1 must have the properties needed to follow through the exact S-matrix approach, 
we already have all we need.) This, of course, still does not answer the question as to why this 
should be true. A better understanding can be gained once the sine-Gordon model is thought of 
as an ATFT, with an underlying Lie algebra structure. It is this structure that endows it with 
the symmetry that the charges flow from. 



2.3 The theory on the half-line 

To introduce a boundary into the model, we must impose a boundary condition on the field, 
implemented through the addition of a boundary term to the action, i.e. 

/oo 
dtB(<p B ) , (2.9) 
-oo 

where <fB(t) = y(0, t), and so depends only on the value of the field at the boundary. The term 
-4buik is 

/oo roo 
dx dt \{d^f -4V(<p), (2.10) 
-oo J — oo 

and we are assuming that the bulk potential has been chosen so as to make the bulk theory 
integrable. 

The equation of motion is the same as before (though restricted to apply only to the half-line) 
but the new term introduces the boundary condition 

d x (f\x=o = -B'(<p B ) ■ (2.11) 

Clearly, not all the conservation laws from the bulk model can still apply now that we have 
introduced a boundary (momentum, for example) but it still turns out to be possible to keep 

2 In practice, the existence of an infinite number of conserved charges was proved via the inverse scattering 
method Q. 
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a (possibly infinite) subset. Working by analogy with the argument for the bulk, the problem 
arises because (12.31) is modified to 



dQs 
dt 



d_ 

dt 



s+l 



@s-l) 



[t s+1 + e,_i]_ c 



(r s+ i + e 



s-l)\x=0 



(2.12) 



The only way this can be saved is to demand that the rhs is a total i-derivative, allowing it to be 
incorporated into the lhs to give a new quantity which is conserved. (For the quantities found 
so far, the Ts are t-derivatives, whereas the 0s are not.) 



In general, 



^- t (Q s ±Q-s) 



{T s +i T T-s-i =F ©-s+l — @s-i)U=o ■ 



(2.13) 

(2.14) 
(2.15) 



This explains why momentum — Q\ — Q-\ — is not conserved (the final line reading (T 2 + T_ 2 
2@o)U=o)- For ener gy) on the other hand, we find 



d_ 
dt 



(Qi + Q- 



(T 2 -T_ 2 )U =0 

dt<pd x tp\ x=0 

-d t ipB'(ip B )\x=o 

|[-*(w0], 



(2.16) 



by making use of the boundary condition. This gives us 



d_ 

dt 



dxa{d t ip) 2 + \{d xV ) 2 + 4V(<p))dx + B&! 



0. 



(2.17) 



showing that energy is indeed still conserved on the half-line. The next natural step is to ask 
whether a B can be found that allows modified versions of all charges of the form Q s + Q- s 
to still be conserved. From above, this is true if (T s +i + O s _i — T_ s _i — 0- s +i)| x =o is a total 
t-derivative. 



Imposing this restriction on the s = 3 charge of the bulk theory, we find B'" 
whose most general solution is 



B', 



B 

B(ip B ) = Mcos -(ip B - ipo) 



(2.18) 



for some constants M and <pq. The similarity of this solution to the requirement on V in the 
bulk theory makes it reasonable to imagine that, just as the bulk potential allowed conserved 
charges for all higher odd s, this form for B should too. This was finally proved for the classical 
theory through the inverse scattering method 25], where it was found that this is the most 
general integrable solution for B. It was also found that all the charges discussed here — "even 
parity" charges from the bulk theory modified by a boundary term — survived with this boundary 
condition. 



2.4 Particle content 



Having established that both the bulk and boundary theories are integrable, the next step is to 
find out what the theories actually describe. Due to the periodic nature of the potential, the 
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Figure 2.1: Single soliton solution 



sine-Gordon model is unusual in having an infinite number of vacua, at 2™ for any n G Z. The 
"particles" of the theory therefore turn out not to be the usual localised humps in the field, but 



rather a configuration that interpolates between two neighbouring vacua, as shown in figure 2.1. 



This configuration has two useful properties. First, it is a "soliton", which means that it 
preserves its shape over time without dissipation or decay. Secondly, because it interpolates 
between vacua, it cannot be destroyed as that would alter the value of the field as x — > ±oo; 
for this reason, the theory is often called "topological" . The "topological charge" of a state is 
defined as the difference (in units of ^) between the value of the field at — oo and +oo and 
must be conserved. A single soliton state has charge 1, while an anti-soliton state (interpolating 
between the vacuum at — oo and the next lower one) has charge -1. 

If two (anti-)solitons collide, they are simply transmitted through the collision, without 
changing shape, and so can truly be considered as particles (which, the theory being integrable, 
cannot be created or destroyed). If their rapidities are allowed to be complex, however, rather 
than purely real, the situation changes. If a soliton and an anti-soliton are given conjugate 
rapidities, a "bound state" appears. Due to the periodic up-and-down motion of the field, these 
particles are known as "breathers" and are categorised by the imaginary component of their 
rapidity, which determines their period and mass. 

The soliton and anti-soliton both have the same mass, which we shall call m s , while the 
mass of the breather formed by a soliton-anti-soliton pair at a relative rapidity of 9 = iu is 
m = 2m s cos (^). That completes the particle spectrum of the bulk theory as, if two breathers 
are persuaded to bind together, they simply form a third breather. 



2.5 Construction of boundary bound states 

Rather than try to analyse the boundary theory in a similar way to the bulk, it is easier to use 
the method of images. The idea is to find a particular configuration of the bulk theory where 
the value of the field at (p(0, t) happens to obey one of the integrable boundary conditions. The 
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left-hand half line then provides a solution to the boundary theory. 

The vacuum state of the boundary theory simply requires that <p(0, t) = ipQ for all time, 
whatever the value of M, so a suitably-placed stationary soliton is all that is required. By 
analogy with electromagnetism, we might then imagine that the boundary state consisting of 
n particles with rapidities 81,62,... ,9 n corresponds to the bulk state with an "image" set of 
particles behind the boundary (with opposite rapidities) and, again, a stationary soliton near 
the boundary. 



This problem was first tackled by Saleur, Skorik and Warner [26], who found the 3-soliton 
solution. The choice as to whether each particle was a soliton or anti-soliton and their relative 
initial positions selected which boundary condition was obeyed. In addition, in the Neumann 
limit the position of the stationary particle became infinite, reducing the result to a two-soliton 
solution. 

The natural generalisation of this is to consider a 2n + 1-soliton solution (which reduces to 
a 2n-soliton solution in the Neumann limit). This is easier than it might appear as, in the limit 
where the particles are well separated (i.e. t — > ±00), the state of the field at the boundary is 
determined only by the central stationary soliton and the two moving solitons that are closest, 
allowing the 3-soliton solution to be re-used. 

In addition, as SSW found in the 3-soliton case, the absolute positions of each pair of moving 
solitons are irrelevant in the solution of the boundary condition; it is only the phase delay that 
is important. Using this fact, any pair of solitons in the 2n + 1-soliton solution can be moved 
off to infinity (provided their phase delay is preserved), reducing it to a 2n — 1-soliton solution. 
Using these two facts, it is perhaps beginning to become clear that the general solution can be 
built out of the 3-soliton solution with a little cunning. 



2.5.1 Notation 

For the classical problem, it is convenient to re-scale the field and coupling constant to re-express 
the bulk sine-Gordon equation as 

4>tt - 4>xx = - sin(0) , (2.19) 
where 0(p = (p. On the half line, the most general boundary condition then becomes 

^U=o = Msini(^-0 o )U=o- (2.20) 



The classical multi-soliton solution for the whole line has been known for some time, and is 



generally expressed in terms of Hirota's r-functions [27] as 



(x,t) = 4arg(r) = 4arctan 



where the r-function for an iV-soliton solution is 



3Q0 
»(t) J ' 



T (x,t)= ^2 e " ^= ie ^) exp 
/Uj=0,l 



E^i{(*i + j^) a? +( fc i-^)* 



l<i<j<N v 3 



(2.21) 



. (2.22) 
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The parameters fcj are related to the soliton rapidities by fcj = e 1 , so the solitons' velocities are 
given by 

- (frr) ' < 2 ' 23 ' 

The ai represent the initial positions of the solitons (but see below) while the are +1(-1) for 
solitons (anti-solitons). 

For the sake of simplicity, we shall number the particles in decreasing order of rapidity, so 
that particle 1 has the highest rapidity, particle 2 has the next highest, and so on. This ensures 
the logarithm in the r-function is always real. Other orderings give the same result — as they 
clearly must — but it is less transparent that the r-function is real. 



2.5.2 The position problem 



Before going any further, a problem immediately arises with the interpretation of the a, as the 
positions of the solitons. If this was truly the case, for example, a 3-soliton solution as t — > ±oo 
would reduce to a single-soliton solution with the same value of the position parameter. This, 
however, is not true. The one-soliton solution is just 



leading to 



Ti(x, t) = 1 + ie± exp 



4>i(x,t) =4 arctan ( t\ exp 



a 1 



ai 



(2.24) 



(2.25) 



Taking the 3-soliton solution, note that, as t — > oo, the soliton with positive rapidity will 
contribute a highly negative exponential whenever it appears in the sum, whereas the one with 
negative rapidity will contribute a correspondingly positive exponential. From this, it is clear 
that the two dominant terms will therefore be the ones where fi\ = and y 3 = 1. Thus, as 
t — > 00, 



T 3 (x,t) 



e 3 exp 



1 



x + [k 3 



1 — €2 exp 



-x + y + 2 In 



1 



1 

h 



t — do 



leading to 



4>3(x,t) 4arctan ( — €2 exp 



x-^-21n 
2 



l-k 3 



If we now remember that tan (a; ± ^) = — tan(x" 
re- written as 4 arctan — - = y ± 2ir. Thus, we find 



i + h, 

-1 , this implies that 4 arctan x 



03 (x, t) ± 2tt « 4 arctan I ei exp 



(in 

-x + + 2 In 



l-fc 3 
l + k 3 



(2.26) 

(2.27) 
y can be 

(2.28) 



The 2ir on the lhs, which is equal to the spacing of the vacua, just represents the fact that the 
Hirota solution imposes = at —00, while the natural assumption here is that = ±2ir, so 
that the leftmost soliton reduces the field to zero heading in towards the central particle. Thus, 
we do indeed have a single-soliton solution, but with 



a i = o! 2 + 4 In 



(2.29) 
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Repeating this exercise with t — > — oo instead gives the same result, but with k% replaced by 
Since we would like the stationary soliton to solve the same boundary condition in both cases 
(as this condition stays unchanged for all time), it is clear that we are forced to take ks 
as SSW did. 



j_ 

ki ■ 



The reason for this is easy to see once it is realised that, to shift the solution in time, all 
that is needed is to shift each position parameter by velocity x time, irrespective of any collisions 
which may have happened in the interim. The effects of collisions are thus built into the solution, 
and the parameters are only indirectly related to particle positions at any given time. In what 
follows, however, it will be easier to work in terms of "actual" parameters, and transform back 
to Hirota's parameters at the end. 

A more general analysis shows that the 2N + 1-soliton solution with N pairs of solitons with 
opposite rapidities examined at t — ► — oo reduces to 2N + 1 single-soliton solutions as expected, 
but with each soliton position modified by a term involving all rapidities higher than its own. 
This means that the "position" parameters aj only have a genuine interpretation as a position 
for the particle with the highest rapidity. As t — > +oo, the opposite situation arises, with the 
position parameters modified by all lower rapidities. To be precise, let us take x — > +oo and 
t — > — oo with j Vi, to keep ourselves in the neighbourhood of particle i. This means that 
i > N (we are considering the particles with negative rapidities). Then, by the same reasoning 
as before, the two dominant terms will be those with fij = 1, j < i and fij = 0, j > i. This gives, 
for i — 1 even, 



i-1 



T2N+i{x,t) « (-1)* 2 JJej-exp 



i-i 



1 + iei exp 



+2 E 

l<i'<j<i-l 
1 
2 



kj + — J x + ( kj — — | / - a , 



In ' h ' 



leading to 
4>2N+i(x,t) ~ 4arctan I 



+2 £ In 

i<j<»-i 



/■•, + ^-j x+ (ki - — ] / - rt, 



fej ki 
ki -\- kj 



(2.30) 



ki ~\~ kj 



(2.31) 



(2.32) 



+2 ^ In 

l<J<i— 1 

Thus, compared with the appropriate single soliton solution, 

a 1 = a, + 4 £ ln(|^| 

Note that, this time, <p —>■ as x — ► oo is the natural situation, in agreement with the Hirota 
formula. For 2A^ + 1 — z odd, we need to use the same trick as for the three-soliton case, but 
finish up with the same formula. Finally, for i < N we need to take x — > — oo but the result 
remains true. Taking the other limit (as t — > +oo), the analogous result is 

a 1 = a l + A £ Inf^V (2.33) 



«+l<i<27V+l 



ki -\- kq 
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Note also that the only way to ensure the boundary condition stays constant in time is to impose 

h = l/k 2 N+2-i- 

From this, we can calculate the phase delay between any pair of particles with equal and 
opposite rapidities — i and 2N + 2 — i — in terms of their "actual" position parameters a' n as 
t — ► — oo as 



0>i + 027V+2-1 — 0>i + 0-2N+2- 



4^> 



kj k^ 



k% H~ kj 



(2.34) 



assuming i < N. 



2.5.3 Solving the boundary condition 

As a warm up to the general solution, it is useful to consider the simplest possible solution, with 
only one stationary soliton. This corresponds to the ground state of the boundary model. In 
this case, all boundary conditions reduce to the demand that (f>\ x =o = ^o> f° r some constant <£o 
depending on the boundary conditions. (In the Dirichlet case, <3?o becomes simply <pQ.) Putting 
this into ( |2.25 ), we find 



•1 



01 (0, t) = 4 arctan (eie^j = $ , (2.35) 
implying 

ai = 2 In ^ei tan . (2.36) 

Taking advantage of the fact that the 3-soliton solution must tend to this near the origin as 
t — > —00, we can immediately write down the position parameter of the stationary soliton in 
the 3-soliton solution as 

-> tHKstt)-" 1 Ht(!^) 2 ) ^ ^ 

Noting that, in terms of the rapidity variable, ( ^~^ ) 2 = tanh (|) 2 , this agrees with the formula 
for given by SSW in their appendix. They derive this specifically for the case of Dirichlet 
boundary conditions, but it can now be seen to have the same form for all boundary conditions. 



2.5.4 The general solution 

By extension of the above argument, the position parameter of the stationary soliton in the 
2N + 1-soliton solution is 

a N+1 = 2 In ( eyv+itan-?^ J[ tanh (?£\ J . (2.38) 

As has already been mentioned, the general 2N + 1-soliton solution reduces to the 3-soliton 
solution (involving the 3 slowest solitons) as t — > —00, so the phase delay for the slowest pair 

3 Their as is twice the 0,3 which appears in the Hirota formula, accounting for their loss of the factor of 2 in 
front of the logarithm. Also, they consider the left half-line rather then the right. 
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should be given by the SSW formula, which (with our conventions) is 



2 In < 



-6163 tanh 



tanh(0)~ 



tanh \ (0 + irf) tanh \ {9 — irf) 



tanh i(0 + C) tanh i(0-C) 
where r\ and £ are the solutions of the simultaneous equations 



±1 



(2.39) 



Mcos(i< 
Msin(i ( 



2 cosh £ cos 77 
2 sinh £ sin rj . 



(2.40) 



The ambiguity in the sign of ( 2.39|) is simply a vagary of the solution method (due to the fact 
that the bulk vacua are 27r-periodic, whereas the boundary vacua are only 47r-periodic; working 
in terms of the bulk makes the stable and unstable possibilities appear together). We shall 
concentrate on the negative sign, which corresponds to the stable boundary value. 

By virtue of ( 2.34| ), this can be re- written using the "actual" position parameters instead, as 



2 In 



-eie 3 



tanh \ ( 



+ irf) tanh h(6 — irf) 



tanh ~(0 + C) tanh M8 - C) 



±1 



(2.41) 



Turning now to the faster particles, we need to use the fact that, for the slowest particles, only 
the phase delay is important. This means that we can take their actual positions off to ±00 
without affecting the validity of the solution, and essentially reduce the problem to the 2N — 1- 
soliton case. Now, the next slowest particles have gained the mantle of being the slowest, and 
so must have a phase delay of the same form. (Note that, in doing this, we have made the 
slowest particles collide with all the faster ones in turn on their way to infinity, changing their 
positions. The symmetry of the situation, however, ensures that the phase delays between the 
pairs of particles stay intact.) 

Repeating this for all the particles shows that, for each pair, all that is relevant is the phase 
delay, and this always has the SSW form. In terms of the position parameters, we then have 



2 In I 



-eie 2 N+2-i JJtanh 



tanh 2 (0 + irf) tanh ^ ( 
tanh \(6 + Q tanh |( 







±1' 



(2.42) 



where a* = a, + a,2N-\ 



This completes the solution, but for one point: through this argument we have shown that 
if the 2N + 1-soliton solution exists then it must have the given form, but we have not shown 
that it actually exists. For that, we would need to substitute the results back into the Hirota 
formula to check — a cumbersome task, and one for which we lack the energy. In the meantime, 
we content ourselves with the observation that it seems a reasonable assumption, and bears up 
to all the numerical checks we have carried out. 



2.6 Boundary bound states 
2.6.1 Boundary breathers 

The natural progression from this is to consider extending the rapidities to complex values. 
While this can be used give a solution where breathers rather than solitons interact with the 
boundary, it can also be used to construct "boundary breathers" or boundary bound states. 



2.6 Boundary bound states 



31 




Figure 2.2: 3-soliton solution, period 10 

These solutions arise when the pair of particles which are given complex conjugate rapidities 
consist of one in front of the boundary and one behind. 

Due to the requirement that their rapidities must also be equal and opposite, this implies 
that they must be given purely imaginary rapidities. Curiously, in the Dirichlet case (as SSW 
noted) the pair of particles must also consist of two solitons or two anti-solitons, not a particle 
and its anti-particle as in the bulk. 

Simply by continuing all rapidities to imaginary values, we can generate a sequence of bound 
states through solutions with successively greater numbers of solitons. The one subtlety is that 
both members of each pair have to be given the same initial position parameter for the solution 
to still obey the boundary condition. This is a consequence of the way the solution was found: 
the r-function was split into real and imaginary parts, assuming all rapidities were real. Making 
some rapidities imaginary disturbs this in general, but putting both members of each pair at 
the same position allows the split used to remain valid. Since all other solutions to the problem 
with real rapidities can be related to this by a time translation, it is reasonable to assume that 
the same is true for the imaginary case. The only difference is that, with imaginary rapidities, 
there is no movement in real space. 

As with the bulk breathers, the period of a boundary breather is given through the imaginary 
part of the rapidity. For the 3-soliton solution, with the moving pair given a rapidity of = iu, 
the period is 27r/sin(«). Now, however, for the breathers coming from higher solutions, each 
pair has its own period. If there exists a common period, whose length is an integer multiple 
of the periods of all the pairs, then the motion is still periodic, but, in general, it will now be 
aperiodic. 



To demonstrate the form of the boundary breathers, we have chosen periodic solutions by 
giving each pair an integer period. These are shown in Figures 2.2 , |2.3j , and 2A for the 3, 5, and 
7 soliton solutions respectively. 
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Figure 2.4: 7-soliton solution, periods 10, 12, and 14 
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2.6.2 Another bound state 

The breathers mentioned above are not the only bound states in the classical theory. The phase 
delay (|2,41| ) becomes infinite at 6 = (, which must be due to the formation of a stable bound 
state. In the bulk theory, this could not happen (all bound states must be formed at imaginary 
rapidities), so it is further evidence of the changes wrought by the introduction of a boundary. 

Considering the 3-soliton solution, we can imagine the particle with negative rapidity as 
being taken off to infinity at this point, leaving just a two-particle process. The remaining 
moving particle sweeps past the boundary, shifting the stationary soliton on the way past; this 
only appears as a bound state when we restrict ourselves to the half line. Then, the incoming 
particle reaches the boundary and disappears, leaving the boundary state changed. 

The final state can be found by considering the limit of the 3-particle r-function where the 
position parameters of both moving particles are taken to infinity. The result of this is that (f>o 
becomes <p' Q , given by 




(2.43) 



2.7 Predictions 

The first prediction to take across to the discussion of the quantum theory is that there should 
be a hierarchy of excited states. For example, the states formed by binding a soliton to the 
boundary should be analogous to the 3-soliton solution found above. After that, further solitons 
should create the quantum versions of 5, 7, . . . soliton solutions. Furthermore, the introduction 
of a breather should allow the formation of a state which could otherwise have been formed by 
two successive solitonic particles. 

A final, and slightly more subtle point, is that the "actual" position parameter used for a 
given pair of particles (with imaginary rapidity) is monotonically decreasing for u < rj, as shown 



in figure 2.5. This means that, in a given solution, the soliton pair with the least rapidity will be 
positioned farthest from the boundary. If we imagine that such a solution, translated into the 
quantum regime, is built up with the soliton finally positioned nearest the boundary interacting 
first, this means that particles must interact in decreasing order of rapidity. 
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Figure 2.5: Plot of a' versus u for rj 



= 1.0 



Chapter 3 

Quantum Boundary sine-Gordon 
Theory 



"The most exciting phrase to hear in science, the one that heralds new discoveries, is 
not 'Eureka!' but 'That's funny ... 

— Isaac Asimov 



3.1 Introduction 



"Mathematicians are a species of Frenchman: if you say something to them they 
translate it into their own language and presto! it is something entirely different." 

— Goethe 



As we saw in the previous chapter, introducing a boundary into the classical theory brings with 
it a number of new phenomena, and in particular a new set of boundary bound states. In this 
chapter, we will investigate these further in the full quantum sine-Gordon model. 

A major complicating factor in this work is the fact that even simple poles in the boundary 
reflection factors should not necessarily be interpreted as being due to the formation of bound 
states, since many have an interpretation through the Coleman-Thun mechanism. Indeed, this 
model provides a good arena for demonstrating the range of possible explanations this mechanism 
can throw up, in some cases involving bulk and boundary matrices working together to induce 
a cancellation in the naive order of a diagram. 

The two main tasks, therefore, are to find suitable interpretations where required, but also 
to find a method of proving that the remaining poles are indeed associated with bound states. 
Two elementary lemmas — which simply serve to impose momentum conservation on boundary 
processes — will turn out to give us all the ammunition we need for this, and should also be 
readily applicable to other models. 

The groundwork for the study of the boundary sine-Gordon model was laid by Ghoshal and 



Zamolodchikov [13], before being taken further by Ghoshal and Skorik and Saleur [28]. They 
provided the basic ground-state reflection factors, and investigated the first few excited states; 
we will take this forward to provide (hopefully) a full and rigorous solution to the problem. 

After reviewing these results in the first section, we will go on to a detailed investigation of 
the Dirichlet boundary condition (where the value of the field at the boundary is fixed for all 
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time). This displays most of the features of the general solution, and will allow us to extend the 
results straightforwardly to all other integrable boundary conditions. 



3.2 Review of previous results 
3.2.1 The theory in the bulk 



As we discussed in the previous chapter, the classical sine-Gordon model (2.1) is integrable. This 
can be shown to be true at the quantum level as well [3C] and so the exact quantum S-matrix 
can be found through the axiomatic program. The essential difference between the classical and 
quantum theories is that the breather particles, which could be formed classically by a soliton- 
anti-soliton pair at any imaginary rapidity, now become quantised. These states can now only 
be formed at relative rapidities of i (it — nir/2\), n = 1, 2, . . . , < A, where 



A- 8?r 



1 



(3.1) 



and so will be labelled as B n . Their mass is therefore m r , 



2m s sin (mr/2\). 

If we denote the soliton S-matrix as S^(9) for rapidity 9, with a, b, c, d taking the value + (— ) 
if the particle is a soliton (anti-soliton), the non-zero scattering amplitudes are S~\z\_{0) = 
SZZ(9) = a{9) (soliton-soliton or anti-soliton-anti-soliton scattering), S+Z{0) = SZ+(9) = b(9) 
(soliton-anti-soliton transmission), and S_T(9) = Sz1l(9) = c(9) (soliton-anti-soliton reflection). 
Explicitly, 



where u - 
p(u) 



-i9 and 



I TT 

sin(A(u — tt)) 



a{9) = sin[A(7r — u)]p(u) , 
b(9) = sm(Xu)p(u) , 
c(9) = sin(A7r)p(n) , 



r ((2/ - 2)A - r (i + 21 x - 
r((2/-i)A-^)r(i + (2/-i)A-^ 



(3.2) 



■/(« 



(3.3) 



As pointed out in [BT[, this factor can also be written in terms of Barnes' diperiodic sine function 
S2(x\uji,u>2) |3^, |33f . This is a meromorphic function parametrised by the pair of 'quasiperiods' 
(ui,U2), with poles and zeroes at the following points: 



poles 
zeroes 

In terms of this function, 

p(u) = 



x = nitui+n 2 002 »2 = 1> % ■ • ■ ) 

x = mi u>\ + m2 oj 2 (mi, ?7i2 = 0, — 1, — 2 . 



1 ^2 (iT — U 


\l,2n)S 2 (u\ 




sin(A(u — 7r)) 5*2 (tt + u \ 


%2n)S 2 [-u 





(3.4) 



(3.5) 



The amplitudes b{9) and c{9) have simple poles at 9 = z (vr — ^-), 
n = 1, 2, . . . , < A, which can be attributed to the creation of B n in the forward channel. There 
are also poles at 9 = in a(9) and b(9) corresponding to the same process in the cross channel. 
Since all poles that we will be discussing, both in the bulk and at the boundary, occur at purely 
imaginary rapidities, from now on we will use the variable u = —i9 and always work in terms of 
purely imaginary rapidities. 
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3.2.2 The theory with a boundary 



Returning again to the previous chapter, the bulk theory be restricted to the half-line x S 



(— oo,0] while still preserving integrability by adding a "boundary action" term [13] 



dt M cos 







(<Pb ~ ¥>o) 



(3.6) 



where M and ip$ are free parameters, and (f B (t) = cp(x,t)\ 



x=0- 



This does not conserve topological charge in general, so four solitonic boundary reflection 
factors need to be introduced, as well as a set of breather reflection factors. The solitonic factors 



which we quote here were given in [13], while breather factors can be found in [28]. 



Solitonic ground state factors 



The reflection factors for the sine-Gordon solitons off the boundary ground state will be denoted 
by P±(u) (a soliton or anti-soliton, incident on the boundary, is reflected back unchanged) and 
Q±(u) (a soliton is reflected back as an anti-soliton, or vice versa). These are given by 



P+(u) = cos(£ + Xu)R(u) 
p- (u) = cos(£ - Xu)R(u) 
Q±[ u ) = |sin(2An) J R(n), 

where 

R(u) = Rq(u)Ri(u). 
The first factor — Rq(u) — is boundary-independent, and can be written as 

T (1 + X(4k - 4) - ^M) T (AXk - 



i? («) = n 



k=l 



r (A(4Jfe - 3) - r (1 + A(4A; - 1) - ^) 



/(« 



(3.7) 



(3.8) 



(3.9) 



The boundary-dependent term is R\(u), given by 

1 



Ri(u) 



cos^ 



a(rj, u)a{id, u), 



whereQ 



a{x,u) 



n [x, f - u) n (-x, f - u) n ( x , -f + u) n (-x, -§ + «) 



n(x,f)n(x,-f)n(-x,f)n(- 



-X 



and 



(3.10) 



(3.11) 



nrr v) = fr r(i + (2/ + i)A + f-^)r(i + (^ + |)A + f) 
1 ' j firft + (2i + i)A + f-^)ra + (2i + i)A + f)- 

The parameters £,77,$, and A; are real and arbitrary apart from being constrained by 



(3.12) 



cos(7?) cosh($) 
cos 2 (77) + cosh 2 (?9) 



icosC 



1 + 



(3.13) 



x Note that there is a small error in Ghoshal and Zamolodchikov's formula (5.23) for a. This corrected version 
was supplied to Patrick Dorey and the author by Subir Ghoshal. 
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Figure 3.1: Vacuum structure 



The relationship of these parameters — which arise in the course of finding the most general 
solution to the four requirements given in Chapter 1 — to the ones which appear in the action 
was, for a long time, unknown. The problem has only recently been solved by Al. Zamolodchikov; 
further details can be found in Appendix |A.2| . These formal parameters, however, are easier to 
work with in practice than the physical ifo and M, and so we shall continue to use them. 

The theory is invariant under ipo — > ipo + ^ , and also under the simultaneous transformations 
ipo — > —ipo and soliton — > anti— soliton. Introducing the boundary breaks the degeneracy of the 



bulk vacua, and selects the lower line in figure |3,l| as the lowest-energy state, with the upper 
line as the first excited state. Continuing tpo through 5 thus simply interchanges the roles of 
these two states, and selects the upper one as the ground state. 

In light of this, we are free to choose ipo to be in the interval < ipo < 5. Note also that the 
topological charge of the ground state is no longer zero, as in the bulk model, but 

q = ^f dx ^ x ^ = " = §7 ' (3,14) 

with the charge of the first excited state being 1 — We will find — at least for the Dirich- 
let case — that all the boundary states have one of these charges so, for convenience, we shall 
designate them simply as and 1 respectively. 



Breather ground state reflection factors 

For the breather sector, Ghoshal [^] obtained the relevant reflection factors — R^ (u) for brea- 
ther n and boundary ground state |0) — from the solitonic reflection factors using the general 



boundary bootstrap equation |2Cj, |13] 



where u n = ir — and the R^ x ^(u) are the solitonic reflection factors, such that R^ x ^(u) is 
the factor for a soliton to be reflected back as an anti-soliton and so on. The f™ b are the bulk 
vertices for the creation of breather n from (anti-)solitons a and b. These obey = (— l) n /™ + . 
The bootstrap is illustrated in figure |3.2|. 
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Figure 3.2: Breather bootstrap 



3.3 The boundary Coleman-Thun mechanism 



To discover the boundary spectrum, the most natural approach is to look for simple poles in the 
reflection factors, which might be expected to be related to the formation of boundary bound 



states. As we have already mentioned in section 1.2.4, however, a complicating factor is the 
fact that not all simple poles correspond to bound states, some having an interpretation as 
anomalous threshold singularities. 

This problem becomes especially serious once a boundary is involved, due to the increased 
complexity of the on-shell diagrams which become possible. This makes it hard to be sure that 
any given pole really does correspond to a new boundary bound state. In the bulk, a simple 
geometrical argument shows that poles in the S-matrix elements of the lightest particle can never 
be explained by a Coleman-Thun mechanism, and so must always be due to bound states [fTTH - 
We wish to find analogous criteria for the boundary situation. To this end, the following two 
lemmas turn out to be useful. Suppose the incoming particle is of type a, and that its reflection 
factor has a simple pole at 9 = iu. 



Lemma 1 Let U a = minb jC (tt — U^ b ). Ifu < U a , then the the pole at iu cannot be explained by a 
Coleman-Thun mechanism, and so must correspond to the binding of particle a to the boundary, 
either before or after crossing the outgoing particle. 



Proof: All processes must take the form shown in figure |3.3| or the crossed version shown in 



figure 3-4 . Conservation of momentum demands that all rescattering must take place within 
the hatched region, which is drawn from the furthest point from the boundary where either the 
incoming or outgoing particle undergoes any interaction. If neither particle decays, we simply 
have a diagram of the form of figure [T^ or figure |3.9| . Otherwise, momentum conservation 
requires that neither product of the particle which decays on the boundary of the hatched 



region has a trajectory which takes it outside that region. Fixing the notation by figure |3Tq (with 
angles XJ\ C and U£ b defined correspondingly), this reduces to demanding 7r — U£ b < u < U h ac . If 
we introduce U a then we must have U a < u < tt — U a (i.e. just u > U a , as u < 5). Thus, if 



u < U a , then the only possible explanations for the pole are figure 3.8 and figure 



Lemma 2 If the boundary is in its ground state, then lemma |^ can be strengthened, requiring 
that the incoming particle bind to the boundary if u is outside the range U a < u < ^ — U a . In 
addition, if min;, iC U bc > § , the incoming particle must always bind to the boundary. 



Proof: With the boundary in its ground state, all rescattering must take place in the area 
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Figure 3.3: General process, with Figure 3.4: General process, with 
incoming particles uncrossed incoming particles crossed 





Figure 3.5: General process 
when boundary is in ground state 



Figure 3.6: Decay process 



shown in figure |3.5| . Reasoning as before but demanding that both product particles be emitted 
into this more restricted region, we find it — U° b < u < U\ c — \ , or U a < u < ^ — U a . In 
addition, both particles b, c must be emitted into an angle of |, so U bc < | for at least one pair 
of particles b, c. If either of these conditions are violated, then the incoming particle must bind 
to the boundary. 



These two results, between them, will allow the spectrum of the boundary sine-Gordon 
model to be fixed completely, provided it is assumed that no pole corresponds to the creation 
of a boundary state if it has an alternative (Coleman-Thun) explanation. 

For the problem under discussion, writing the rapidity bounds U a as U + /\ for the soliton 
(anti-soliton) and as U n for the B n , we have 



U± 



U, 



TT 

2 ~ 

7T 

2A 

TT 



,-TT 



2A 

n ^ n max 



(3.16) 



where B„ 



2 2A 

is the highest-numbered breather present in the model. To derive these results, 
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note that a soliton (anti-soliton) can only decay into an anti-soliton (soliton) and a breather 
(with vertex U^ n = f + ^f)- A breather can either decay into a soliton-anti-soliton pair 
(U™_ = 7r — or a pair of breathers {U! nm = ir — ^ with n = m + Z or m = n + Z, or 

^nm = 7r( l\ m) with i = " + m). 

These restrictions can also be combined to produce a stronger version of lemma 1 when 
the incoming particle is a soliton. If U + < u < j, decay within the hatched region is only 
possible into the topmost breather and an anti-soliton. One or other of these particles will be 
heading away from the centre of the diagram. If the process in uncrossed, as in figure |3.3] , the 
breather will be created heading towards the centre of the diagram, the anti-soliton away (we 
are being somewhat cavalier with the direction of time; this should be considered as a purely 
geometric argument). The anti-soliton must itself obey our lemmas; if in any further decay 
before it reaches the boundary one of the decay products is heading away from the boundary, 
then there would be no way to close the diagram while conserving momentum at every vertex. 



For a crossed process (figure 3.4 ) the breather is the outermost particle, and is again restricted 



in its decay by our lemmas for the same reason. 

The anti-soliton created by the uncrossed process heads for the boundary with a rapidity 
less than U- and so, by lemma |], may not decay. By the same token, the breather of the 
crossed process cannot decay either so, if the initial soliton is not to form a bound state, the 



only possible alternative processes are figure A. 2 and figure |A,3j. If these are found not to occur 



(for example, if the necessary boundary vertices are not present) then the pole must correspond 
to a bound state for any u < j. 

3.4 The Dirichlet case 
3.4.1 The soliton sector 

The Dirichlet case is exceptional in that topological charge is conserved and so Q± = 0. The 
remaining factors can be rewritten as 



P ± (u) =R (u)l\ 



i=i 



r | + 2ZA ± I + ^ r | + (21 - 2)A T J + ^ 



-I 



Y (I + (21 - 1)A + | + r (I + (2Z - 1)A - i + £ ) 

(u -u)\ , (3.17) 



where Rq(u) is as before and £ = r] = ^P-. Taking tpQ to lie in < </?o < § , £ is in the range 



* tt(A + 1) 

< C < ^ ■ (3.18) 



These factors can again be written in terms of Barnes' multiperiodic functions, as 

^2(tT| + 7 r + u|f,2 7 r)5 2 (tTj-«|f,2 7 r 
P ± ( U ) = floM— 7 ; : ( 7 : f » (3-19) 



with 



S 2 (l-u 


^,2vr)5 2 (^ + u 




S 2 (% + u 


^,2vr)5 2 (^-u 


2T' 2vr ) 



Mu) = S)l ; i A \ • (3.20) 
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3.4.2 The breather sector 

In the Dirichlet case, with topological charge conserved, the bootstrap equation reduces to 

/^S + y) fljft (au)j^ (« - Y) = fh,%) («)• ( 3 - 21 ) 



Ghoshal found that, for the boundary ground state, the breather reflection factors were 

?" 

l '|0) 



R^( U ) = R^(u)RP(u), (3.22) 



where 



v( n )r„.\ _ (j) (2" + 1 ) TT (A) (2A + 1 ) /o 



v2A + i) j=i (2k + I) 

and 



= i§2 ( AlF + 3 + 2A") 

This makes use of the notation 



sinh C - -I- 



which will also be helpful later. 



3.5 Initial pole analysis 

3.5.1 Solitonic ground state factors 

The Ro(u) factor is insensitive to the boundary parameters, and so all its poles should be 
explicable in terms of the bulk. The only poles are at u = ^y, where N = 1,2,3,..., with no 
zeroes. These can be explained by the creation of a breather which is incident perpendicularly 



on the boundary, as shown in figure 3.7. Here, as in all subsequent diagrams, the time axis 
points up the page, and the x axis points to the right. Solitons and anti-solitons are drawn as 
solid lines, while breathers are drawn as dotted lines. 

Turning now to ^-dependent poles and zeroes, we find zeroes at 

£ (2n+l)7r 

U = "A + ^T-' (3 ' 26) 

where n = 0,l,2,..., for P + , and at the same rapidities but with £ — > — £ for P~ . There are 
also poles in P + only at u = u n , with 



£ _ (2n+l)7r 
A 2A 



(3.27) 



A soliton can only decay into an anti-soliton and a breather, with a rapidity difference between 
the two °f f + 2T f° r breather b. Thus, by lemma ^, all these poles must correspond to bound 
states, as shown in figure ED]. For reasons which will become clear in a moment, we shall depart 
from the convention of |29f| and, rather than labelling the state corresponding to pole v n as f3 n , 
will label it according to topological charge and n as 
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3.5.2 Solitonic excited state reflection factors 

Using the boundary bootstrap equations given in [13] — which come from considering figure |3~To| — 
solitonic reflection factors can be calculated for this first set of bound states. In our case, these 
equations read 

K)^ =E J ^>( u ) s rf( w - Q 2x)Sta(«+«L) > ( 3 - 28 ) 

c.d 

where a,b,c, and d take the values + or — and is the (imaginary) rapidity of the pole at 
which particle a binds to boundary state \x) to give state \y). The mass of state \y) — m y — is 
given by 

m y = m x + m s cos a v ax . (3.29) 

Taking x to be the ground state |0) and y to be one of the set of excited states |l;n), this 
gives 

P \i:n)(. u ) = P\Q)(u)a(u - v n )a(u + v n ) 

P \l;n)( U ) = P \Q)( U ) b ( u ~ Vn)b{u + V n ) + P^(u)c(u - V n )c{u + V n ). (3.30) 
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Note that P^. Q Ju) = P^(u), where P±(u) is P^(u) under the transformation £ — > 7r(A+ 1) — £. 



The reason for this is clear if we look back at figure 3.1 ; this transformation is equivalent to 



reflecting the diagram in the horizontal axis, interchanging the ground and first excited states. 
Perhaps the neatest way to write the new reflection factors is 

P± n) (u) = P^)a 1 n (u), (3.31) 



where 



The factor a\(u) simplifies to 



if \ a(u + v n )a{u - v n ) 

a„[u) = —. —. r . (3.32) 

nV ' a(u + v )a{u - u ) y ' 



_L_lJ__£1 (A L_£ 

i , , -i-r \Xtt^~2X X J \ Xtt 2X X i 

4(« = II 77^ fyi ; ( 3 - 33 ) 



l + l_£_Ll) f J J__£ + 1 

Xtt ~ 2X X ' 1 1 \ Xtt 2X X ' 1 



Looking at the pole structure, we find that the functions P ± (u) have common simple poles 
at vq and v-n where N = 1,2, 3, . . . . In addition, P+(n) has simple poles at u = wn>, where 

(3.34) 



£ 7r(2iV'-l) 

w N , = vr - ^ ^ — i = z^at/ 



and iV' = 1, 2, 3, . . . , and simple zeroes at —w^n for appropriate values of N" . Finally, a*(u) 
has simple poles at vq and v n , and double poles at Vk, k = 1, 2, . . . ,n — 1. 

Before proceeding to a more rigorous discussion, we shall now digress to give an outline 
of how the bootstrap might be expected to work. If one of these poles does correspond to 
a new bound state, factorisability leads us to expect that moving the soliton and anti-soliton 
trajectories past each other (so that the anti-soliton is incident on the boundary first) should 
also create the same state. The most obvious explanation for this would be for the anti-soliton 
to bind to the boundary first, followed by the soliton, to form the state. From above, however, 
only solitons can bind to the ground state, so we must look further. 

The next most obvious way this could happen is via the soliton and anti-soliton forming 
a breather, either before or after the anti-soliton has reflected from the boundary. The poles 
required to allow the first process (of the form ir + 4 — ox ) are n °t present, whereas 
those necessary for the second (of the form w m ) are. Our candidate process therefore becomes 



figure 3.11, where the soliton and anti-soliton bind to form a breather, which then creates 
the state in one step. It is quite difficult to imagine any further alternatives, so let us — for 
the moment — take the existence of such a process as a necessary condition for a pole to be 
responsible for the formation of a boundary state. 

The consequence of this is that the wn poles are selected as the only possible candidates, 
and it appears that new bound states can only be formed by anti-solitons. Such states hence 
have charge (agreeing with the idea that they can also be formed from the ground state by 
the action of a breather). In addition, it is also clear that only those wn such that wn < v n can 
be considered, as, otherwise, the breather version of the process would see the breather created 
heading away from the boundary, rather than towards it. 



Designating such a new state as \0;n,N) and bootstrapping on it leads to 
^0;n,JV> ( U ) = P \l;n) ( u ) a ( u ~ W N )a(u + W N ) 

P \0;n.N) ( U ) = P \l;n) ( U )K U ~ W N )b(u + W N ) + P^ (u)c(u - W N )c(u + W N ) . 
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Figure 3.11: States can be created either by breathers or solitons 



Substituting in ( p. 31 ) and taking advantage of the fact that wn = (so a(u db wn) 



a(u ± vn)), this becomes 



P \o-,n,N)( u ) = a li( u ) P \o)( u ) a ( u ~ v N )a{u + v N ) (3.36) 



P \0;n,N)( U ) = a l( P \ )( U ) b ( U ~ "jv) 6 (« + »n) + P^(u)c(u - V N )c(u + U N )) , 



which (apart from an extra factor of a* (u)) is just the first bootstrap ( 3.30| ) under the transfor- 



mation £ — > 7r(A + 1) — £ and with solitons and anti-solitons interchanged on the lhs. Thus, the 
pole structure follows naturally from the above. This can also be written as 



P \0;n,N)( U ) = P^iuKiuy^u). (3.37) 

Repeating the factorisation argument shows that now we should focus on v n i poles such that 
v n i < wj\f. These are present now in the solitonic factor, though (due to the extra factor of 
a^(u)) only for n' > n. However, since any such state obeys v n > wn > v n < in any case, this 
restriction is not relevant. The resultant state must now have charge 0. 

A pattern is emerging, and it is not hard to see how the process would continue. Starting 
from the ground state, and taking the broadest guess (given our assumptions) for the spectrum, 
states can be formed by alternating solitons and anti-solitons, the solitons having rapidity v ni 
and the anti-solitons having rapidity w n . (for some sets n and N_) . An schematic pole structure 



is shown in figure 3.12, in terms of which the criterion for a state to be in the spectrum should 
be that we begin with one of the v n and then, as we move along the index list, move down the 
diagram, switching from side to side as we go. If we finish on a v m (indicating that the most 
recent particle to bind was a soliton) the state has charge 1 while, if we finish on a w m (meaning 
an anti-soliton) the state has charge 0. 

Annotating such a state by its topological charge, c, and the sets n and N_ as |c; n\, N±, ri2, 
N%, . . . ) (noting v ni > wn 1 > v n2 > W N 2 > ■ ■ ■ ) ■> the solitonic reflection factors should be 



p i ni ,N lt ...)( u ) = p tc)( u X( u ) a k( u ) • • • ' ( 3 - 38 ) 



with Pq(u) = P^(u) and Pf(u) = P^(u). From now on, however, it will be more convenient 
to consider a single index list, and denote a^(it) as a^(u), giving 

P± Ju)=P±(u)al(u)a°(u)alJu)...a c „(u), (3.39) 

\c;ni,ri2,.. . ,rifc) v ' (c) v > n\\ ) n 2 \ ) ) rifeV /> v / 
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Figure 3.12: Location of poles. (Note that, in this case, W2 can never participate in bound 
state formation as it is above vq.) 



where k is odd if c is 1 and k is even if c is 0. We will call this a level k boundary bound state. 

,2 fizi' 

2A 



If we choose the ground state mass to be m s sin 2 



, the mass of this state is 
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(3.40) 

(3.41) 
(3.42) 



This choice is convenient in that, as £ passes 7r//3, the masses of the ground and first excited 
states interchange, in line with the idea that the states themselves swap at this point. An 
important point to note is that, in deriving all this, we have simply been considering the soliton 
sector. However, we will see that allowing breather processes as well does not give rise to any 
further states, merely additional ways to jump between states. The Dirichlet boundary condition 
is also special in that either the soliton or the anti-soliton can couple to a given boundary, but 
not both, as might be generically expected. 

Although we have built up the states by applying the solitons and anti-solitons in this 
alternating fashion, precisely how this happens in a given situation will of course depend on the 



impact parameters of the incoming particles. In figure p.ll| we already gave an example of the 
complicated way in which a process may be rearranged as these impact parameters vary, and 
the particular choices that we have adopted are mainly motivated by a desire to assemble the 
full spectrum in the simplest possible way. 



3.5.3 Breather ground state reflection factors 

We now return to the pole analysis, and examine the breather ground state reflection factors 
( 3.22] ). Again, the factor Rq is boundary-independent, and so all its poles should have an 
explanation in terms of the bulk. There are (physical strip) poles at 5 , ^ , f — §f , an d double 
poles at - — ^x, with I = 1, 2, . . . , n — 1. There are no zeroes. The pole at \ is simply due to 
the breather coupling perpendicularly to the boundary, while the poles at |r are explained by 
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Figure 3.13: Breather triangle Figure 3.14: Breather double 
process triangle process 



figure 3.13; . Next, the pole at ^ — Sf comes from a breather version of figure 3.7, i?2n being 



formed. Finally, the double poles at S 



2A 



are due to figure 3.14 



Moving on to the boundary-dependent part, there are poles at 



and zeroes at 



£ 7r lir 
X~2 ± 2X , 



£ tt lir 

U = —X + 2 ± 2X 
£ 7r lir 

U = A + 2 ± 2A' 



(3.43) 



(3.44) 



where, for a breather n, I 



n 



l,n-3,... ,1 > 0. 



The set of poles can be re-written by noting that, for breather m, there is a simple pole of 
the form \{y n — wn) for all n, N > and n,N E Z such that m = n + N . This ties in with 
the discussion in the previous section, since these are the rapidities predicted for the single-step 
process which is equivalent to a soliton binding at an angle of u n followed by an anti-soliton at 

WN- 

3.5.4 Breather excited state reflection factors 



Following the discussion of the solitonic excited state reflection factors, we can introduce corre- 
sponding breather reflection factors: 



- n '|c;ni,n2, 



nk) (u) = R^(u)a^(u)a^(u)a^(u) . . . a^(u 



where Rff(u) = RZ(u) and Rf(u) = R^{u). We have also defined 
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m \ c 
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(3.45) 



(3.46) 
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with an m (u) 



a n (u). 



For i?S (u) , there are poles at 
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and zeroes at 
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(3.49) 



For the other factors, a n ,m (u) has physical strip poles/zeroes at 



S. i ZL _l_ EZE 



1/ 



1/ 



+ 



2A 



p7T 

2A 



i _l ZL _i_ 2ZE 
A 2 2A 



while arl m {u) has them at 



u 



L _i_ 3ZE _i_ EL 
A "•" 2 "T" 2A 



£ 4- ZL _i_ El 
\ ' 2 ' 2X 



£ _ ZL _i_ EZL 

A 2 "r 2A 

£ _l_ ZL _i_ EH 

A ~r 2 "•" 2A 



poles : 


P = 


2n — m + 2x ± 1 


zeroes : 


P = 


— m + 2x db 1 


poles : 


P = 


m — 2x dz 1 


zeroes : 


P = 


— 2n + m — 2x ± 1 


poles : 


P = 


— 2n + m — 2x ± 1 


zeroes : 


P = 


m — 2x db 1 


poles : 


P = 


-2iV-m + 2xd=l 


zeroes : 






poles : 


P = 


— m + 2x dz 1 


zeroes : 


P = 


-2iV - m + 2x ± 1 


poles : 


P = 


2N + m - 2z dz 1 


zeroes : 


P = 


m — 2x zb 1 


poles : 


P = 


m — 2x zb 1 


zeroes : 


P = 


2A^ + m - 2x zb 1 



These poles will be further discussed in section 3.7 below 



(3.50) 



(3.51) 



3.6 An example 



To get an idea of the full picture, and which processes are responsible for the remaining poles, 
we will now look at one particular example in more detail. If we select £ = 1.67T and A = 2.5, 
then we have the first two breathers in the spectrum, with the solitonic poles taking the form 
v n = 7r ( 2 - 2 ~ 2n ) an d Wn = 7r(2.8-2JV) ^ Thus, for this case, only the poles at Vq,Vi and W\ are on 



the physical strip, and so figure 3.12 is simplified to figure 3.15. This is the simplest case which 



requires a broader spectrum than that postulated in [29|. First, let us turn to the soliton sector 
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u 







Figure 3.15: Location of poles in the example 



3.6.1 Boundary ground state — soliton sector 

As argued above, the soliton can bind to the boundary at all rapidities v n which are in the 
physical strip, here just comprising vq and v\. This introduces the states |1; 0) and |1; 1). 



3.6.2 Boundary ground state — breather sector 



The only breather poles are at % — 5 + 



(m—l)n 
2A 



for breather m. In addition, breather B2 has a 



zero at — | + | + 



By lemma [l], the pole for B\ must correspond to a new bound state, the rapidity being less 
than t£-. From figure 3,ll| , it is clear that B\ creates the state which was labelled |5o,i) m [Hj; 
and which we have called |0; 0, 1). 



The pole for the second breather can be explained by figure A. 5 , with the state |1; 0) being 
formed. The anti-soliton is reflected from the boundary at a rapidity of 4 — 7r + 1? — a zero of the 
1 1 ; 0) reflection factor — reducing the diagram to first order through the boundary Coleman-Thun 
mechanism. 



3.6.3 First level excited states — soliton sector 



From before, Pu. Q \ just has a simple pole at i/q, which can be explained by the crossed process in 



figure |3.9| , reducing the boundary to the ground state. For P\t.iy the pole at v\ can be explained 
this way while, for the double pole at z/q, figure A. 4 is required, the first breather being formed 



while the boundary is reduced to the vacuum state. 

For P\i. n \(u), we have the additional job of explaining simple poles at wn, for all N such 
that this pole is in the physical strip. Here, this is only w\. For |1;0), this is appropriate for 
the formation of |0; 0, 1) which, from the previous section, must be present. For |1; 1), however, 



figure [A.3| is invoked, the second breather being created, and the boundary reduced to the 
vacuum state. The breather is incident on the boundary at an angle of \(w\ — v\) = ir — 4 — ^ 
which, looking at the above breather reflection factors, is a zero, ensuring the diagram is of the 
correct order. 
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3.6.4 First level excited states — breather sector 



The pole structure of R™. } can be found from R^y and is 



B\ : pole at ^ 
B 2 : poles at | 



X ' X 

£,37r 7T g i 7T 

A "i" 2A' 2 A ~r 2A 



(3.52) 



By lemma |], the second pole for B 2 must correspond to a new bound state; by the previous 
arguments, this is |1; 0, 1, 1). This state is not in the spectrum given in [29], but lemma |l| shows 
that there is no way to avoid its introduction. Considerations such as this will open the door to 
a much wider spectrum in the general case. 

The B\ pole is suitable for the creation of |1; 1). The first pole for B 2 can be explained by 



figure A. 7, with the boundary being reduced to the ground state by emission of a soliton. 



For Ru.y, the above poles are supplemented by additional poles from b^^u) to give the 



poles shown in table 3.1. 
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Table 3.1: Breather pole structure for |1;1). Entries are the values of p for which there 
is a pole in the location given in the column heading. The power of the entry gives the order of 
the pole, so e.g. 3 2 indicates a double pole when p = 3. There are no physical strip zeroes for 
either breather. 



The pole at 4 + § — fr can be explained by figure A.8| , with the boundary being reduced 
to the ground state by emission of a soliton. The pole at 4 — ~ for B\ can be allocated to 
the creation of |1; 0, 1, 1), while the pole at 4 — | + ^ for B 2 is due to figure A.9| , where the 
boundary emits B\, being reduced to |1;0). The pole at — 4 + | + |? for B\ is responsible 
for this reduction to 1 1 ; 0) , while the double pole for B 2 comes from an all-breather version of 
figure A. 6, the boundary being reduced in the same way. 



3.6.5 Second level excited states — soliton sector 

For -Pi frQj ) (u), the only poles are simple, at u and W\. The pole at w\ can be explained by 



figure |3.9| while, for vq, we need figure A. 2. The second breather is emitted by the boundary, 
reducing it to the ground state, while a soliton is incident on the boundary at a rapidity w\. 
For the ground state, this is neither a pole nor a zero, but the diagram contains a solitonic loop 
which can either be drawn to leave a soliton or an anti-soliton incident on the boundary. Adding 
the contributions of these two diagrams gives an additional zero. 

For P^.q iAu), we have additional poles at all u, i.e. a simple pole at u\, with vq becoming 
a double pole. By lemma [l], v\ must correspond to the creation of a new bound state, namely 



|1; 0, 1, 1), while, for vq, figure A. 3 should be considered. Again, the second breather is created, 
the boundary is reduced to the ground state, and the breather is incident on the boundary at a 
rapidity of \{vq — w\) = £/A — ir/2 — a zero of the reflection factor. 



3.6.6 Second level excited states — breather sector 



For |0; 0, 1), we have the pole structure given in table 3.2 
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Table 3.2: Breather pole structure for |0;0, 1). 



The poles at — t + + |r are due to figure |A.8| , while the poles in the second column are 
due to figure A.9 . For all these, the boundary is reduced to |1; 0). The pole at 4 — j + for 
2) is due to figure A.12| , while for 
- - + ~ 



S m (m 
The pole at 



m 



1 it is due to a breather version of figure 3.£. 



2tt 
2A 



for £?i is due to figure A. 7. 



3.6.7 Third level excited states — soliton sector 

The only third level excited state is |1; 0,1,1). For P^. 011 y there are simple poles at Wi,vq 
and v\ . Again, the pole at w\ comes from the crossed process figure [T{^. For ui, figure 
suffices while, for uq, figure |A.4| is required, the boundary being reduced to |0;0, 1) while the 
first breather is incident on the boundary at § — 4 + f> another zero. 



3.6.8 Third level excited states — breather sector 



Here, the only possible boundary state is |1; 0, 1, 1) and we find the poles given in table 3.3 
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Table 3.3: Breather pole structure for |1; 0,1,1). 



Comparing this with the structure given above for |1; 1), it can easily be seen that, whenever 
the two both have a pole at the same rapidity, essentially the same explanation can be used. For 
the remaining poles, — | + 4^ + |? can be explained by figure |A.9| , the boundary being reduced 
to 1 1; 0) , while that at — ^ + | + ^ for B 2 is due to figure reducing the boundary to |1; 0), 
and that at 4 — § + |r for B\ is due to an all-breather version of figure A. 7, again reducing the 
boundary to |1; 0). 



3.6.9 Summary 

The above has shown that, by introducing only the states which are required by lemmas 1 and 2, 
the complete pole structure can be explained. Below, we shall find that this is a general feature. 
In addition, the spectrum of states is broader than that introduced in |^] (containing, in addition 
to their states, |1; 0, 1, 1)). It should be noted that the mass of this extra state corresponds to 



mix °f [29] , the mass of a boundary Bethe ansatz (l,l)-string whose apparent absence from 
the bootstrap spectrum was described in that paper as "confusing". This does at least show 
that the Bethe ansatz results of [^] are not incompatible with the bootstrap. However, in more 
general cases it turns out that the bootstrap predicts yet further states, beyond those identified 
in the boundary Bethe ansatz calculations of [^], so a full reconciliation of the Bethe ansatz 
and bootstrap approaches remains an open problem. 
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3.7 The general case 

From the above, we might imagine that the boundary state |c; ni,n,2, TI3, . . . , n m ) exists iff c is 
or 1 and n\, n 2 , n 3 , . . . are chosen such that ir/2 > u ni > w n2 > v nz > . . . > 0. This turns out to 
be correct, and will be proved in two stages. Firstly, we need to show that all these states must 
be present, before going on to show that, given this, all other poles can be explained without 
invoking further boundary states. 



3.7.1 The minimal spectrum 

The argument proceeds as follows: starting with the knowledge that the vacuum state |0) and 
all appropriate states |1; n\) are in the spectrum, we use breather poles to construct all the other 
postulated states. 

These poles are of the form \{wn — u n ) for breather n + N incident on a charge state (or 
\(y n — wn) for a charge 1 state). If v n — wn < j, lemma [l] shows that they must correspond 
either to the formation of a new state, or the crossed process. From figure |3.11 , this corresponds 



either to adding indices n and iV if they are absent or — if they are already present — removing 
them. (Note that any other option would give rise to a state with a mass outside the scheme 
given by ( |3.40| ), and therefore outside our postulated spectrum.) The condition v n — wn < j 
is always satisfied if v n > wn and v n and wn are as close together as possible, i.e. if |0;n, N) 
exists, but \0;n,N — 1) does not. 

The only subtlety in this argument arises when considering the topmost breather. If n + N = 
n max, lemma [l] on its own is not strong enough to require the presence of the state we need, 
and we must invoke the stronger version introduced at the end of of section |3.3| . This makes 
use the idea that there must be a corresponding two-stage solitonic route to the same state, i.e. 
a soliton with rapidity u n followed by an anti-soliton with rapidity wn- Considering these two 
processes instead, the stronger lemma shows that both form bound states, as v n and wn must 
be the lowest poles of their type — and so have rapidity less than j — for n + N to equal n max . 
This shows that the state exists, and hence that the breather pole is due to its formation. 

Since the arguments for the two sectors are analogous, let us focus on the charge sector here. 
The challenge is to create any state |0;x) — for some set of indices x = (ni,ri2, ■ ■ ■ ,n 2 k) — from 
the ground state using just these poles. As a first step, consider creating |0; 771,772)- If v ni and 
w<n2 are as close together as possible, we simply make use of the pole at ^(w^ ~ v n 1 )- Otherwise, 
introduce the set mi,m2, ... ,mt such that u ni > w mi > z/ m2 > w mz > . . . > v mt > w n2 , with 
each successive rapidity as close to the previous one as possible. Now, we can successively create 
|0; x, ni, mi), then |0; x, n±, mi, 1712,1713} and so on, up to |0; x, n\, mi, 1172,1713, . . . , m t , 112). 

By now invoking the crossed process, a suitable breather can be used to removed the indices 
mi,m2, followed by 7713,7714 and so on, until all the m indices have been removed to leave 
\0;x, 77,1,77,2). 

Repeating this procedure allows |0; ni, 172, 773, 774) to be created, and hence |0;x). Note that 
this allows any state in our allowed spectrum to be created, but no others, as the condition 
v ni > w n2 > . . . is imposed by the existence of the necessary breather poles. Charge 1 states 
can be created analogously by starting from a suitable state |l;77i). 

One remaining point is to check that all the necessary breather poles do indeed exist. How- 
ever, starting from ( |3.45 ), they occur in the R^(u) factor, and it is straightforward to check 



that they are never modified by the other a factors. 
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3.7.2 Reflection factors for the minimal spectrum 



The boundary state can be changed by the solitonic processes given in table 3.4 



Initial state 


Particle 


Rapidity 


Final state 


10;%,. . . ,n 2k ) 


Soliton 




\l;ni, . . . ,n 2k ,n) 


|l;m, . . . ,n 2k -i) 


Anti — soliton 


W N 


|0;m, . . . ,n 2 fc-i,iV) 



Table 3.4: Solitonic processes which change the boundary state. 



The breather sector is more complex, as indices can be added or removed from any point in 
the list, and not just at the end, as for solitons. In addition, processes exist which simply adjust 
the value of one of the indices, rather than increasing the number of indices. For breather m, 



these are given in table 3.5. This should be read as implying that any index can have its value 



Initial state 


Rapidity 


Final state 


1 0/1; . . . n 2x ,n 2x +i ■■■) 
|0/1; . . . n 2x -i,n 2x ...) 
\0/l;...n 2x ...) 
|0/1; . . . n 2x -i ■■■) 


\(y n - WN),n + N = m 
^(wn - v n ),n + N = m 

2 { v -n 2x ~ W rl2x J rm ) 
2 {W-n2x-l — W n2x _ 1+rn ) 


|0/1;... n 2x ,n,N, n 2x+1 ...) 
|0/1; . . . n 2x -i,N,n,n 2x ...) 
\0/l;...n 2x + m...) 
|0/1; . . . n 2x -i +m...} 



Table 3.5: Breather processes which change the boundary state. 



raised, and that a pair of indices can be inserted at any point in the list, including before the 
first index and after the last (providing the resultant state is allowed). Both these tables have 
been derived on the basis that, whenever assuming that a pole corresponds to a bound state 
leads to a state with the same mass and topological charge as one in our minimal spectrum, the 
assumption is taken to be correct. As with our earlier assumption (that, if a pole has another 
possible explanation, it is not taken as forming a bound state), this is intuitively reasonable but 
not necessarily rigorous. It does, however, lead to consistent results, and there is no conflict 
between the two assumptions: we have been unable to find any alternative explanation for any 
of the poles listed above. 

It is vital for what follows that, for all the above processes, there is very little dependence on 
the existing boundary state. For the solitons, the topological charge of the state and the value 
of the last index in the list are all that matter. Any two states which have the same topological 
charge and last index can undergo processes at the same rapidities to add an index. Similarly, 
for the breathers, provided either the relevant two indices can be added at some point in the 
list to create an allowed state, or that the index to be adjusted is present in the list, the other 
characteristics of the state are irrelevant. 



3.7.3 Solitonic pole structure 



This turns out to be relatively straightforward. All poles are either of the form u n or wn- 
Looking at a charge state with 2k indices, and labelling this as x = (ni,n 2 , . . . , n 2 k), we find 
the results shown in table 3.6 for P^^(u). These poles come from the a factors so, for P + , 
there is an additional pole at all v. 

For the charge 1 states, the picture is very similar, and, considering P + first, we find the 
pattern given in table 3.7 for a state with 2k — 1 indices. For P~ there are additional poles at all 
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Pole 


Order 


Pole 


w 1 ... w n2 _i 


2k 


V x ..-Vnx-\ 


W n2 


2k - 1 




W n2 +1 ■ ■ ■ W n4 -1 


2k -2 


V ni +l . . . 1^3 — 1 




2k -3 






2 


u n 2k _ 3 +l • • • ^nak-i-l 


Wn 2k 


1 





Table 3.6: Pole structure for P, Q , (u). j4n eniry o/ ; /or example, w\ . 
i/icrf i/iere zs a po/e o/ order 2fc toi, U72, ^3, • • • , w n2 -i ■ 



itf n2 -l indicates 



w. (For the charge case, there are poles at w x for x < 0, but none of these are in the physical 
strip.) 



Pole 


Order 


Pole 




1 


VQ,V-U . . . 




2k 


v\... v ni -l 




2k -1 


Vm 


W\ . ..W n2 -i 


2k -2 


u ni+l • ■ ■ u n 3 -l 


W n2 


2k -3 


v n 3 


w n 2k -4+l ■ ■■ w n 2k - 2 -l 


2 


u n 2 k-3+l ■ ■ ■ h 'n 2k _ 1 -l 


W n 2 k-2 


1 


u n 2 k-i 



Table 3.7: Pole structure for P^.^u). 

An important point to note is that, comparing |0;rai,?i2, . . . , ?i2fc-i> n 2fc) ( a general level 2k 
state) with the level 2 state |0; ^2fc-i) n 2k), we find no additional poles, though the order of some 
poles has increased. In the example above, all level 2 states were explained by diagrams where 
the boundary was reduced either to the vacuum by emission of a breather, or to a first level 
excited state by emission of an anti-soliton. The same processes turn out to be present for any 
level 2k state to be reduced to a level 2k — 1 or 2k — 2 state. Thus, we might imagine explaining 
the poles in the level 2k reflection factor via similar processes to the ones which explained them 
in the level 2 factor. At times, however — as we shall see — parts of these processes will need 
to be replaced by more complex subdiagrams to allow for the fact that the boundary is in a 
higher excited state, explaining the differences in the orders of the poles. Considering the level 2 
processes so far introduced as "building blocks" , this can be considered as an iterative process: 
level 4 states can be explained by replacing parts of level 2 processes with building blocks, while 
level 6 states can be explained by similarly replacing parts of level 4 processes with building 
blocks, and so on. A generic process of the type we will examine can therefore be viewed as a 
cascade of building blocks, each appearing as a subdiagram of the one before it. 

A similar argument applies to level 2k+l states and level 3 states, drawing the same diagrams 
with all rapidities transformed via £ — ► ir(X + 1) — £. We will concentrate on the charge sector 
below, and consider a generic level 2k state. 



For poles of the form u n , consider figure A.13j . The boundary decays to the state |0; ni,n,2, • • • , 
n 2k~2) by emission of breather ri2k + "2fc-l a t a rapidity of \{vn 2k _ x — w n2k ). This then decays 
into breather n-ik-l — n heading towards the boundary at a rapidity of ^(w- ri2k _ 1 — v n ) and 
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breather ri2k + n heading away from the boundary at a rapidity of v n — (j— ^ TO2fc +") 7r ) . This 



2 2A 

then decays to give the outgoing particle and one heading towards the boundary at a rapidity of 
w n2k . For n < n<2k~i, it is straightforward to check that all these rapidities are within suitable 
bounds. 

This diagram is naively third order. However, breather ri2k-i — n , which is drawn as simply 
reflecting off the boundary, in fact has a pole, meaning that the diagram should be treated as 
schematic and the appropriate diagram from the next section inserted instead. In addition, 
as noted in the discussion of the example, the soliton loop contributes a zero for an incoming 
anti-soliton through the Coleman-Thun mechanism. When this is taken into account, we obtain 
the correct result. 



For V2k-i, the slightly simpler figure A. 2 suffices. The remaining v poles are only present in 



the soliton reflection factor, and can be explained by figure A. 3, with the boundary decaying by 
emitting an anti-soliton at w n2k , which then interacts with the incoming soliton to give breather 
n + n2k, heading towards the boundary at a rapidity of \(y n — w n2k ). Looking ahead again, the 
interaction of this breather with the boundary contributes the required zero. For v n < w n2k , 
this diagram fails, the breather being created heading away from the boundary; this is the point 
when the pole is to be considered as creating the bound state |1; n±, . . . , ri2k,n). 

For the wn poles, the story is very similar, this time being based on figure [A.4| (requiring a 
suitable pole/zero for B^^ n2k on state |l;rti, . . . , n,2k-i) at 4 — § + 7r ( Ar+ '^fc-i) ) for N < ri2k 
and figure |3^ for n 2 fc. As noted above, all charge 1 state poles can be explained by the same 
mechanisms, with the rapidities transformed according to £ — ► 7r(A + !) — £. 



3.7.4 Breather pole structure 

This is considerably more complicated. However, with a bit of work it turns out that, for 
breather n on the state |0; n%, ri2, ■ ■ ■ , n^k), the pole structure is as given in table 3 



In explaining all this, we can begin with the diagrams found previously. For the first line, 



consider an all-breather version of figure A. 5, where the breather decay is chosen to produce 
breather n + x — ri2q' on the left, which then binds to the boundary to raise index n^q' tow + i. In 
some cases, this is not possible, the appropriate state not being in the spectrum, but, then, we 



can consider an all-breather version of figure A. 10, where the boundary decays so as to remove 



the indices U2 q ' and r^'+i, with the same initial breather decay, and the additional breather 
reflecting from the boundary contributes a zero. This diagram becomes possible just as the 
other fails. In either case, the other breather from the initial decay (which is drawn as simply 
reflecting from the boundary), is breather y = n 2g ' — x at rapidity 4 — | + ZliJd^LJi. This has 
a pole of order 2 less than the initial breather. If this order is less than or equal to zero, the 
diagram stands as drawn while, otherwise, the simple reflection from the boundary should be 
replaced by a repeat of this argument, iterating until the result is less than or equal to zero. For 
the next line, precisely the same argument can be used. 

The next three lines can be explained by a similar argument, based on either increasing the 
value of index n,2 q '-i or removing indices n2q'-x an d n^n' ■ 

For | + § + l 

r(n+2z i) ^ we invoke a similar process. This time, however, the outer legs have 
rapidity vi nJrX \ (where —(n + x) is actually a positive number if the initial pole is to be in the 
physical strip), and so we need to substitute in the explanation of soliton poles of this form from 



before, leading, in simple cases, to figure A. 9 



Finally, for — 4 + ^""oa" > we begin with figure A.J. This time, the reflection factor 



for the central soliton always provides a zero, while the outer soliton has rapidity w n+x . If 
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Pole 



I 7r | n(n+2x-l) 
A 2 2A 



£ , TT , 7r(n+2z+l) 
" A 2 



2A 



£ , 7T I 7r(n+2a;-l) 
A 2 2A 



£ , 3tt _i_ tt(?i+2z+1) 
A 2 2A 



Range 



n 2q < x < n 2q+2 
n 2q ' < n + x < n 2q '+2 
x < 0,n 2<? _i < \x\ < n 2q+1 
n 2q i < n — \x\ < n 2q i +2 

n 2q -i < x < n 2q+ i 
n 2q i_i < n + x < n 2q i + i 
x < 0, n 2q < \x\ < n 2q+2 , 
n 2q '-i < n — \x\ < n 2q t + i 
x < -n,n 2q < \x\ < n 2q+2 
n 2q ' <\x\ - n < n 2q > +2 

a f _ 7T , n(n+2x+l) 
A 2 "r 2A 

with poles <-> zeroes 
As -!■ + - + 7r ( n+2a; ~ 1 ) 

A 2 2A 

with poles <-> zeroes 



Pole/zero order 



2(</ - q) + y 
2(g' - g) + y + i 

2{q' -q) + y 
2{q' -q)-i + y 
2(q'-q) 



Table 3.8: Breather pole structure for a generic charge state. The variable x takes 
integer and half-integer values within the allowed ranges. An entry in the third column represents 
a pole of that order if it is positive, and a zero of appropriate order if it is negative. (Thus an 
entry of +1 is a first-order pole, and an entry of -1 is a first-order zero.) Also, for convenience, 
y is 1 if x (or \x\) attains the lower limit, -1 if n + x (or n — \x\) attains the lower limit, and 
zero otherwise, while i is 1 if x is integer, and otherwise. 



n + x = n 2 k, the diagram is as drawn while, otherwise, we need to replace the two outer anti- 
soliton legs with the explanation of the appropriate pole in the anti-soliton factor. The first 
iteration of this is shown in figure [A.1C . 



3.8 Number of states 

In this section, we examine how the size and content of the boundary spectrum changes with 
variation in £ and A. Since any state can be formed by a suitable sequence of solitons and 
anti-solitons, we will focus on the solitonic sector. 

The relevant poles, v n and w n , both have the same spacing — ? — but, interestingly, the range 
of n for which v n is in the physical strip is independent of A, while that for the u>-type poles 
is not. For v n , < n < | — | while, for w n /, A — | — ^<n'<-| — § + §• Designating the 
lowest-rapidity ^-type pole as i/„, , there are n* + 1 z/-type poles, and either n* or ra* + 1 relevant 
w-type poles, depending on whether the lowest-rapidity pole is z^-type or w-type. (Note that 
any u;-type pole with a rapidity greater than can never be relevant in forming a bound state.) 

Recall now that the criterion for a state |c; n±, n 2 , n^, . . . } to be in the spectrum is that 
u ni > w n2 > u n3 > . . . , corresponding to moving down figure 3.12 , alternating from side to side. 



Since movement must be strictly downward, there are two cases to consider: when the w and v 
poles occur at the same rapidities, and when they do not. 

The first case is the simplest to deal with, as enumerating the states in the spectrum becomes 
equivalent to calculating the number of ways of making an ordered selection of an arbitrary 
number of objects from a set of n* + 1. However, to simplify the rest of the argument, we shall 
formulate it as a recursion relation. 

We shall consider the situation where w x = v x (realised when £ = 7r ( A + 1 ) ). Clearly, all other 
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cases are similar, with the even indices uniformly increased by | — but with the overall 
spectrum size unchanged. 

Consider first a subset of the spectrum, with all indices less than, say, m, leaving m poles 
to play with in each sector. Denote the number of charge and 1 states in this part of the 
spectrum as co(m) and ci(m) respectively. Now consider extending this to m + 1 poles; all the 
states previously present are still there, together with new states involving the extra index. For 
each sector, a new state can be formed by taking an existing state in the opposite sector and 
adding the new index, m (provided the vacuum state is included in the list of charge states to 
allow for the possibility of forming |1; m)). 

Overall, then, c\(m + 1) = ci(m) + co(m) = co(m + 1). Solving this gives co/i(n* + 1) = 
2 n *ci(0) + 2 ra *co(0). Without allowing any poles, the spectrum consists of only the ground 
state, so co(0) = 1 and ci(0) = 0. Thus Co/i(n* + 1) = 2 n *, as expected from the combinatoric 
approach. 

Moving to the case where the w and v poles do not coincide, the argument changes a little. 
Consider the case where w x lies between v x and v x +i, noting that, as before, all other cases 
simply involve a uniform adjustment of the even indices. Again, we can look at the subset with 
all indices less than m, and compare it with that with all indices less than m + 1. The difference 
now is that we can potentially add two extra indices to an existing state, one from each sector, 
since their rapidities no longer coincide. 

A new charge 1 state can only be formed by the addition of the index m to an existing state, 
but a charge state can either be formed by adding m to an existing charge 1 state, or m, m to 
an existing charge state. Thus, ci(m + 1) = co(m) + ci(m), but co(m + 1) = 2co(m) + c\(m). 
To solve these, it is useful to think of writing out the list cn(0), cj.(1), cq(1), ci(2), cn(2), . . . and 
note that the relation for co(m + 1) can be rewritten as co(m + 1) = c\(rn + 1) + co(m). These 
relations then demand that each element of the list is the sum of the previous two. Since 
cq(0) = ci(l) = 1, this is just a Fibonacci sequence, and we can take advantage of the standard 
formula for the n th term of a Fibonacci sequence, F(n): 

F n = -j= (3.53) 

where (ft is the so-called "golden ratio" (ft = 1+ 2 V ^ . From this, co(m) = F(2m + 1) and ci(m) = 
F(2m). 

One small complication is that, once m = n*, w nt is not necessarily in the physical strip. 
This means that, while the total number of charge 1 states must be ci(n* + 1), the number of 
charge states will either be cn(n* + 1) or cn(n*) depending on whether or not w Ut is present. 
It is perhaps easiest to note that, with n'^ + 1 relevant w-type poles, the number of charge 
states is co(n'^ + 1). 

A plot of the spectrum size against A and jry^ is shown in figure [3.16 , Three sets of curves 



are shown: the points where given v and w poles enter the physical strip, and the points where 
the two sets of poles coincide. As drawn, a given v pole will be in the spectrum above the 
appropriate line, and a given w pole will be present to the right of its line. Note, however, that, 
while v poles will never subsequently leave the spectrum, w poles will; crossing a coincidence line 
to the right, the relevant w pole with the highest rapidity passes vq and ceases to be relevant, 
reducing the number of relevant w-type poles by 1. The number of states in each sector has 
been quoted in terms of Fibonacci numbers, so that a x,y" implies a charge sector of size F x 
and a charge 1 sector of size F y . On the equality line, of course, each sector has size 2 n * +1 . 



This is an integer when the two sets of poles occur at the same rapidities 

Since £ lies between and ,r( - A 2 +1 - > , it is more convenient to work with x+T' w hich lies between and ^. 
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Figure 3.16: Boundary bound state spectrum size. The number of states present in- 
creases as v-type and w-type poles enter the physical strip, but changes also occur as the two 
sets of poles pass through coincidence: moving in the direction of increasing X, the topmost 
relevant w-type pole passes and ceases to be relevant, reducing the spectrum. (Notation x,y 
implies F{x) charge states and F(y) charge 1 states.) 



Finally, note that the top of the diagram represents £ = — -, i.e. the coincidence case 
considered above, and the region just below this represents the other case considered. Moving 
diagonally down and to the right from there, the even indices receive successively greater uniform 
additions but the spectrum size merely oscillates, as the z^-type and u>-type poles take turns at 
having the lowest rapidity. 



3.9 Other boundary conditions 

Surprisingly, the extension of the Dirichlet results to encompass more general boundary condi- 
tions does not require much more work. A hint as to why can be gained from the fact that the 
general ground state reflection factors can be rewritten in terms of those for Dirichlet multiplied 
by terms which introduce no new poles at purely imaginary rapidities. (They do, however, in- 
troduce poles at complex rapidities, but these have an interpretation as resonances rather than 
bound states and will be discussed separately.) 

Despite the fact that the reflection factors appear to depend on four parameters — £,,k,r] and 
& — it is clear that essentially only two are independent, the other two being determined by ( |3.13| ). 
If we note that ( 3.13] ) also implies 



(-l) a 

sin(r/) sin(i?) = — sin£, (3.54) 
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with a either or 1, we can re- write cos(^ + Xu) as 



cos(£ + Xu) = cos(£) cos(Au) — sin(£) sin(Au) 

= —k [cos(r/) cosh(-#) cos(Ait) — (— l) a sin(r/) sinh($) sin(Au)] 
_k 
~2 



(3.55) 



cos (77 + (-l) a An) + e - * cos (77 - (-l) a Xu) 



Denoting the reflection factors for the Dirichlet boundary condition on the vacuum boundary 
state as Pfj, Ju,£), (|3.7| ) can be re- written as 



P ± (u) 
Q ± (u) 



Rq(u 



2 cosh ■& 



e ± ^P^ 0) (u, V )+e^-^P i 



D\0) 



-Ro(u 



sm(\u)a(ii9 ,u) 



0l (i J2cosh(i9)cos(r;) 



D\0)' 



(3.56) 



Since the transformations $ — ► a — ► a + 1, and 7/ — ► —77 are all equivalent to soliton — ► 
anti — soliton, we shall set a = and 1? > for simplicity. The Dirichlet case corresponds to 
1? — ► 00, in which case 77 — > ^ and we recover the expected factors. 

In this form, it is clear that we will be able to re-use much of what we have already found 
about the Dirichlet pole structure in the general case. The one important difference is the factor 
of a(i~&,u). This only has poles at complex u, however, and so will not contribute to the bound 
state structure. We can thus ignore this factor for the present. 

All the reflection factors have the same pole structure at purely imaginary rapidities as P^\ \, 
though based on rj rather than £. Arguing as before, these must be responsible for the formation 
of a first set of excited states. We will continue to use the notation v n to label these poles, on 
the understanding that it is more generally defined as 



77 vr(2n + l) 
A 2X 



(3.57) 



Unlike the Dirichlet case, however, where these poles appeared only in one reflection factor, 
they now appear in all four. While time-reversal symmetry argues that the poles in Q + and 
Q~ must form the same state, we must now deal with the possibility that those in P~ and P + 
potentially form different states, degenerate in mass. This cannot be so, however, since e.g. a 
soliton, incident on the boundary, cannot yet "know" whether it will ultimately be reflected 
back as a soliton or an anti-soliton, meaning that the states formed by Q + and P + must be 
the same. A similar argument holds for anti-solitons, and so for Q~ and P~. Since the states 
formed by Q + and Q~ must be the same, all four states must, in fact, be the same state. This 
also means, for the solitonic sector at least, that all states must be non-degenerate. 

This degree of similarity with the Dirichlet case makes it a reasonable guess that the entire 
structure should also be similar, with reflection factors given by 



cr(jl9,lt) 



(u) 



^0(^)2^10 
. sin(A 
2 cosh(i9) cos(r) c ) 



e ±(-l)^p+^ (Ui??)+e T(-l)^p. 



-R {u)- r ■■ - 



D 



P i\x)( u ^) + P Dlx) {u,v) 



x) 



(3.58) 



where t]q = rj and 771 = 7r(A + 1) — 77. The breather factors, in turn, should be given by 

R$(u) = /4 B) (u)li}^(i7,«)^(t*i) , (3.59) 



L lD\x)' 



L lD\x)' 



where P^i£ ) ] {X ^{ u ) is the boundary-dependent part of the Dirichlet factor (i.e. without the Rq 
term) . 



u 
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One difficulty that might be raised with this idea is that, since topological charge is not in 
general conserved, the two charge sectors might not translate into the general case. As we shall 
see in a moment, the above reflection factors are correct as long as the bound state poles at each 
level match the Dirichlet results. The argument given before for deciding whether or not a pole 
is due to a bound state works just as well here, indicating that this is indeed the case, so the 
conclusion must be that there are still two sectors. In one sector solitons bind with rapidities 
u n and in the other they bind at w n , as is necessary for continuity with the Dirichlet limit. The 
only difference is that the sector label now does not correspond to topological charge; in fact, 
it does not appear to correspond to anything other than the number of labels the state carries. 
For this reason we shall now call them "odd" and "even", rather than charge 1 or 0. The other 
difference is that, at all stages, the poles appear in all four factors, allowing either a soliton or 
an anti-soliton to form a bound state. 

The formula ( |3.58| ) is most easily proven by induction. Since we already know it is true for 
the ground state, all that remains is to show that it is consistent with the bootstrap. 

In its full glory, the boundary bootstrap equation reads 

= £^{6 - ^af 5 )<I(e)5S> + ia&). (3.60) 

Given that we are taking all boundary states to be non-degenerate, and assuming that all states 
can be created by either a soliton or an anti-soliton, we are free to take the incident particle to 
be whatever we please. For convenience, then, we shall set a = 6, leading to 

w = &ss&v - <) R m mf> + <)■ (3-6i) 

(Making the other choice — a ^ b — can be shown to produce an equivalent set of bootstrap 
equations, reinforcing the idea that all reflection factors produce the same boundary state.) 

The boundary couplings g x& can be found from 

0+0 P-0 0+0 

p+ i 9+5^5 p— ^ i 9_s9 s f ->.± ^ i g±$9 s (o m\ 



This means that, using our assumed form for the reflection factors 

^ = <7-(-l) n e ( - 1)C * 

for pole u n or w n as applicable. 

Overall, then, the bootstrap reads 



(3.63) 



a(u - 4) P±(u)a(u + 4) + (-1)^-^0* (u)c(ti + of) 



Qf 3) {u) = a{u-4)b{u + 4)Qf 5) {u). 



Applying this to a state of our assumed form does indeed give (after some cumbersome 
algebra) the requisite result. The other point which remains is to show that, at each step, 
the spectrum is the same as before. However, looking at the breather factors given above, it 
is clear that their pole structure at imaginary rapidities is always the same as for Dirichlet. 
The argument to determine the states which are required in the model depends exclusively 
on breather poles, and so must go through precisely unchanged. The only danger is that the 
remaining enumeration of the explanations for the other poles might run into problems. 

The solitonic factors have poles whenever either of the Dirichlet factors do, the order being 
the higher of the two. Similarly, there are zeroes whenever both Dirichlet factors have zeroes, 
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the order being the lower of the two. This turns out to mean that the explanations used before 
still apply, with the difference that the extra boundary vertices allow solitons and anti-solitons 
to be interchanged within the diagrams in ways not possible in the Dirichlet case. 

This allows a diagram which previously explained a soliton pole to be re-used to explain an 
anti-soliton pole at the same rapidity. In addition, the difference in the order of a pole between 
the soliton and anti-soliton factors was due to loops which allowed a cancellation between dia- 
grams for one but not the other (as in e.g. figure |A~2l ) . Altering the factors from their Dirichlet 
values destroys this delicate cancellation, raising the order to the higher of the pair. With this 
borne in mind, the discussion is completely analogous to that given previously, and so we shall 
not repeat it here. 

Finally, it is also worth noting that the general factors can still be written in the form 

P tcm,Ni,... > («) = P tc) ("X (")•■■ ' t 3 ' 64 ) 
with P^(u) = P^Ju) and P^(u) = P, Q Ju). An analogous expression holds for the Qs. 



3.9.1 Resonance states 

We now return to the extra factor of ct(m?, u). This provides extra complex poles, found from 
the imaginary poles we have been discussing by replacing rj with i$. Thus, the most notable 
poles (and the ones we shall concentrate on) are at k = y - ^jr~^ ■ 

A feature of these poles is that they never fall into the physical strip. Those which fall into 
the unphysical strip immediately below the physical one (as the poles just given do) however, 
do have an explanation as resonance states ||. In bulk QFT, a resonance state is an unstable 



bound state, and a similar idea applies here. From the Breit-Wigner formula [34|, we can find 

i£ 

2 



the mass M and decay width T of the state using the usual formulae with M — > M + For 



the bulk, this becomes 



iT\ 2 

M + — J = mi + m,2 + 2miwi2 cosh(<r — iQ) , (3.65) 



for the binding of particles with masses mi and mi. In our case, we find 



M + — =m s cosh(ci - »0) , (3.66) 



or 



M = m cosh a cos (3.67) 
r = — 2m s sinhcrsin© . (3.68) 

The lifetime of such a particle is r oc j?; to compare this with the discussion in the previous 
chapter, this can be converted into a phase delay by multiplying it by the real velocity v = tanh a 
to get a oc — 2m s cosher sin 0)" 1 . For the poles ^ — o\ ' then becomes 

a oc \2m s cosh ^ sin ^ J . (3.69) 



A 2A 

In the classical limit — ► 0, taking m s = Jj, these become simply 



F oc (3.70) 
a oc (271 + 1)- 1 . (3.71) 
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This means that, in this limit, the resonance states become stable, though the phase delay 
remains finite. The poles collapse onto the real axis, though at an infinite distance from the 
origin. In the classical calculations of the previous chapter, however, due to the re-scaling of the 
field, the poles collapse at a finite distance from the origin with an infinite phase delay, as we 
have already found. 



3.10 Summary 



"No doubt aardvarks think that their offspring are beautiful too." 

— John Ellis 



We have found that the spectrum of boundary bound states of the boundary sine-Gordon model 
can be characterised in terms of two "sectors" . With Dirichlet boundary conditions, these have 
topological charges 4^ and 1 — 4^ (which we labelled as "0" and "1" respectively). Otherwise, 
if topological charge is not conserved, the sectors remain, but lose this interpretation. It is still 
useful to label them as "0" and "1" , but this is best thought of as "even" and "odd" , since they 
require even and odd numbers of solitonic particles for their creation. 

A boundary state can be described in an index notation as \c;m,ri2, ■ ■ ■ ,rik) for sector c, 
with c = for k even and c = 1 for k odd. For the Dirichlet case, such a state can be created 
by a succession of alternating solitons and anti-solitons, beginning with a soliton. With other 
boundary conditions, this requirement is eased, and any selection of solitonic particles becomes 
possible. To create a state in the odd sector, the necessary rapidities are of the form, 

while for the even sector they are 

Wm = «-l- < 2m - l \ (3.73) 

These are interchanged by the transform rj — > 7r(A + 1) — rj. Any such state can be formed, 
provided the rapidities involved are monotonically decreasing, i.e. v ni > w n2 > v n . A > . . . , and 
its mass is given by 



m nun2! ... = m s sm 



2 / V 



2 



2A 



J + Yl m s cos ( u nJ + Yl m sCos{w nj ) (3.74) 



odd j even 

m, sm 



] even 



odd 

rj {2rij — l)ir 
msCOS{ A + 2A 



This spectrum is considerably larger than that suggested in [29], though all the states introduced 
are required to satisfy our lemmas. It is worth pointing out that a second part of the analysis 
of |2l| involved an examination of the (boundary) Bethe ansatz for a lattice regularisation of 
the model. Some of the masses which emerged in the course of that study — those of the (n, N)- 
strings — were outwith the spectrum proposed in [29], but are now included as the masses of the 



states |1;0, n, N). It remains to be seen, however, whether the other masses in our spectrum 
can be recovered in the Bethe ansatz approach. 
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The number of states present in the spectrum clearly depends on the boundary parameters, 
as illustrated in figure |3.16| . It is convenient to express this in terms of Fibonacci numbers, F(x). 
If there are n v-type poles, and m relevant to-type poles, there are, in general, F(2n) charge 1 
states and F(2m + 1) charge states. Explicitly, these are given by 



n 



+ 1 and m 



77 1 

i + 2 



A 



77 1 

i + 2 



(3.76) 



where the square brackets denote the integer part of the number. This changes when the two 
sets of poles coincide, in which case there are 2 n_1 states in each sector. 

Finally, we note that the general method used to derive the spectrum, via the simple geo- 
metrical argument leading to the two lemmas given in Section [D^, can be applied equally well 
to any two-dimensional model. Using this to deduce the existence of as many states as possible 
led — in our case — to the full spectrum. In other cases, we may not be so fortunate, but using 
it as a starting point should make the derivation of the full spectrum a finite (though possibly 
lengthy) task. 



Chapter 4 

Affine Toda Theory 



"It was here that the thaum, hitherto believed to be the smallest possible particle of 
magic, was successfully demonstrated to be made up of 'resorts' (Lit.: 'Thing-ies') or 
reality fragments. Currently research indicates that each reson is itself made up of a 
combination of at least five 'flavours', known as 'up', 'down', 'sideways', 'sex appeal' 

and 'peppermint'." 

— Terry Pratchett, Lords and Ladies 

4.1 Introduction 

"Once upon a time and a very good time it was. " 

— James Joyce 

The sine-Gordon model, which has occupied us for the previous two chapters, is a member of 
the larger family of affine Toda field theories ( ATFTs) , and it is to these that we will now turn 
our attention. These theories are, in general, not as well understood as the sine-Gordon model, 
even in the bulk. 

ATFTs are also integrable, and rely on a Lie algebra structure built into their Lagrangian 
to provide the necessary conserved charges. A tantalising problem with them — and one which 
will provide the basis for this work — is that the underlying structure shows up again in their 
S-matrices, among other places, though it is not at all clear how it arises. The difficulty is 
that the exact S-matrix program, while it provides a good method for obtaining a result, is 
totally disconnected from the original Lagrangian. For the boundary sine-Gordon model, this 
caused problems in relating the parameters in the reflection factors back to the parameters in 
the Lagrangian, while here it hides the path of the Lie algebraic parameters into the S-matrix. 

Obtaining a better understanding of this is still an unsolved problem, but we will find a neat 
method of constructing S-matrix elements through rules based on the Lie algebra, generating a 
number of new identities in the process. 

Before we plunge into the full quantum theory, a preliminary discussion of the classical version 
will serve to introduce much of the structure, as it did for the sine-Gordon model. In addition, 
building an exact quantum S-matrix cannot begin without making some initial assumptions; 
looking at the classical theory will help us make a more educated guess. Classical theories are 
also more intuitively comprehensible, so the more that can be gleaned from them and transferred 
across to the full quantum case, the more tractable it becomes. 
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We shall only attempt a relatively brief introduction to the topic here, sufficient for our 
needs. For a more detailed review and further references see e.g. 



4.2 The Lagrangian 

Affine Toda field theory (ATFT) is a massive integrable 1+1-dimensional theory with a number — 
which we shall call r — of scalar fields 4> a , and with a Lagrangian of the form 

1 m? r 

C = -Wipadpta -^M^2 n j exp(/?oy • <f>) , (4.1) 

1 P 3=0 

where m determines the mass scale (though it does not equate to the mass of any individual 
particle) and (3 is a dimensionless coupling constant. 

The ctj can, in principle, take any values, but it turns out that the resulting theory 
is only integrable if, for j = 1 . . . r, they can be considered as the simple roots of a rank-r 
semi-simple Lie algebra g. 

This is because leaving qo out of the sum gives a conformal theory, known as conformal 
Toda theory or just Toda theory. The possession of conformal (or scale) invariance naturally 
gives such theories an infinite number of symmetries, and hence any conformal field theory must 
be integrable. However, because this means that the theory cannot depend on any fixed length 
scale, all the particles in it must be massless (as, otherwise, the inverse of the mass would 
provide such a scale). Including the extra root to form the "affine" theory can be considered 
as a perturbation which breaks the conformal invariance — and so provides its particles with 
mass — while still retaining an infinite number of symmetries. (Taking cto to be the affine root 
is purely a conventional choice of labelling.) Interest in these theories was initially stimulated 
by this connection to perturbed conformal field theories [pH , and the fact that, through the 
breaking of the conformal symmetry, the particles acquired mass. 

An important feature of such algebras is that they can be conveniently classified [37] in terms 
of a Cartan matrix — C — defined by 

C, = p^-, (4.2) 

where (ai,otj) denotes an inner-product on the roots aj and a~. This matrix encodes the 
relationships between the simple roots, and is particularly simple in that it is composed entirely 
of integer entries. The content of the matrix is often described by a Dynkin [p^j diagram, 
where each simple root is drawn as a "spot" and the spots corresponding to roots cn and ay are 
connected by n lines if CV,- = n. In the case where CV,- ^ Cji, an arrow is drawn on the lines 
pointing from the long root to the short root. 

In Cartan's classification, there are two infinite sets of untwisted "simply-laced" algebras 
(where all roots are of the same length) known as and , with three exceptional cases, 
eg 1 ^, &p and eg . There are also "nonsimply-laced" algebras (where one root has a different 
length to the others) divided into two infinite sets, br and cp , and two exceptional cases, 
and . A listing of their Cartan matrices, together with their Dynkin diagrams, can be found 



in Appendix [B.5| . A good introduction to the topic can be found in |3' 

While it appears on an equal footing with the other simple roots, and it can be drawn as an 
extra spot on the Dynkin diagram to describe its inner products with the other simple roots, «o 
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is not itself simple, in that it can be described as a linear combination of the other simple roots: 

r 

a = - 22 n i a i ■ ( 4 - 3 ) 



8=1 



The m — usually called marks or Kac labels — are integers, chosen to make ^ rnai = 0. (The 
value of cto is prescribed by demanding that this be true with uq = 1.) Two other useful pieces 
of notation are the Coxeter and dual Coxeter numbers, h and h y defined by 

r r 

h = i + ^2 n i and hV = 1 + Yl n ^ ' ( 4 - 4 ) 

i=l i=l 

where the co-marks, n( ', are related to the marks through = nia 2 /2. The Coxeter and dual 
Coxeter numbers will arise frequently in the context of the periodicity of poles of the S-matrix 
or the number of distinct conserved charges. 

Finally, the so-called "incidence matrix" G deserves a mention. This is just the negative of 
the Cartan matrix with all the diagonal terms set to zero, meaning that it simply encodes the 
relationships between the roots. 

Without wishing to delve too far into the development of the theory (further details can be 
found in |T^]), it contains r massive particles, which can be associated with spots on the Dynkin 
diagram of g. Their masses and couplings can be easily extracted from a perturbative expansion 
(in small 0) of the the potential term of the Lagrangian fl4.1| ) : 



2 r 2 r 2 r 

V ^ = 1^2 S ll n 3 exp(/?«i ■(/>) = Jp^^ + ^Yl nirfatj^ 

1 j=0 1 3=0 j=0 

2 a r 

Yl njdjdjofjFtfP + ■■■ (4-5) 

3=0 

This allows us to read off a (mass) 2 matrix 

r 

{M 2 ) ab = m 2 Y^ njtfa) , (4.6) 

3=0 

and a set of three-point couplings 

r 

C abc = 13m 2 , (4.7) 

3=0 

as well as infinitely many higher couplings, at successively higher orders in j3. 

If a basis of fields is chosen so as to make the bare propagator diagonal, (M 2 ) becomes 
diagonal also, allowing the classical masses to be read off as eigenvalues. Finding such a basis, 
and especially computing the three-point couplings in it, is too long a task to be attempted here, 
but it can be done, and closed- form answers obtained 17]. These results, together with other 
relevant Lie algebraic data, can be found in Appendix p.5| . 



4.3 The Quantum Theory 

To find the S-matrix of the quantum theory through the bootstrap approach, we need to begin 
with a suitable guess at one or more of its elements. If, after working through the bootstrap, 
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the result is consistent — each three-point coupling must introduce poles in all three relevant 
S-matrix elements — then the guess could be said to be good. Otherwise, corrections need to be 
made until a consistent result is achieved. 

From the earlier classical results, we might guess that the same couplings transfer across 
to the quantum case, and so predict a minimal pole structure for the S-matrix. One potential 
problem with this approach is that the classical case is the (3 — > limit of the quantum theory 
so, as (5 moves away from zero, the mass ratios, and hence the pole positions, would be expected 
to change due to renormalisation. As luck would have it, moving away from this limit in 
simply-laced cases does not change the position of the poles we have considered so far (one- 
loop calculations showing that the masses renormalise in such a way as to leave their ratios 
unchanged). In an intuitive sense, the bootstrap equations determine the algebraic structure so 
precisely that any continuous change in the parameters (such as the coupling angles) disturbs 
the way the pieces fit together and destroys the solution. Thus the classical mass ratios remain 
even in the full quantum theory. For simplicity, we will go through this case in more detail, and 
just quote the results for the nonsimply-laced cases. 



4.3.1 Simply-laced cases 

The next logical step is to construct a putative S-matrix element with a suitable pole structure. 
A good "building block" for this is provided by 

S inhl(fl + ^) 
{X) sinhi(0-^)- (4 ' 8) 

As mentioned in Chapter [l], this automatically enforces unitarity. It also has only one pole 
(at 9 = and one zero (at 9 = — t^), making it easy to form a suitable product. (In the 
nonsimply-laced cases, the poles are no longer always multiples of t£, so a different block is 
needed.) Crossing symmetry is enforced by demanding a suitable pole structure, and — ATFTs 
being elastic scattering theories — we need not worry about the Yang-Baxter relation, leaving just 
the bootstrap to be satisfied. Building e.g. S\i in this way and working through the bootstrap, 
we do indeed find that it is consistent. 

While this turns out to encode the bound state poles correctly, there is no mention of the 
coupling constant, so it is unlikely to be the complete story. Trying to introduce a dependence 
on the coupling constant leads to the idea that the full S-matrix elements are the elements 
found so far (usually termed "minimal" since they are also the complete S-matrix elements of 
certain perturbed conformal field theories known as minimal models) multiplied by a suitable 
factor. This factor is firstly determined by the fact that the resultant S-matrix must still respect 
unitarity and crossing symmetry, making it natural to also build it out of the (x) blocks. In 
addition, all the necessary bound state poles are already encoded in the minimal S-matrix, so 
the extra factor should not introduce any more physical strip poles, at least for j3 small, though 
it must still respect the bootstrap. 

Finally — and this is the reason why an extra factor is not just an aesthetic invention — at 
(3 = 0, (^7j) shows that all the classical three-point couplings disappear, so the extra factor 
should provide zeros to cancel all the physical-strip poles in the minimal elements, tending to 
them as (3 — ► 0. This means we might be tempted to build the full S-matrix out of blocks of the 
form [x)/{x ± B), where B is a coupling-constant dependent constant. One final complication, 
however, is that the sign of the residue at a pole determines whether it corresponds to a forward 
or cross-channel process. This can be found to be correct for the minimal elements, and must 
be kept so, determining the sign above. 
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Following this through motivates the introduction of an extended block 

{x} = + (4 9) 

W (x-l + B)(x + l-B)' 1 j 

from which the S-matrices of all the simply-laced ATFTs can be built. The S-matrix elements 
are usually written in the form 

h 

Sab(P) = II W ro " b(!e) . (4-10) 
x=l 

where the non-negative integers m a b(x) denote the multiplicity of the block. 

An interesting property of this block is that {x}b = {x}2-B, heralding a duality. Ensuring 
no extra physical poles for real (3 means that < B{(3) < 2, and we have constructed B to 
vanish at j3 = 0, so we might imagine that B — * 2 as (3 — > oo. This would set up a strong- weak 
coupling duality, the theory becoming free in either limit. Determining the precise form of B(/3) 



turns out to be difficult, but it is conjectured to be [41, 17, 42, 44] 



implementing the duality as 



B(^)=2-B(P). (4.12) 







4.3.2 Nonsimply-laced cases 

For these, the S-matrix can be written in a product form [45] as 



h r v h v 

S ab (0) = U H{x,y} m ^y\ (4.13) 

x=l y=l 



where the {x, y} are of the form 



with 



and 



(x- l)(x + 1,2/- 1) 



. , .(2-B)a; By , 



, v sinh (i (0 + iirxS) 

(*)= ... ■ • ( 4 -!6) 
smn (2 (" — nrx)) 

The m a b(x,y)s are again non-negative integers, serving to encode the Lie algebraic structure of 
the model. This time, B((3) is conjectured E3] to be given by 



(4.17) 
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A difference now, however, is that, while there is still a strong-weak coupling duality present, 
it relates the strong coupling regime of one theory to the weak coupling regime of a different 
theory. For example 




(4.18) 



For this reason, the algebras c n and c£„+i are termed a "dual pair" . The simply-laced algebras 
(and a^n) are self dual, and all the other algebras fall naturally into the dual pairs (bit , °2n-i)' 
(cn \ d^+i), (<72 > d^), and (/| , e 6 ). The S-matrices for each member of a dual pair are 
the same except for the interchange of h and h v . In light of this, we will concentrate on the 
untwisted algebras from now on, and drop the superscripts. 



4.4 Lie algebra structure 

S-matrices were first found through this approach for the simply-laced cases EL fij 



and later for nonsimply-laced cases [46, 47]. (The results are summarised in Appendix 
This was all accomplished on a case-by-case basis, and, although there were many hints of the 
underlying Lie algebra in the results, it was not clear how that had been transferred across from 
the Lagrangian. This is frustrating as, apart from the general demands of unitarity, analyticity 
and crossing symmetry, the ATFT S-matrix is principally shaped by the Lie algebra. 

These results were put on a uniform Lie algebraic basis for the simply-laced cases by 



Dorey [48|. He considered the Weyl reflection Wi corresponding to the simple root a%, defined 
by 

2 

Wi(x) = x ir(ai,x)ai. (4.19) 

rv. 

i 

From this, he set w = w\W2 ... w r to be a Coxeter element, with (w) the subgroup of the Weyl 
group generated by w, and defined roots fa by 

fa = w r w r -i . . . w i+1 (oti) . (4.20) 

The other crucial ingredient was a two-colouring of the spots on the Dynkin diagram, where 
each spot was labelled as either "black" or "white" such that no two adjacent spots had the 
same colour. Then, the integers m a i,(x) turned out to be just 

m ab (2p + 1 + u ab ) = (X a ,w~ p fa), (4.21) 

where A a is the fundamental weight corresponding to root a, and u ab = (c(a) — c(6))/2, with 
c(a) = ±1 encoding the colour of the roots. In addition, if we define Tj as the orbit of <j) a under 
(w), then C* jfc ^ iff there exists am G Tj, au\ G Tj, and anA S such that a^ + a^+a^ = 
0. 



These results were initially found by observation. However, Freeman []60| showed how to 
diagonalise the mass matrix in a Lie algebraic way, allowing them to be re-derived more rigor- 
ously. 



Similar results were later found for the nonsimply-laced algebras by Oota [49] through a 
deformation of the Coxeter element. Oota also produced an integral formula for the S-matrix, 
explicitly built from the Cartan matrix, which we shall discuss below. This formula was later 
reproduced by Fring, Korff and Schultz M, while a similar result was conjectured by Frenkel 
and Reshetikhin [50] in the course of a general study of W-algebras. 
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The starting point for our discussion will be the processes shown in figure 4.1. When two 
identical initial particles have a relative rapidity of 29 ^ + t a 0H, these comprise all the possible 
diagrams which can result. The interesting point about them, however, is that it turns out 
that the particles present in the middle of each diagram are always those adjacent to the initial 
particle on the Dynkin diagram. All such particles are present, but no others. 

For the simply-laced cases, only the first three diagrams are relevant, applying to the cases 
where the initial particle has one, two and three adjacent particles respectively. Note that all 
the intermediate particles are parallel and have zero rapidity in the centre of mass reference 
frame. 

For the nonsimply-laced cases, the situation is complicated slightly when the adjacent par- 
ticles are associated with roots which are longer than that of the initial particle. The fourth 
and fifth diagrams describe this for the case where there is only one adjacent particle, and the 
Cartan matrix entry is 2 or 3 respectively. For more than one adjacent particle, the relevant 
vertical line must be replaced with this more complex pattern, as shown in the last diagram. 

The precise makeup of these diagrams is as follows. The first diagram speaks for itself, 
while, in the next (for two adjacent particles), the unspecified particle is always the lightest in 
the theory. The case with three adjacent particles only occurs for d n and e6-8; for which the 



particles are given in table 4.1. 



The next two only occur for c n and g2 , with the particles as shown. The last diagram occurs 
in b n and f±, with the particles given in table 14.21. 



Theory 


a 


b 


c 


d 


e 


f 




n-2 


n — 3 


n — 1 


1 


n 


n 


ee 


4 


3 


2 


1 


6 


5 


e7 


7 


5 


3 


2 


1 


6 


es 


8 


6 


4 


2 


1 


7 



Table 4.1: Diagrams for the cases with three adjacent particles 



Theory 


a 


b 


c 


d 


e 


K 


n — 1 


n-2 


1 


n 


n — 1 


h 


3 


1 


1 


3 


4 



Table 4.2: Diagrams for nonsimply-laced cases with two adjacent particles 

These results come from a case-by-case analysis; it should be possible to derive them from 
the Lie algebraic rule given above, but, for the moment, we have not attempted to do this. 
However, if we take them as axiomatic of how to encode the Lie algebra into the S-matrix, there 
are many consequences. 



4.5 The consequences 

The first important consequence of this result is that, as with the bootstrap relations, another 
particle can be introduced, on a trajectory which either crosses the two incoming particles before 



they interact, or afterwards. This is shown in figure 4.2. Due to factorisation, the amplitudes 



for these two processes should be the same, giving rise to what might be called the "generalized 
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Figure 4.1: Processes which impose Lie algebra structure on the S-matrix 
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bootstrap" 
Sij{6 + 9h + ti0H)Sij(0 — Oh - UOh) = 

\{\{S jl (9 + (2n-l-G il )9 H ). (4.22) 



r Gi, 

-2me(B)Gij 



1 = 1 71=1 

For conciseness, we have denned Oh = ^gir^ and 0// = J^x v , with the integer r v being the 
maximum number of edges connecting any two vertices of the Dynkin diagramF]. The integers 
ti are defined by ti = r v ( ai £*^ ; where the squared length of the short roots is normalised to 2[]. 

This formula was first discovered for simply-laced cases by Ravanini, Tateo and Valleri- 
ani [|5l|] , and was independently derived for the nonsimply- laced algebras in unpublished work 



by the author and P. Dorey []52[ (see also 0) and by Fring, Korff and Schultz ||. 

A subtlety is the exponential factor on the rhs, involving the step function, G, defined by 

if x < 0, 

\ if x = 0, (4.23) 



Q(x) = lim 

€-*0 



1 1 x 
— I arctan — 

2 vr e 



1 if x > 0. 



Due to the periodicity of the exponential, this term has no effect unless 9 = 0, and accounts 
for the fact that, at this point, the additional particle cannot really be said to cross either the 
incoming particles or the intermediate particles. In applications such as the thermodynamic 
Bethe ansatz, it is important that the formula nonetheless continue to make sense at 9 = 0, 
so the extra term is introduced to keep the equation correct. A more detailed discussion and 
derivation can be found in Appendix [B|. 



Another form of this result was used by Oota 4£ ] in his derivation of an integral formula for 



the S-matrix. In Appendix |B.1| , we show that it can be re-stated as 

ml b (x + l)q-^+m q ab (x-l)q^=Y,< c ^)[Gca} q , (4.24) 

c 

where 

m lb( x ) = m ^(x, y)q y , (4.25) 

and the standard notation [n] q = (q n — q~ n )/{q — q~ l ) has been used. 

1 This is 1 for the a, d and e series, 2 for fn and 3 for g2- 
2 Thus ti = 1 for short roots and ti = r v for long roots. 
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This restatement as a recursion relation makes it clear that, with the input of m^ b (0) and 
m g ab (l), all other m a b(x,y) follow. 

These two inputs turn out to be 

m«(0) = m«(l) = q ta [t a ] q 5 ab . (4.26) 

The first of these, m 9 (0) = 0, is trivial for simply-laced cases (as {0}=1). Otherwise, it amounts 
to requiring that poles which are on the physical strip at one value of the coupling stay there for 
all values, which is necessary on physical grounds. (The existence of the three-point couplings 
is not dependent on the coupling, and hence processes which are possible at one value of the 
coupling must be possible for all values.) If {0,y} was present in the S-matrix, for example, it 

would lead to a pole at 9 = (y — 1)9 h — 9h = % ~^p$T^ ~ ^27t • ^ or sufficiently small B, this 
becomes negative. 

The other condition, m q (l) = q ta [t a ] q S a b, implies that the part of the S-matrix coming from 
these blocks is just 

<2,2t a ) (427) 



(2,0X0,2^) • 

The only pole is thus at 29h + 2t a 9n, which is precisely that required for our special processes. 
Thus, this is just the statement that these processes should exist when the two incoming particles 
are identical, but not otherwise. 

In sum, postulating the existence of these special processes, together with basic physical 
requirements, serve to completely fix the minimal S-matrix. These processes, in other words, 
seem to completely encode all the Lie algebraic information contained in the S-matrix. 



The significance of (4.22) can further be seen if we take its logarithmic derivative, and use 



the fact [17] that we can identify the resulting elements with conserved charges. This gives 

E. . 7 / 2 — B Bt{ \ 

[Gu]q {l7Ts) q s = 2 COS ITS I — — + — — I 

1=1 v ' '' 



(4.28) 



The ith. component of the conserved charge with spin s is denoted by q\. The forward-backward 
shifts on the rhs of ( 4.22j ) have been absorbed into the deformation of G, where we have defined 



q(t) (and q(t), to be used later) as 



,'(2-B)t\ , / Bt . 

q(t) = exp ( and q(t) = exp f __ ) . ( 1.20) 



In simply-laced cases, since [n] q = n for n = 0, 1 (as all entries of the incidence matrix are 
in these cases), and we have all U = 1 and h = r v h v , this reduces to the eigenvector equation 

r 

Y^G tl q l s = 2cos{^)qi, (4.30) 
l=i 

a well-known but curious result |j53[| . For nonsimply-laced cases, however, note that the t{ in the 
cos term prevents this from being a proper eigenvalue equation. 



4.5.1 An integral formula 

As well as the product form for the S-matrix elements introduced above, Oota also found an 
integral form, which explicitly builds in the dependence on the Cartan matrix. The proof of this 
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exploits (§31). Since, in consequence, it relies on little other than the Lie algebraic structure 



in the particle couplings, it is perhaps not surprising that Frenkel and Reshetikhin [p(J also 
conjectured a very similar result in their more general study of W-algebras. 

Further details can be found in Appendix (Q), but, for reference, the formula is 

(fOO JL 
4 J —e %ke {smke h ■ smkB H ■ M^qfrk), q{*k)) 

+ S -f}), (4.31) 

The matrix M introduces the dependence on the Cartan matrix and is defined by 

M ij (q,q) = ([K] q7} )7 j L [t j ] 1 , (4.32) 

where K is the "deformed Cartan matrix", given by 

[Kij] m = (qq U + q-'q-^ij ~ [Gij]g . (4.33) 

In the limit q — > 1 and q — ► 1, this recovers the standard Cartan matrix. In some sense this 
can be understood as a quantum deformation, since taking the classical limit (B — > 0) enforces 
q — ► 1. In the simply-laced cases, at least, this reduces the deformed Cartan matrix just to 
the ordinary Cartan matrix with an additional factor proportional to the identity matrix. We 
should also note that our "eigenvector" result can be neatly restated using this, as as 

r 

^2[ K il}q(ms)q(ms)Q l s = ■ (4.34) 
1=1 

To understand what M represents, think that, for the simply-laced cases (where all the ti are 
1), it is just the inverse deformed Cartan matrix. The consequences of the extra factor, which 
modifies it from this, will be seen later. 



The formula given by Frenkel and Reshetikhin [50] is similar to ( 4.3l| ), but without the factor 



of {—l) &ab exp ^/^j ^-e^dabj ■ For real 6, this is 1 except when 6 = 0, in which case it becomes 

-1 for a = b. Including the factor or not thus amounts to selecting the value of S aa (0); with 
the factor, S a b(0) = 1, but without it S a b(0) = (— l) Sab . This second is the value taken by the 
product form S-matrix, and so will be the version we adopt here. 

With a little more work, this can be put into an slightly simpler form. If we first define a 
new matrix <j) as 

<l>ab = -i-ttegS ab (0), (4.35) 
do 

we find 

/oo 
dke ike {4sinfc0 fc • smk9 H ■ M ab {q(Trk), 5(7rfe)) + 5 ab } . (4.36) 
-oo 

Defining its Fourier transform (f) as 

Mk) = — / d60 ab (6)e- ik9 (4.37) 

then leads to 

4> ab (k) = -27r(q(irk) - qinky^iqiirk) - g(vrA;)- 1 )M a6 (( ? (7rA : ), q(nk)) + 2vr5 afe . (4.38) 
(Note that using Frenkel and Reshetikhin's form would have removed the final 5 a b term.) 
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4.5.2 A formula for the conserved charges 

An interesting consequence of the integral formalism — and our reason for introducing it here — 
is that it can be used to find a formula for the conserved charges of the theory, by taking the 
logarithmic derivative of ( [4.31 ) and again identifying it with the conserved charges. Doing this, 
and noting that the resulting integral can be re-expressed as a contour integral over the upper 
half-plane, the problem is reduced to finding the poles of the expression. The only poles are in 
the matrix M, so, before we can continue, we must find a formula for this. The easiest route to 
the information we need is to compute <fi a b(k) for the product form and compare with the above. 

The first step — calculating (p a b — is straightforward, and yields 



h r v h x 



b ab = - 9 m ab{ x iV) 



x=l y=l 



£ {tanh^ + s^S/)))}- 1 

s(x,y)eSi 



where 



{tanh(±(# + *'(*, y)))}- 1 

s'{x,y)eS 2 



, (4.39) 



S x = {(x-l)e h + (y-l)8 H ,(x + l)6 h + (y + l)9 H ,(l-x)6 h -(y + l)8 H , 

(1 - x)6 h -(y + l)d H , ~(x+ l)9 h + (1 - y)9 H } , (4.40) 

S 2 = {(l-x)e h + (l-y)6 H ,-(x + l)0 h -(y + l)e H ,(x-l)e h + (y + l)e H , 

(x - l)9 h + (y + 1)9 H , (x + l)6 h + (y- 1)9 H } . (4.41) 



The Fourier transform of these terms is given in Appendix as 



tanh I — ha 



e -ike M = _ 4nie -ik(2a+inn) 



cosh(7rfe) 
sinh(27rA:) ' 



(4.42) 



where care must be taken to choose n such that there are no poles between the real axis and 
the line 2a + inn. Working this through finally gives 



'Mb 



h r v h v 

(k) = 27r5 ab - 27r ^2 m <*b( x > J/XsC 71 "* 5 ) ~ 9(P^)~ 1 )(9( 7rfc ) - Qinky 1 ) 
x—\ y—1 



q x q y — q x q y 



l —± 1 1 (4 43) 

1 _ q 2hq2rVhV ' V*- 1 * ) 



and hence 



h r v /i v 

M a b(q,q) = ^2^2 m «*( aJ ' y)^7 
x=l y=l 



q x q y — q x q y 

q2hq2r v h v 



(4.44) 



This shows that the only poles present are at k = im, m being any integer, so the result 
is that we can re-express the integral in the form of a Fourier expansion, and thus read off a 

(s) 

relation between ip ^ and M as 



^Pab = 2sin7rs • smhsOh ■ smhsOu • M(q(iirs), q(iirs)). 



(4.45) 
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Of course, to find an expression in q^q h s: we need to include a scaling factor. Noting that 
X)I=i Is Is = dab, where 8{ is the ith component of a rank-r algebra, we could use q®q h s = 

Combining this with the expression for M, we get 

h r v h v 

ip^ = 2 sinh s9 h • sinh s6 H • ^ ^ ?7iafe(a:, y) sin 

X=l J/=l 

Prom this, it is straightforward to see that the matrix ip y V is non-zero for generic B by simple 
case-by-case analysis. (This is different from this minimal case where, as noted by Klassen and 
Melzer |53|| , we can get a zero matrix for s = § in simply-laced cases, even if that exponent 
is present.) Had there been cases where was zero for some s, then taking the logarithmic 
derivative of an S-matrix identity would sometimes have resulted in a trivial conserved charge 
identity. As it is, however, we can always derive a non-trivial conserved charge identity from an 
S-matrix identity and vice versa. 



(2 



B)x + By 



r v /i v 



(4.46) 



4-5 The consequences 



77 



4.5.3 Multi-linear Identities 

"Life must be understood backwards; but ... it must be lived forwards. " 

— Soren Kierkegaard 

The RTV result and its generalisation allow us to perform a simple trick and generate a large 
number of S-matrix identities. Interchanging i and j in ( 4.22p does not change the lhs if ij = tj 



- the two roots are the same length - due to the symmetry of the S-matrix, so we can equate 
the rhs before and after interchanging to get 

r Gu r G jV 

n n + (2n - 1 - g U )0h) = n n s *< e + »' - 1 - (4.47) 

!=ln=l l'=\n'=l 

(Note that the presence or absence of an exponential factor does not affect this, as ti = tj ensures 
Gij = Gji.) If i and j are such that the corresponding rows of the incidence matrix consist of 
entries no greater than 1, this reduces to 

r r 

Y[S tt (pfv = Y[S jV (e) G «, (4.48) 
1=1 i'=i 

and we can obtain identities for products of S-matrix elements, all evaluated at the same rapidity. 
The existence of such identities was first discovered by Khastgir Q, though without such a 
systematic method for describing them. In addition, we also have identities in which not all 
rapidities are equal. 

To generalise the connection between S-matrix product identities and conserved charge sum 
rules to this case, we can take logarithmic derivatives to find that if 

[] S ob {9 + ifl b )= SM6 + iftv), (4-49) 

a,be{i,j} a',b'e{i',j'} 

for some sets {i,j} and {i',f} then 

E e-<&'tfg*= £ e-^V^'gf. (4.50) 

a,b£{i,j} a',b'<={i',j'} 



Applying this to ( 4.47]) gives 

r r 

^2[Gii]q( i7TS )q l s qi = y^,[Gjl']q(iTr S )Q l s<ll, (4-51) 



1=1 l'=l 



where it should be noted that the sums over n and n' in (4.47) have been absorbed by the 
introduction of the [G a &]g(7rs) notation. 

To give a simple example of this result, in the br algebra we have, for 1 < i < r — 1 

S(r-l)(i-l){8)S( r -l)(i+l){0) = Si( r _2)(0)Sir(9 + 0H)Si r (9 - 6 H ), (4-52) 

and 

+ flT 1 ^ 1 = lWs~ 2 + \ MB • ills: (4-53) 

(2) 

with (through the duality transformation B —> 2 — B) corresponding identities for flt2r_l" 
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It is still an open question as to whether we have found all such identities, or merely a 
subset, but there is good reason to believe that these represent all that can be found. From the 



multi- linear identities ( 4.47 ) come all possible identities involving shifts only depending on 6h 
while bringing the full machinery of the generalised bootstrap into play ultimately allows the 
proof or disproof of any identity. 

In the first situation, case-by-case analysis shows that the first row of all the S-matrices 
consists of linearly independent elements, as each has at least one pole which is not found in 
any of the others. If our identities provide a way to re-write all the other S-matrix elements in 
terms of this set, it can be used as a basis. Any other identity can then be proved or disproved 
by expanding it in the basis, and comparing terms. 

In general, this idea works very well. The only difficulty arises for d n due to the pair of 
degenerate particles. For n odd, the elements S n r n -i)(0) and S aa (6) (for a = n — l,n) cannot 
be separated, and the best that can be done is to say 

n-2 

Saa(e)S n(n _ 1) (6) = J] S lp (6). (4.54) 

p=l step 2 



For n even, this separates into 



Sn(n-1)(0) = J] S ^ ( 4 ' 55 ) 

p=3 step 4 
n—A+x 

Saa(9) = [J Sl P^> ( 456 ) 
p=l step 4 

(where x = 1 for n divisible by 4, and x = 3 otherwise), but this cannot be done for n odd. 
However, in this case, e.g. S n ( n ^u(9) becomes linearly independent of the first-row elements, 
and so can be added to the basis, allowing the argument to still be used. 

For the more general situation, the generalised bootstrap (in common with the usual boot- 
strap) allows the entire S-matrix to be built from an initial knowledge of one element, usually 
Sx\. This means that any other element can be written in terms of Su (with various forward- 
backward shifts) by repeated use of the bootstrap. Inserting this into any identity to be proved 
then reduces it to a product of elements Su with a variety of rapidities. If these are linearly in- 
dependent of each other (as seems reasonable) then simply comparing terms would be sufficient 
to prove or disprove the identity. 

Neither of these arguments is as rigorous as we would like, but they do hold out the reasonable 
possibility that the claim might be true. This would reinforce the idea that all the structure in 
the S-matrix is due to the underlying Lie algebra. 



4.6 Summary 



an ill-favoured thing, sir, but mine own ... 

— William Shakespeare 



The aim of this chapter was to find a concise way of encoding the Lie algebraic information into 
the S-matrix of all ATFTs. This was achieved by looking at the processes responsible for poles 
at 29h + 2U9h whenever the incoming particles were identical. These could be explained by 
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figure where, crucially, the intermediate particles consisted of those adjacent to the initial 
particles on the Dynkin diagram on the algebra. 

Sending in a third particle either before or after the interaction, and using the principle of 
factorisation to equate the results led to the "generalised bootstrap" 



Sij(9 + e h + Midship -e h - ue H ) 



-2ra-e(0)Gi 



r Gu 



Iin^+(2n-l-W- (4.57) 



1=1 n=l 



This, together with demanding the existence of these processes (and their associated poles) 
completely fixes the minimal S-matrix. The remaining question, however, is how the processes 
arise from the initial Lagrangian formulation. 

Taking the logarithmic derivative of ( |4.57| ) then leads to an equation for the conserved charges 
of the theory, namely 



^2[Gil]q(i7Ts)Q l s = 2cos 



ITS 



2 - B Bti 
+ 



2h 2r v /i v 



(4.58) 



which reduces to a simple eigenvector equation in simply-laced cases. These charges can also be 
written as 



h r v h v , 
q a s q b s oc 2 sinh sO h ■ sinh s6 H ■ ^ ^ rn ab (x, y) sin ( — 

x=l y=l 



(2-B)x | By 



r y h y 



(4.59) 



Since the S-matrix is symmetric, the lhs of ( [4.57 ) is unchanged by interchanging i and j, 
whereas the rhs is not, leading to the identities 

n n s 3 i(o + (2n - 1 - guWh) = n n + ^ - 1 - g ^h) , (4.60) 

(=ln=l i'=ln'=l 

which probably describe all identities with shifts only involving just as the generalised 
bootstrap contains enough information to prove or disprove all possible identities. 



Chapter 5 

Conclusions 

'"Good morning, ' said Deep Thought at last. 

'Er ... Good morning, O Deep Thought,' said Loonquawl nervously, 'do you have ... er, 

that is 

'An answer for you?' interrupted Deep Thought majestically. 'Yes. I have.' 
'To Everything? To the great Question of Life, the Universe and Everything?' 

'Yes. ' 

'Though I don't think, ' added Deep Thought, 'that you're going to like it. ' 

'Doesn't matter!' said Phouchg. 'We must know it! Now!' 

'Alright, ' said the computer and settled into silence again. The two men fidgeted. The 

tension was unbearable. 

'You're really not going to like it, ' observed Deep Thought. 

'Tell us!' 

'Alright, ' said Deep Thought. 'The Answer to the Great Question 

'Yes ...!' 

'Of Life, the Universe and Everything said Deep Thought. 

'Yes ...!' 

'Is ...' said Deep Thought, and paused. 
'Yes ...!!!...?' 

'Forty-two,' said Deep Thought, with infinite majesty and calm." 

— Douglas Adams, The Hitch Hiker's Guide to the Galaxy 

5.1 Introduction 

The aim of this study was to investigate the fundamental objects of ATFTs: the S-matrices of 
the theory in the bulk, and the reflection factors of the theory with a boundary. 

For the bulk theory, the form of the S-matrices and the particle structure were already well- 
known; the intriguing question was how the Lie algebraic structure built into the Lagrangian 
manifested itself in the S-matrix. For the boundary theory, on the other hand, even for the sim- 
plest possible ATFT — sine-Gordon — the reflection factors for all except the ground and lowest 
excited states of the theory were unknown, as was the boundary bound state structure. 
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The focus for both pieces of work could therefore be said to be their bootstrap structure: 
tying it in to the underlying Lie algebraic structure in the bulk; and finding a rigorous way to 
identify the bound states hidden in the boundary reflection factors. 



5.2 Bulk ATFTs 

This work was based on the observation that, for any ATFT, two identical particles (say i) 
colliding at a relative rapidity of 29^ + Itfin results in the production of all the particles which 
are adjacent to it on the Dynkin diagram, and only those. Taking these processes as a starting 
point, a "generalised bootstrap" was constructed, which explicitly related the structure of the S- 
matrix elements to the Cartan matrix. By using these equations, together with the requirement 
that no more couplings than necessary be introduced, it was found that the complete minimal 
S-matrix could be derived. 

The weak link is that these processes have been introduced as axiomatic, rather than via 
a derivation from Dorey's Lie algebraic coupling rule. An important open problem is whether 
our simple tree-level argument will stand up to a perturbative verification to higher loops in 
the Feynman diagrams. We can, however, gain some measure of confidence from the fact that 
there is substantial evidence for the validity of the S-matrix formulae, which are successfully 
reproduced. 

With this in place, it should then be possible to tie in all the other results which have been 
found by observation on a more rigorous basis. 



5.3 Boundary sine- Gordon 

The task here was more basic: a determination of the bound-state structure and reflection factors 
for all integrable boundary conditions. This was achieved, principally with the help of two rather 
general lemmas which showed that poles at sufficiently small rapidities could not correspond to 
anything other than a bound state without violating momentum conservation. By taking the 
spectrum to consist of just the states which were required to satisfy the lemmas, we could then 
show that all the other poles had an explanation through the Coleman-Thun mechanism. 

Since the lemmas are quite general, they apply to all theories with a boundary, integrable 
or not. An interesting open question is whether, as here, they are strong enough to completely 
determine the spectrum, or merely provide a starting point. 

The natural way to continue the work would be to generalise it to other ATFTs. It has 
been found [17| that, at the so-called "reflectionless points" (which occur at integer A) the full 



S-matrices for sine-Gordon become the minimal matrices for the d\ theory. The soliton and 
anti-soliton correspond to the two mass-degenerate particles, while the breathers correspond 
to all the others. This might make the extension of the results found here to d\ relatively 
straightforward. However, results for a\ have already been found [59|, and indicate that the 



coupling plays a bigger role in the boundary spectrum than in the bulk. Thus, while all the 
coupling information is contained in the minimal S-matrix for the bulk, the story is probably 
not so simple with a boundary. However, it might still provide a good starting point. If this 
could be achieved, only the exceptional cases would remain to complete the ADE series. 

At this point, the position for the boundary theories would be analogous to that for the 
bulk, in that the next logical step would be to put everything on a manifestly Lie algebraic 
footing. While a unified discussion of all boundary ATFTs is perhaps still some way off, it 
should nonetheless be an attainable goal. The theories could then be said to be under complete 
control, at least from this point of view. 



Appendix A 

Boundary sine-Gordon Details 



"This is a one line proof... if we start sufficiently far to the left." 

— peter@cbmvax.cbm.commodore.com 



A.l Infinite products of gamma functions 

The products which arise in the course of this work are of the form 

T(kl + a-xu)T(kl + b-xu 



P(«) = n 



/(« -»• -u) 



^ \ Y{kl + c - xu)Y{kl + d-xu) 
Rather than examine this product directly, we take logs and use the standard formula 



(A.l) 



In T(z) = z ln(z) - z - - \n(z) + ln( y/2) + — + 0(z~ 3 ) (A.2) 

Assuming that the sum over I and the expansion in z can be exchanged, potential divergences 
arise from terms of the form F w ^ n a / an d n < 2. To begin with, we will consider the 

terms arising from the block of four terms explicitly shown. 

Firstly, there is a contribution of Yl^i a + b — c — d from the z terms, which can be set to 
zero by demanding a + b = c + d. For the l/12z terms, the overall contribution from the four 
terms is 

El / a — c b — d \ 

, x 12 V (kl + a - xu)(kl + c - xu) + (kl + b - xu)(kl + d-xu)J 

which can be seen, for a + b = c + d, to be of the form l/l 2 and hence convergent. 

A similar argument applies to the —^ln(z) terms, showing they also provide a convergent 
contribution. This breaks down when considering the zln(z) terms, however, and their contri- 
bution formally reduces to 

E(=£W.)). <A.4) 

1=1 

which is divergent unless a = c or b = c, both of which are trivial cases. However, repeating this 
argument on the other block (with u — > —it) can be seen to yield the same result, allowing the 
two divergent terms to cancel, and leaving a product which is convergent overall. 
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For comparison with other results, it is useful to write P(u) in other ways. Firstly, it can be 
written in terms of Barnes' diperiodic sine functions using the expansion as given in pip: 



S2(x\ux,u)2) = exp 



(wi + uj 2 - 2x) (7 + log(2vr) + 2 log ' ^ 
2^ 



p J W 1 +U12— x 

r(-) 2 



0)1 



L02 



, (A.5) 



where 7 denotes the Euler constant. For blocks of the form we are interested in, this simplifies 
to 



5*2 (x 3 \U>1 , W 2 ) 5 2 (x 4 1 UJi , L0 2 ) 



n 

n=l 



p / raa;2 1 ^1—^2 



2wi 



2£Jl J 



p / nti>2 1 1^1—^2 



x 3 ) p I n^2 



p / nu!2 1 tJi — 1£2 

2uji y \ ui 2wi 2oji / 



/ (*^m * x rri) 



(A-6) 



(where x'^ 



Xr, 



lo\ — U2) provided xi + X2 = X3 + X4. Comparing with (|A.1|) we have 



P(u) 



52(wi(l — a + xu)\u>i,uJik)S2(iOi(l — b + xit) wife) 



52(^1(1 — c + uik)S2(u>i(l — d + xu)\{jj\,u)\k) ' 

where w\ is arbitrary. In section [3.2.2 we took oj\ = x -1 for simplicity. The identity 



52 (wi + 0J2 — x\u\,U2) 



1 



S 2 (x\uJl,UJ2) 

was also used. 

These products can also be written in an integral form, through 



(A.7) 



(A. 



iogr(c) 



dx 
— ( 

x 



C-1 + 



-(c-i)* _ 1 

1 - e~ x 



, Re C > 0. 



(A.9) 



Since, for the expressions we consider, not all the T-functions have arguments with positive real 
part, it is not possible to give a general formula for P solely in these terms. Instead, we give 
expressions for the reflection factors. To simplify matters, define 



I\u) 



I 2 (u) 



2A f + °° , , {2\ux\ 
— / dx cosh 

J-00 V t / 



2A f +QO , , f2Xux\ 
— / dx cosh 

tt I 1 / 



sinh ( A — — 



2 sinh x cosh Ax 
sinh ( ^1 — 277.* - 



2 x 



2 A f + °° , , /2A«xA 
dx cosh ( I 



7T 



sinhx 

2 cosh x sinh (A + l + 2n — — )x 



(A.10) 



2A f + °° , , /2Aux\ 
/ ax cosh 

71" J-00 V 7T / 



2 sinh x cosh Ax 
2 cosh x sinh ( — + 2n — A — l)x 



2 sinh x cosh Ax 



A. 2 Relation of M and (po to r\ and i? 
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(where i^(tt) and I^(u) are related to each other through £ — > 7r(A + 1) — £). The constant n* 
is the number of z^-type poles in the physical strip, which we recall can be written as 

7T 2 



The reflection factors can then be written as 



d 


log 




du 


Ro(u) 


d 


log 




du 


Ro(u) 



i odd j even 

/1( U ) - (1 - c) /2 (n) + £ J* („) + £ £ (f|) > 



odd 



j even 



(A.ll) 



(A.12) 



for topological charge c and x = (ni,n,2, ■ ■ ■ ,n2k+c)- These factors were given in [29] for the 
first two levels of excited states (the extent of the spectrum they found); the above is simply a 
generalisation of this to the whole spectrum. 



A. 2 Relation of M and (p to r] and d 



For the action defined as 



/U poo j roo 

dx / dt —{d^f + 2^cos(2/?v?) + 2fi B / dtcos /%(0,t) - <? ) , (A.13) 
-00 J— 00 J— 00 

Zamolodchikov [54] has claimed that 



cosh 2 (/? 2 (i?±M7)) 



Mb 



sin(-7r/3 )e 



2n ±2i/3 W 



(A.14) 



where this should be read as two equations, one with the positive signs, and one with the negative. 
To match our conventions, we need to re-scale this according to ip — > \/2vr<^, po ~^ V^^o and 
j3 -» /3/2 v / 2vf. Then we need to identify \i = mg/2/3 2 and \i B = M/2. This means that 
Zamolodchikov's formula becomes 

2 



P 

cosh 2 ( — (1? ± ir/) 

87T 



2 V mo 



8111 



(A.15) 



This result agrees with earlier results for special cases |57|] . To get a better idea how the two 
sets of parameters are related, it is useful to deconstruct ( |A.14| ) into equations for M, <po, rj and 
•& individually, giving each in terms of the other set of parameters. For the first two, we get 



¥0 



M 



cosh (§=('& + iff) 



cosh (£ ~ ir l)) 
2ml cosh (£ ^ + ir] >) cosh 



ff 1 



sm v 



(A.16) 



where the choices of sign must match, and are determined by requiring M to be real. (If we 
take 77 and t? real, this means we must take the positive sign.) For n and 1?, the task is made 
much easier by introducing the change of variables given by 



cos 



8tt 



cosh 



cos 



0hL 

8tt 



+ cosh 




(A.17) 
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where, in the limit j3 2 — > 8tt, £g —* £ and kp ^ k . In terms of these, we find 

cos 2 £3 = ~ (1 ± cosO^o)) (A.18) 

where the choice of sign is as above. The parameters r\ and i? are then determined by 

cosh O^p) > cos (l^) = ± \/*$ ~ 2 ^ COS 2 ^ + 1 ' ^ A ' 20 ^ 

where, in principle, either rj or $ can correspond to either choice of sign, and the sign here is 
unconnected to the earlier choice. 

For sine-Gordon theory, (3 is taken to be a strictly real parameter. The boundary parameters, 
M and (fo, must also be real to keep the boundary potential realQ. The important point to note 
is that this means that the rhs of ( A.20| ) is purely real, forcing tj and $ to either be real or 



purely imaginary. In addition, the choice with the negative sign has modulus less than or equal 
to 1. This means that there is always a choice of rj and where both are real. The symmetry 
between rj and ii9 makes the choice where both are purely imaginary equivalent. The remaining 
two choices — where one is real, the other imaginary — make the lhs's of ( A.14| ) real, while the 



rhs's are complex conjugates of each other, and so are untenable. Thus, we can take rj and i9 
real without loss of generality. 



In the Dirichlet limit, i.e. M — > oo, we have — > ±oo also, reducing ( A.20| ) to 



This gives $ — ► ±oo and 



f3 2 T] 

cok]i(^) = 2A$ - cos 2£ ,• . (A.21) 



2 v 

^ = nvr ± (3tp , (A.22) 



for any neZ. All of these choices for r\ correspond to the same physical reflection factors, so 



we can take rj = -^tp, recovering the result conjectured in | 13[| . 



A. 2.1 Comparison with other results 

The sine-Gordon theory can be considered as the continuation to imaginary coupling of the sinh- 
Gordon model. For this model, an independent proposal for the relation between the parameters 



in the lagrangian and the reflection factors was made by Corrigan [55] and was futher studied 



in [ 56 1 . The field equation used there was 



with boundary condition 



Ofy - d 2 x 4> + ^f- sinh(V2/3<^) = , (A.23) 

ft 



V2m ( -4#o,t) r „-^0(o,t) 



Mo = ^ ( e e~^> - £l e^^ ) , (A.24) 



1 Allowing them to be complex (M = M r + iMi and 950 = (for + ifoi) leads to the demand M r sinh 
Mcosh = if the potential is to be kept real. The only solution to this is Mi = ipoi = 0. 



2 
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where the boundary parameters are eo and e±. The parameters in the reflection factors were 
then found to be 

Bn iB-Q 

-J- = (a + ai)(l - B 12) and = (a - oi)(l - B 12) , (A.25) 

7T 7T 

where B was related to the coupling constant by B = 2/3 2 /(4-7r + /3 2 ), and oo, ai was given by 

eo = cos7rao ei = coS7rai. (A. 26) 

Their conventions differ from ours in the bulk by the transformations (f> — > 4>/2, [3 — > V2(3, 
and m — > mo / y/2/3- Applying these to the boundary condition gives 

Mo = ™o//? 3/2 (e<fl-i*W> - 6 ie l^ -*)) . (A.27) 
To finally turn this into a suitable form for comparison, we need the trigonometric identity 

ae b + ce~ b = ^fac cosh(b + d) if cosh d = ^-^£ . (A.28) 

2y/ac 

After a little algebra, and continuing (3 to i/3, their results then become 

?2„ 



cos(^j = e e 1 -^l-el-e\ + (A.29) 
cosh(^) = eoei + ^/l-eg-ef + ege?, (A.30) 



for the boundary condition 



2m 0v / e^eT { 0<p(O,t) . » . . 

^1° = — ^2 — cos ( 2 za ) ' ( A - 31 ) 



where the value of a is 



cosha = J— =. (A.32) 
2 v /£ oei 

To match this boundary condition to ours, we need to identify 

6061 =^^7 J = 4sin(/378) ^ • (A - 33) 

More algebra then shows that this gives 



cos 



47T 

2 



where 



= Kp — - 2k 2 C os 2£/3 + 4jz 2 p + 1 (A.34) 
coshf^-^j = + y^/cj - 2^2 cos 2^ + 4k 2 + 1 , (A. 35) 



This is very similar to ( A.20p , but it is not quite the same. This does not necessarily mean that 



either is wrong; the differences could simply be down to e.g. implicit choices of renormalisation 
scheme in the derivation of the respective formulae. The resolution of this question is still open. 
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A. 3 On-shell diagrams 

In this appendix we collect together some of the on-shell diagrams used in the main body of the 
thesis. All boundaries are initially in the state \ni,n.2,--- ,ri2k), where k can be any integer, 
and we have suppressed the topological charge index (which is zero). Analogous processes for 
charge 1 states can be found by applying the transformation £ — * tt(X + 1) — £ to all rapidities 
shown. 

In addition, where the boundary is shown decaying through emission of a breather, only 
the process where this removes the last two indices is given. Similar processes always exist to 
remove any other adjacent pair of indices, or to simply modify an index; see section |7rJ for 
the appropriate breather boundary vertices. 





Figure A. 3: Incoming soliton, soliton Figure A. 4: As [A.3| with incoming 
boundary decay soliton crossed 
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rv 



|ni, . . . ,n 2k ,m) 



2 n.-^ n ~ 2\ 




Wn 2k +n + 7T 

m,... ,n 2 fc_i) 



Figure A. 5: Incoming breather, soliton Figure A. 6: As [A.q with outgoing 
bound state soliton crossed 




m,... ,n 2 fe-i) 

~ w n 2k +n 



n2k 




> I 2 k 



ni, . . . ,n 2 k-i) 

~ w n 2 k+n 



Figure A. 7: Incoming breather, soliton Figure A. 8: As A. 7 with incoming 
boundary decay breather crossed 



7T 




W0-V2„ 2k _ 1 +n 
2 

I = n 2 k + «2fc-l 

_Z// n 2 fc_i+n 

ni, . . . ,n 2k -2) 

2k-l ~ Wn 2k 



2-IT+V2n 2k _ 1 +n-WQ 



7T 



■ . n 



n 2A 




Figure A. 9: As |A.5| , outer legs replaced Figure A. 10: As |A.8| , outer legs 
by A.l replaced by |A.2 
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a = 2lT 
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'2k-l~ Wn 2k 
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2 

ni, . . . ,ra 2 fc-i) 
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Figure A. 11: As A. 6, outer legs 



replaced by A. 2 



Figure A. 12: As A. 5, outer legs 
replaced by [A. 3 




Figure A. 13: As 



A.2 



, outer legs replaced by all-breather version 



Appendix B 

Miscellaneous Proofs 



"The trouble with facts is that there are so many of them." 

— Samuel McChord Crothers 



"Basic research is what I am doing when I don't know what I am doing." 

— Wernher von Braun 

In this appendix, we present various proofs which are subsidiary to the main text, but serve 
either to fill out the bare bones of it, or provide cross-checks on the results presented. 

B.l Oota's starting point 

Oota, in deriving his integral formula for the S-matrix, began by defining 

m lb( x ) = ^2 mafe ( x ' y)v V > ( B - 1 ) 

as well as the matrices 

(D q ) ab = q ta 5 ab , (T q ) ab = [t a ] q 5 ah , (I q ) ab = [G ab ] q . (B.2) 
He then stated (after a case- by-case analysis) that the matrices m q {x) satisfied 

m q (0)=0, mi(l) = D q T q , (B.3) 

as well as the recursion relation 

D- x m q Jx + 1) + D q m q ab (x - 1) = I q m q (x) . (B.4) 



As we shall see, the recursion relation follows from the generalised bootstrap ( [4.22 ). Exam- 
ining first the recursion relation, note that it can be re-written as 



m l b 



(x + 1)<T** + m q ab {x - l)q ta = ^2[G ac ] q ml b (x) . (B.5) 
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Turning now to (4.22), we can use the product-form notation to re-write the rhs as 



h r v hv ( { x ,y} mb ^ x 'y\G ac = 1 
rhs = H H H I ({x,y-l}{x,y + l}) m »^>y\G ac = 2 . (B.6) 

c,G a ^0x=i y =i ( ({x,y-2}{x,y}{x,y + 2}) m >>^y\G ac = 3 

This comes about because the forward-backward shift on the rhs has the effect of shifting y for- 
ward and backward by 1 or 2, though it should be noted that this is not quite as straightforward 
as it seems and, for example, a forward shift on its own does not have the effect of producing 
any neat shift in y. To put this another way, 

h r-h- ( {x,y} m ^y\G ca = l 
rhs = IT nil) {x,y} m »^ x >y-V +mb c ( - x >y +1 \G ca = 2 . (B.7) 

c,G ca ^0x=l y=l [ ^ x ^ y jm bc (x,y~2)+m bc (x,y)+m bc (x,y+2) = 3 



Looking at the lhs of ( [4.22 ) and writing it the same way, we find 



h r v ft. v 

lhs = Yl [J {x,y} m ^ x -^ y - ta)+m ^ x+ ^ y+ta) . (B.8) 

x=l y=l 

Comparing block multiplicities, this reduces to 

m ab (x - 1, y - t a ) + m ab (x + 1, y + t a ) = 

{m bc (x,y),G ca = 1 
m b c(x,y - 1) + m bc (x,y + l),G ca = 2 . (B.9) 

m bc (x, y - 2) + m bc (x, y) + m bc (x, y + 2),G ca = 3 

Multiplying through by q y , this can be rearranged to 

m ab (x - 1, y - t a )q y - ta q ta + m ab (x + l,y + t a )q y+ta q- ta = 

{m bc (x,y)q y ,G ca = 1 
m bc (x,y-l)qy- 1 q + m bc (x,y + l)qy+ 1 q~ 1 ,G ca = 2 . (B.10) 

m bc {x, y - 2)<f " V + m bc (x, y)q y + m bc (x, y + 2)<f +V 2 , G ca = 3 

Summing both sides over all integer y, we can then re- write this in terms of the matrices m q (x) 
as 

{m q hc {x),G cb = 1 
m g bc (x)(q + q- 1 ),G cb = 2 , (B.ll) 

ml(x)(q 2 + l + q- 2 ),G cb = 3 

which, noting that [n] q = q n ~ x + q n ~ 3 + . . . + g^" -1 ), is just (|B.5|). 



B.2 The generalised bootstrap at 9 = 

The generalised bootstrap is naively 

r G a i 

s ab (o + e h + t a 9 H )s ab (9 -e h - t a e H ) = J] JJ s bl (e + (2n - 1 - G al )e H ) . (B.12) 

1=1 n =l 



B.3 Check that generalised bootstrap follows from integral formula 
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A subtlety arises when we consider the case 9 = 0, since we can either consider the lhs as 

lim lim S ab (9 + H)S ab (9 - H) (B.13) 

e-+oH->8 k +t a e H 

or as 

lim lim S ab (9 + H)S ab {9 - H) . (B.14) 

H^e h +t a e H e^o 

For 8 ^ 0, this distinction makes no difference, since we are not near a pole of S, but at 9 = 
we are potentially considering a pole, and hence way the limit is taken is important. By leaving 
9 arbitrary and fixing H from the beginning, we have been implicitly using the first form, but 
is perhaps more sensible to consider H as a shift from S(9) — as would be the case in the Bethe 
ansatz approach — in which case the second form would be more appropriate. 

To see the difference between these two forms, we can consider one basic block of S, (x) for 
some x. Shifted forward and back by H, this becomes 



sinh \ (9 + H + i ttx) sinh \ {9 - H + in x) 



sinh t;(9 + H — iirx) sinh ^(9 — H — iirx) 



(B.15) 



For the sake of argument, we shall take H positive. It is clear that as long as H ^ iirx, all the 
arguments of the sinh functions are non-zero (noting that we need not worry about periodicity 
as x < 1) and thus the result is the same in either limit. However, at H — > iirx, a discrepancy 
arises as, if we take this limit first, we find 

sinh±(0 + 2il) sinh±(0) 

P V- 2 - (B.16) 

sinh±(#-2F) sinh±(#) 

which reduces to -1 if we take the 9 limit as well. Taking the limits in the other order, however, 
we get 

sinh I (H + iirx) sinh i(— H + iirx) 



sinh \{H — inx) sinh ^(— H — inx) 



(B.17) 



which reduces to 1 even before taking the H limit. 

For 9 = 0, then, if we want to take the 9 — > limit first, we must modify (|B.12) by a factor 



of -1 for every "problem" S-matrix block, i.e. for every block of the form (x,y) = (l,t a ) in Sab. 
Going to the larger block, {x,y}, this turns out to mean a factor for every block {2,t a ± 1}. 

The easiest way to go from here is to appeal to ( |B.5| ), with x = 1, which gives 

q~ ta ml b (2)=q t nt b ] q [G ba ] a (B.18) 

or 

m q ab {2) = ^(gfc-l + q^ 3 + ■■■+ q l -^)[G ha ] q . (B.19) 

We are now looking for terms in g* a±1 in this expansion. If G ba = 1, then the lowest term is q ta+1 , 
meaning we need one minus sign. If G ba = 2, we introduce a factor of (q + q^ 1 ), leaving us with 
terms like q ra and q Ta+2 , but none of the right form. If G ba = 3, it is simplest to note that we 
must therefore be looking at G2, and that r a = 3, r b = 1, showing m q ab {2) = q 8 '+2q 6 '+3q 4 +2q 2 '' + 1, 
with us searching for powers of q 2 or q°. Thus, G ab = 3 leaves us needing to introduce a net 
minus sign as well. 

To summarise, we need to introduce a minus sign to one side of (|BTl2|) for G ab oddQ in the 



case 9 = 0; the term used in ( f4.22| ) is perhaps as good a way as any, and turns out to be useful 
in further calculations. 



1 Note that G a b is necessarily odd if Gba is, though the two need not be equal. 
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B.3 Check that generalised bootstrap follows from integral for- 
mula 

The most straightforward method (and the one we shall use) is to propose an identity of the 
form 

r Gu 

SijiO + h + U0H)Sij(6 -6 h - UB H ) = e y II S d e + (2n - 1 - Gu)6 H ), (B.20) 

1=1 n=l 

and aim to find y by substituting in the integral formula for the S-matrix. 

Since equation (|4.22| ) applies to the case where the limit is taken first, we need a prescription 
for taking the other limit. It turned out to be simplest to replace Oh and Oh in ( B,20| ) by Oh + ie 



and Oh + ic : and take the limit e — > last. Substituting in ( 4.3l| ) and simplifying, we find 



e y = lim exp ^ / ^e ik0 {q(irk) - q{-nk)} {q(irk)^ - q(-irk)^ } • 



J = l 



[Kajq'^q'^Mijiqiirk), q{irk))) (B.21) 



where q'(t) = q(t)e-* and q'(t) = q(t)e^ . Looking back to ( [4.44 ), we can see that when the 



integrand in ( |B.21| ) is expanded out, all the terms are of the form t(x, 0) = ^. e lke e x \ k \ ^ with 
x real, which is divergent if x is positive. It is, however, implicit in Oota's formulation that any 
terms which are naively divergent must be analytically continued. For x negative, t(x, 0) is just 
a standard Fourier transform which has the result 2i arctan — -. Thus the analytic continuation 

X J 

x — ► —x to x positive should just introduce a minus sign, so each divergent term of this type 
with x positive should be replaced by the same term with x — > —x and an additional minus sign. 

If = 0, each term t(x, 0) becomes 0, unless x = 0. If there is no i(0, 0) term, the rhs must 
therefore reduce to 1. If 9 ^ 0, the limit ordering does not matter, so we can take the e limit 
first and reduce X^[=i[^]gg^i( { ?( 7r ^)' Qi^k)) to 5ij[tj\q^ k y Each t(x,9) is then matched by a 
t(—x,9), so the rhs again reduces to 1. The only way the rhs can come to anything other than 
1 overall for any 9 is if there are terms like i(0, 0) present. 

For this to happen, we require [Kn] ql gi = (q' , (f -independent part) + (terms in q',q'). From 
the definition ( 4.33| ), and the fact that [n] q can be expanded out as q n_1 + q n ~ 3 + . . . + q~ 



for n integer, this reduces to requiring Gu odd, in which case the q' , g'-independent part is -1. 
We also need Mij(q(irk),q(-7rk)) = q(-Kk)q(irky^ +(teicms in higher powers of q,q) for the same 
/. Expanding out (4.44), the lowest power of q a q b present is q x q y , for the smallest x and y 
such that mij(x,y) ^ 0. The second condition can thus only be satisfied if the product-form 
S-matrix Sij(0) contains a block {l,tj}, and flB,3j) sh ows that this occurs iff I = j. This should 
be compared with the discussion of equation (|B.12[) , where the discrepancy between the two 
possible limit prescriptions was caused by a pole from this block; we are essentially approaching 
the same pole here. 

Overall, then, we find that there is one block of the form t(0, 0) if Gij is odd, but none 
otherwise. In this case, we find y = for Gij even and 

e y = exp (- f°° Y eike ) ' G v odd - ( R22 ) 

This is just the x — > case of the previous Fourier transform, so we find y = —2iir®(0) for Gij 
odd or y = otherwise. This is equivalent to y = —2iirQ(0)Gij, showing that we have indeed 
found a generalisation of the RTV formula. 



B.4 Fourier transforms 
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To complete this section, we must discuss the exceptional Being self-dual, the 

S-matrix for this theory cannot be found from the above. Following Oota, however, we note 
that the necessary prescription is to replace each reference to r v fe v by h v = h = 2n + 1, take 
all t a = 1, and replace the incidence matrix by the "generalised incidence matrix" ||53[|, which is 
obtained from the incidence matrix of an by replacing the last zero on the the diagonal by a 
one. Doing this, we obtain the correct integral S-matrix, and hence a generalised RTV identity, 
for this case. 



B.4 Fourier transforms 

Here we attempt to find 

**>=f (=t#t4W < b - 23 > 

J-oo \smh (| + a) J 

To do this, we need to use the Convolution Theorem, which states that, if F{a) and G{a) 
are the Fourier transforms of f(x) and g(x) respectively, then 

1 /'OO poo 

— / F(a)G(a)e- iax da = / f(u)g(x - u)du. (B.24) 

^ J —oo J — oo 

This, together with the standard results 

,. , cosh(ax) n , . zvrsinh(^) 

f(x) = , : , < a < b -> F(a) = — r -. — , , — rr (B.25) 

J K ' sinh(6x) ' y ' b [cosh {?f) + cos (^)] V ' 

f(x) = 5{x) -» F{a) = 1 (B.26) 
allow us, setting a = it and b = 2tt in |B.25| to find 

— f ( .—to _ f - (B.27) 

4-7T J.^ ycosh (I) / J_ 00 smh(27rn) 



Returning now to B.24, if we make the change of variables 9' = 9 + 2a + m7r, where n is an odd 



integer, we find it becomes 



<f,{k) = e Ut(?a+tom) I I \ZS I jiurjff ( B _ 28 ) 




where we have implicitly moved the contour of integration from the real axis to a line 2a + inn 
above it. We can do this provided there are no poles of the function between the real axis and 
this line, and, when we make use of this result, we will pick the arbitrary constant n to make 
sure this happens. If we were to take n such that there were m simple poles in this region, we 
would incur a correction term of i2nm, being a contour integral of the function with the contour 
going along the real axis to infinity, up to 2a + inn, back along this line to minus infinity and 
then back down to the real axis and back to the start. 

We are now in a position to connect the above together, and find, finally 

kk) = 4nie lk ^ 2a+m ^ cosh ^ k ) (B.29) 
Yy ' sinh(27r£;) v ; 



(being careful over the sign, due to the discrepancy in the sign of the exponential between B.27] 
and [R28D . 
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B.5 Dynkin diagrams 



Where there are roots of different lengths, the filled in spots refer to short roots. 
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B.6 Cartan matrices for simple Lie algebras 

Here, we give explicitly the Cartan matrices for all the untwisted simple Lie algebras, with the 
root ordering and normalisation we have used. 



c (1) 



d (1) 



( 2 -1 • 
-1 2 -1 • 
0-12- 



2-10 
-1 2 -1 
0-12/ 



/ 2 -1 
-1 2 -1 
0-12 



2-10 
-1 2 -1 
0-22/ 



/ 2 -1 
-1 2 -1 
0-12 



2-10 
-1 2 -2 
0-12/ 



/ 2 -1 
-1 2 -1 
0-12 



2 


-1 








-1 


2 


-1 


-1 





-1 


2 








-1 





2 J 
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e (1) 



/ 2 -1 \ 

2 -1 

2 -1 -1 
0-1-12-10 

-10 0-12 

\ -1 2/ 

e (1) 
e 7 
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-1 
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-1 
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-1 
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-2 
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-1 


-1 
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B.7 S-matrices 



For the self-dual cases, the S-matrices were originally found in [41, 42, 17]. The non-self-dual 



of 1 45] and write the S-matrix as 



cases took a little longer but were finally obtained in [46, 47, 45 1. We adopt the general notation 

(B.30) 



h r v h v 

S ab (0) = U H{x,y} m ^ 
x=l y=l 



where the {x, y} are of the form 



{x,y} 



(x-l,y-l){x + l,y + l) 
{x-l,y + l)(x + l,y-l) 



with 



and 



, {2-B)x , By 

{x,y} = ( — - + 



2h 



(x) 



sinh 



2r v fr v/ ' 

+ ITTX)) 



sinh (i (6 — ittx)) 



(B.31) 
(B.32) 
(B.33) 



For convenience, this notation can be extended to include 

a [x, y] b = a {x, y} b x crossing = a {x, y} b x {h - x, r y h y - y} b 

with the subscript being omitted if it is equal to one and 

{x,yh = {x,y - l}{x,y + 1} 

(x,yh = {x,y - 2}{x,y}{x,y + 2} 

2{x,y}2 = {x -l,y} 2 {x + l,y} 2 

= {x-l,y- 1}{x -l,y+ l}{x + l,y- l}{x + l,y+ 1}. 

Whenever an entry appears to the power n below, this means that m a b(x, y) should be taken 
to be n rather than 1 for that entry. 



an h = n + 1 and r v /i v = n + 1 



Sab(0) 



a+b-1 

n 

p=\a— fe|+l step 2 



b { n ] h = 2n and r v /i v = An - 2 

S ab (9) = 

San{0) = 



a+b-1 

Yl [P,2p] 2 , a,b< 

p=\a— b\+l step 2 
a+b-1 

]{ [p, 2p] , a <n 

p=\a— b\+l step 2 
n-1 

JJ {n - p, 2n - 1 - 2p} 

p=l—n step 2 
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Cn ] h = 2n and r v ft v = 2n + 2 



a+6-1 

Sab= [J 

p=|a— 6|+1 step 2 



dfc h = 2(n- 1) and r v ft v = 2(n - 1) 



S ab (6) = 

s ab (0) = 

S ab (0) = 

Sn(n-i)(0) = 



a+b-l 

p=|a— 6|+1 step 2 
n-2+6 

1^ {P) p} ) b < n — 1, a = n — 1 or n 

p=n— 2 step 2 
2n-4+a; 

~Y\_ {p,p} j a = b = n or a = b = n — 1 

p=l step 4 
2n-4— x 

n fr'^ 

p=3 step 4 



(In the above, x = 1 for n even and x = — 1 for n odd. 



ei 1} ft = 12 and r v /i v = 12 



(In this and the subsequent sections, x listed on its own should be taken to mean {x}.) 



Block 



1,7 
[4] 
4,6,10 
[3], [5] 
2,6,8 

5,11 
[1U5] 



Block 



[3], [5] 
[2],[4],6 2 
[3]- [5] 
[4] 

1,[3],5,7 2 
[2],[4] 2 ,6 2 
[3],5 2 ,7,11 



Block 



2,6,8 
[1],[3] 2 ,[5] 3 
[2],[4] 2 ,6 2 

[3], [5] 
1,[3],5,7 2 
4,6,10 
1,7 



ei 1} ft = 18 and r v ft v = 18 



m a b(x,y) = 1 for (x,y) = (p,p) where 



a 


b 


Block 


a 


b 


Block 


a 


b 


Block 


1 


1 


W,9 


2 


5 


[2], [6], [8] 


4 


6 


[2],[4],[6],[8] 2 
[3],[5] 2 ,[7] 2 ,9 2 


1 


2 


[6] 


2 


6 


[4], [6], [8] 
[3],[5],[7],9 2 


4 


7 


1 


3 


[5], 9 


2 


7 


5 


5 


[1],[3],[5],[7] 2 ,9 2 
[3],[5] 2 ,[7] 2 ,9 2 
[2],[4] 2 ,[6] 2 ,[8] 3 


1 


4 


[2], [8] 


3 


3 


[1],[5],[7],9 


5 


6 


1 


5 


[5], [7] 


3 


4 


[4], [6], [8] 
[3],[5],[7],9 2 
[3],[5],[7] 2 ,9 
[2],[4],[6] 2 ,[8] 2 
[1],[3],[7],9 2 
[4],[6] 2 ,[8] 


5 


7 


1 


6 


[3], [7], 9 


3 


5 


6 


6 


[1],[3],[5] 2 ,[7] 2 ,9 3 
[2],[4] 2 ,[6] 3 ,[8] 3 
[1],[3] 2 ,[5] 3 ,[7] 4 ,9 4 


1 


7 


[4], [6], [8] 


3 


6 


6 


7 


2 


2 


W,[7] 


3 


7 


7 


7 


2 


3 


[4], [8] 


4 


4 








2 


4 


[5], [7] 


4 


5 
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ePh = 30 and r v h v = 30 



m a b(x,y) = 1 for (x,y) = (p,p) where 



Si 


K 
u 


1 

1 


1 
1 


1 

± 


o 
z 


1 

± 


Q 



1 
± 




1 
1 





1 
± 


D 


1 
1 


7 


1 
± 


Q 
O 


o 
L 


9 


o 

/ 


Q 

o 


2 


4 


2 


5 


2 


6 


2 


7 


2 


8 


3 


3 


3 


4 


3 


5 



Block 



Block 

[5],[7]M9],[11],[13]V5^ 
[3],[5],[7],[9] 2 ,[11] 2 ,[13]M5 2 
[4],[6] 2 ,[8] 2 ,[10] 2 ,[12] 2 ,[14] 3 
[1],[5],[7],[9],[11] 2 ,[13],15 2 
[4],[6],[8] 2 ,[10],[12] 2 ,[14] 2 
[3],[5],[7],[9] 2 ,[11] 2 ,[13] 2 ,15 2 
[3],[5],[7] 2 ,[9] 2 ,[11] 2 ,[13] 3 ,[15] 2 
[2],[4],[6] 2 ,[8] 2 ,[10] 3 ,[12] 3 ,[14] 3 
[1],[3],[5],[7],[9] 2 ,[11] 3 ,[13] 2 ,15 2 

[4],[6] 2 ,[8] 2 ,[10] 2 ,[12] 2 ,[14] 3 
[2],[4],[6] 2 ,[8] 2 ,[10] 3 ,[12] 3 ,[14] 3 
[3],[5] 2 ,[7] 3 ,[9] 3 ,[11] 3 ,[13] 4 ,[15] 4 
[1],[3],[5],[7] 2 ,[9] 2 ,[11] 3 ,[13] 3 ,15 2 
[3],[5] 2 ,[7] 2 ,[9] 3 ,[11] 3 ,[13] 3 ,15 4 
[2],[4] 2 ,[6] 2 ,[8] 3 ,[10] 4 ,[12] 4 ,[14] 4 
[1],[3],[5] 2 ,[7] 3 ,[9] 3 ,[11] 4 ,[13] 4 ,15 4 
[2],[4] 2 ,[6] 3 ,[8] 4 ,[10] 4 ,[12] 5 ,[14] 5 
[1],[3] 2 ,[5] 3 ,[7] 4 ,[9] 5 ,[11] 6 ,[13] 6 ,15 6 



[!].["] 
[7], [13] 
[2], [10], [12] 
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Epilogue 

"The Road goes ever on and on, 
Down from the door where it began. 
Now far ahead the Road has gone, 
And I must follow, if I can, 
Pursuing it with eager feet, 
Until it joins some larger way 
Where many paths and errands meet. 
And whither then? I cannot say." 

— J.R.R. Tolkien, The Hobbit 
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